기본 콘텐츠로 건너뛰기

통계관련 함수와 메서드 사전

A B C d E F G H I K L M N O P Q R S T U V W Z A statsmodels.ap.stats.anova_lm(x) statsmodels.formula.api.ols 에 의해 생성되는 모형 즉, 클래스 인스턴스(x)를 인수로 받아 anova를 실행합니다. np.argsort(x, axis=-1, kind=None) 객체 x를 정렬할 경우 각 값에 대응하는 인덱스를 반환합니다. Axis는 기준 축을 지정하기 위한 매개변수로서 정렬의 방향을 조정할 수 있음(-1은 기본값으로 마지막 축) pandas.Series.autocorr(lag=1) lag에 전달한 지연수에 따른 값들 사이의 자기상관을 계산 B scipy.stats.bernoulli(x, p) 베르누이분포에 관련된 통계량을 계산하기 위한 클래스를 생성합니다. x: 랜덤변수 p: 단일 시행에서의 확률 scipy.stats.binom(x, n, p) 이항분포에 관련된 통계량을 계산하기 위한 클래스를 생성합니다. x: 랜덤변수 n: 총 시행횟수 p: 단일 시행에서의 확률 C scipy.stats.chi2.pdf(x, df, loc=0, scale=1) 카이제곱분포의 확률밀도함수를 계산 $$f(x, k) =\frac{1}{2^{\frac{k}{2}−1}Γ(\frac{k}{2})}x^{k−1}\exp\left(−\frac{x^2}{2}\right)$$ x: 확률변수 df: 자유도 pd.concat(objs, axis=0, join=’outer’, …) 두 개이상의 객체를 결합한 새로운 객체를 반환. objs: Series, DataFrame 객체. Axis=0은 행단위 즉, 열 방향으로 결합, Axis=1은 열단위 즉, 행 방향으

R DataFrame(데이터 프레임)

내용

Data frames

데이터 프레임은 다른 열에 다른 데이터 모드(숫자, 문자 등)가 포함될 수 있다는 점에서 행렬보다 더 일반적입니다. 데이터 프레임은 R에서 다룰 가장 일반적인 데이터 구조입니다. 여러 데이터 모드가 존재하는 경우 행렬로 나타낼 수 없습니다.data.frame() 함수로 데이터 프레임이 생성됩니다.

 ID<-c(1,2,3,4)
age<-c(25, 34, 28, 52)
diabetes<-c('type1','type1','type2','type1')
stats<-c('poor','improved','excellent','poor')
data<-data.frame(ID,age, diabetes, stats) 
data
 ID age diabetes     stats
    1  1  25    type1      poor
    2  2  34    type1  improved
    3  3  28    type2 excellent
    4  4  52    type1      poor
하나의 열에는 하나의 모드로만 구성되지만 열마다 모드는 다를 수도 있습니다. 일반적으로 각 열은 각 변수에 대한 값으로 열(column)을 변수(variable)로 나타내기도 합니다.
데이터 프레임에서 요소를 호출하는 다양한 방법이 있습니다. 먼저 벡터나 행렬과 같이 인덱스를 사용합니다. 데이터프레임 역시 행렬구조이므로 각 요소의 인덱스는 행과 열을모두 포함합니다. 그러나 하나의 인덱스만을 표시할 경우 열을 나타냅니다.
data[1,2]
[1] 25
data[1,]
ID age diabetes stats
    1  1  25    type1  poor
data[,2]
[1] 25 34 28 52
data[2]
age  
    1  25
    2  34
    3  28
    4  52
data[1:2]
ID age
    1  1  25
    2  2  34
    3  3  28
    4  4  52
열이름으로 호출할 수 있습니다.
data[c('diabetes','stats')]
diabetes     stats
  1    type1      poor
  2    type1  improved
  3    type2 excellent
  4    type1      poor
