기본 콘텐츠로 건너뛰기

벡터와 행렬에 관련된 그림들

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

댓글

이 블로그의 인기 게시물

[Linear Algebra] 유사변환(Similarity transformation)

유사변환(Similarity transformation) n×n 차원의 정방 행렬 A, B 그리고 가역 행렬 P 사이에 식 1의 관계가 성립하면 행렬 A와 B는 유사행렬(similarity matrix)이 되며 행렬 A를 가역행렬 P와 B로 분해하는 것을 유사 변환(similarity transformation) 이라고 합니다. $$\tag{1} A = PBP^{-1} \Leftrightarrow P^{-1}AP = B $$ 식 2는 식 1의 양변에 B의 고유값을 고려한 것입니다. \begin{align}\tag{식 2} B - \lambda I &= P^{-1}AP – \lambda P^{-1}P\\ &= P^{-1}(AP – \lambda P)\\ &= P^{-1}(A - \lambda I)P \end{align} 식 2의 행렬식은 식 3과 같이 정리됩니다. \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)\\ &= \textsf{det}(I)\end{aligned}\end{align} 유사행렬의 특성 유사행렬인 두 정방행렬 A와 B는 'A ~ B' 와 같...

[sympy] Sympy객체의 표현을 위한 함수들

Sympy객체의 표현을 위한 함수들 General simplify(x): 식 x(sympy 객체)를 간단히 정리 합니다. import numpy as np from sympy import * x=symbols("x") a=sin(x)**2+cos(x)**2 a $\sin^{2}{\left(x \right)} + \cos^{2}{\left(x \right)}$ simplify(a) 1 simplify(b) $\frac{x^{3} + x^{2} - x - 1}{x^{2} + 2 x + 1}$ simplify(b) x - 1 c=gamma(x)/gamma(x-2) c $\frac{\Gamma\left(x\right)}{\Gamma\left(x - 2\right)}$ simplify(c) $\displaystyle \left(x - 2\right) \left(x - 1\right)$ 위의 예들 중 객체 c의 감마함수(gamma(x))는 확률분포 등 여러 부분에서 사용되는 표현식으로 다음과 같이 정의 됩니다. 감마함수는 음이 아닌 정수를 제외한 모든 수에서 정의됩니다. 식 1과 같이 자연수에서 감마함수는 factorial(!), 부동소수(양의 실수)인 경우 적분을 적용하여 계산합니다. $$\tag{식 1}\Gamma(n) =\begin{cases}(n-1)!& n:\text{자연수}\\\int^\infty_0x^{n-1}e^{-x}\,dx& n:\text{부동소수}\end{cases}$$ x=symbols('x') gamma(x).subs(x,4) $\displaystyle 6$ factorial 계산은 math.factorial() 함수를 사용할 수 있습니다. import math math.factorial(3) 6 a=gamma(x).subs(x,4.5) a.evalf(3) 11.6 simpilfy() 함수의 알고리즘은 식에서 공통사항을 찾아 정리하...

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$$ 실수...