[SAS]
DATA a1;
INPUT gender $ wei hei age;
CARDS;
F 65 171 23
F 66 172 24
F 69 176 38
M 67 173 43
M 68 177 40
M 72 178 42
;
PROC REG;MODEL wei=hei;
RUN;
[R]
gender <- c("F","F","F","M","M","M")
wei <- c( 65,66,69,67,68,72)
hei <- c(171,172,176,173,177,178)
age <- c(23,24,38,43,40,42)
health <- data.frame(gender,wei,hei,age)
lm(wei~hei,data=health)
model_lm <- lm(wei~hei,data=health)
summary(model_lm)
out_aov <-aov(wei~hei) # 분산분석표 구하기
summary(out_aov)
[Python]
import pandas as pd
gender=["F","F","F","M","M","M"]
wei = [65,66,69,67,68,72]
hei = [171,172,176,173,177,178]
age=[23,24,38,43,40,42]
health = pd.DataFrame({'gender':gender,'wei':wei,'hei':hei,'age':age})
health
from scipy import stats # 회귀분석을 위해 Scipy 패키지의 stats 모듈을 사용
model_out = stats.linregress(hei,wei)
print(model_out)
'R & SAS 300제' 카테고리의 다른 글
비모수(non-parametric) 검정 - SAS, R, SPSS 비교 (0) | 2021.11.25 |
---|---|
단계적 회귀분석(mtcars 이용)-SAS, R 이용하기 (0) | 2021.11.21 |
상관분석 - SAS, R, Python (0) | 2021.11.20 |
t검정-두 그룹간의 평균치 차이 검정 - SAS, R, Python (0) | 2021.11.18 |
행과 열의 합계 구하기 - apply계열(sapply, lapply, tapply), aggregate (0) | 2021.11.17 |