'$'와 함께 열이름으로 호출할 수 있습니다.
data$ID
[1] 1 2 3 4
data$diabetes
[1] "type1" "type1" "type2" "type1"
하나의 열을 인덱스로 호출할 경우 호출된 객체 역시 dataframe 형을 유지합니다. 그러나 "$열이름"으로 호출된 객체는 요소자체의 atomic type 즉, number 또는 character 형이 됩니다. 결과적으로 atomic이 됩니다.
위 변수(열) diabetes, stats는 각각 2개, 3개의 수준으로 그룹화할 수 있습니다. 즉, Factor로 만들 수 있습니다. 이 두 변수에 대한 교차표(cross table)은 table()함수를 사용하여 작성할 수 있습니다.
diaFac<-factor(data$diabetes); diaFac
[1] type1 type1 type2 type1
      Levels: type1 type2
statfac<-factor(data$stats);statfac
[1] poor      improved  excellent poor     
      Levels: excellent improved poor
table(diaFac, statfac)
statfac
    diaFac  excellent improved poor
      type1         0        1    2
      type2         1        0    0
위의 교차표는 데이터프레임의 열을 호출하여 직접적으로 작성할 수 있습니다.
table(data$diabetes, data$stats)
excellent improved poor
type1         0        1    2
type2         1        0    0
데이터프레임에서 열을 탐색할 경우 attach(), detach(), with()를 적용하면 코드가 간단해 집니다.
 summary(mtcars$mpg)
Min. 1st Qu.  Median    Mean 3rd Qu.    Max.
  10.40   15.43   19.20   20.09   22.80   33.90
plot(mtcars$mpg, mtcars$disp)
위 코드는 attch()와 detach() 함수를 사용하여 다음과 같이 사용합니다.
attach(mtcars) 
summary(mpg)
Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  10.40   15.43   19.20   20.09   22.80   33.90
plot(disp, hp)
detach(mtcars)
with()문을 사용하여 데이터프레임의 작업들을 그룹화할 수 있습니다. 그룹화는 중괄호 내에서 이루어지며 여러개의 명령을 실행할 경우 개행을 해야 합니다.
with(mtcars, {summary(mpg, disp, wt)                  
  + plot(mpg, disp)})
R의 중괄호 내에 있는 코드의 결과는 print() 함수에 의해 출력하거나 개체에 할당하는 것으로 반환됩니다. 그러므로 위 코드의 summary() 함수의 결과는 반환되지 않습니다. 다음 코드에서 이 결과를 반환하기 위해 개체에 할당하였습니다. 개체화에 사용되는 할당연산자는 '<-'이 아닌 '<<-'을 사용해야 합니다.
with(mtcars, {stats<-summary(mpg, disp, wt)
    + plot(mpg, disp)})
    stats
에러: 객체 'stats'를 찾을 수 없습니다
with(mtcars, {stats<<-summary(mpg, disp, wt)
    + plot(mpg, hp)})
     stats
Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
      10.40   15.43   19.20   20.09   22.80   33.90
data.frame()에서 인수 row.names에 특정한 열이나 특정한 벡터를 전달하여 행 이름을 지정할 수 있습니다.
data1<-data.frame(ID,age, diabetes, stats, row.names=age); data1 
ID age diabetes     stats
     25  1  25    type1      poor
     34  2  34    type1  improved
     28  3  28    type2 excellent
     52  4  52    type1      poor

댓글

이 블로그의 인기 게시물

matplotlib의 그래프 종류

1. 산포도(scatter plot) plt.scatter(x, y) >>> import matplotlib.pyplot as plt >>> import numpy as np >>> data=np.random.rand(1024, 2) >>> data[:3, :] >>> plt.scatter(data[:,0], data[:,1]) >>> plt.show() 2. 막대그래프(bar chart) plt.bar(x, hight, width, align='center') 매개변수중 width에 인수를 전달하여 막대의 두께를 조절할 수 있다. 또한 align의 인수는 'center'와 'edge' 이다. 기본값은 'center'이다. 이 값은 x축의 레이블이 막대의 중간에 위치(center) 또는 왼쪽 가장자리에 위치(edge)시킨다. 코드에서 np.random.randint 는 특정한 범위내에서 지정한 갯수의 랜덤수를 생성 np.unique(배열, retrun_counts=False, axis=None) : 객체 내의 중복되지 않은 수들을 반환한다. return_counts=True이면 각 수에 대한 빈도수를 반환한다. axis를 통해 행(1), 열(0)을 선택한다. >>> x=np.random.randint(1, 6, size=100) >>> uni,count=np.unique(x, return_counts=True) >>> uni array([1, 2, 3, 4, 5]) >>> count array([25, 17, 23, 16, 19], dtype=int64) >>> plt.bar(uni, count) >>> plt.show() 위의 막대그래프의 막대의

유사변환과 대각화

내용 유사변환 유사행렬의 특성 대각화(Diagonalization) 유사변환(Similarity transformation) 유사변환 n×n 차원의 정방 행렬 A, B 그리고 가역 행렬 P 사이에 식 1의 관계가 성립하면 행렬 A와 B는 유사하다고 하며 이 변환을 유사 변환 (similarity transformation)이라고 합니다. $$\begin{equation}\tag{1} A = PBP^{-1} \Leftrightarrow P^{-1}AP = B \end{equation}$$ 식 1의 유사 변환은 다음과 같이 고유값을 적용하여 특성 방정식 형태로 정리할 수 있습니다. $$\begin{align} B - \lambda I &= P^{-1}AP – \lambda P^{-1}P\\ &= P^{-1}(AP – \lambda P)\\ &= P^{-1}(A - \lambda I)P \end{align}$$ 위 식의 행렬식은 다음과 같이 정리됩니다. $$\begin{align} &\begin{aligned}\textsf{det}(B - \lambda I ) & = \textsf{det}(P^{-1}(AP – \lambda P))\\ &= \textsf{det}(P^{-1}) \textsf{det}((A – \lambda I)) \textsf{det}(P)\\ &= \textsf{det}(P^{-1}) \textsf{det}(P) \textsf{det}((A – \lambda I))\\ &= \textsf{det}(A – \lambda I)\end{aligned}\\ &\begin{aligned}\because \; \textsf{det}(P^{-1}) \textsf{det}(P) &= \textsf{det}(P^{-1}P)\\ &= \t

sympy.solvers로 방정식해 구하기

sympy.solvers로 방정식해 구하기 대수 방정식을 해를 계산하기 위해 다음 함수를 사용합니다. sympy.solvers.solve(f, *symbols, **flags) f=0, 즉 동차방정식에 대해 지정한 변수의 해를 계산 f : 식 또는 함수 symbols: 식의 해를 계산하기 위한 변수, 변수가 하나인 경우는 생략가능(자동으로 인식) flags: 계산 또는 결과의 방식을 지정하기 위한 인수들 dict=True: {x:3, y:1}같이 사전형식, 기본값 = False set=True :{(x,3),(y,1)}같이 집합형식, 기본값 = False ratioal=True : 실수를 유리수로 반환, 기본값 = False positive=True: 해들 중에 양수만을 반환, 기본값 = False 예 $x^2=1$의 해를 결정합니다. solve() 함수에 적용하기 위해서는 다음과 같이 식의 한쪽이 0이 되는 형태인 동차식으로 구성되어야 합니다. $$x^2-1=0$$ import numpy as np from sympy import * x = symbols('x') solve(x**2-1, x) [-1, 1] 위 식은 계산 과정은 다음과 같습니다. $$\begin{aligned}x^2-1=0 \rightarrow (x+1)(x-1)=0 \\ x=1 \; \text{or}\; -1\end{aligned}$$ 예 $x^4=1$의 해를 결정합니다. solve() 함수의 인수 set=True를 지정하였으므로 결과는 집합(set)형으로 반환됩니다. eq=x**4-1 solve(eq, set=True) ([x], {(-1,), (-I,), (1,), (I,)}) 위의 경우 I는 복소수입니다.즉 위 결과의 과정은 다음과 같습니다. $$x^4-1=(x^2+1)(x+1)(x-1)=0 \rightarrow x=\pm \sqrt{-1}, \; \pm 1=\pm i,\; \pm1$$ 실수