소개글
아래 문제에 대한 R 코드 소스입니다.
-----------------------------------------------------------
1. 5개의 그래프 재현 (그래프는 레포트에 첨부)
- 히스토그램, CLT, 카이제곱 density plot, Sine 함수 등
2. "surveyH.txt" 자료 불러오기 (데이타는 레포트 마지막 페이지에 있습니다.)
(1) sal 변수에 대해 plotting
(2) sal에 대해 plotting하는데 gen(0/1)별로 색깔,점 모양 변화 주기
(한 플랏 안에 색깔 달리 표현하기)
(3) sal에 대해 plotting 하는데, 자료 위치(점 위치) 에 교육기간이 보이도록.
교육기간 ( 5종류)
(4) (1) 을 하고 sal 중앙값에 빨갛고, 굵은선 긋기
3. x~normal(10,3) 따를때, 샘플 100개 뽑아서 신뢰구간 구하고,
구한 신뢰구간을 lines를 이용하여 plotting
(이런 신뢰구간 100개를 그리고 모수의 선 (mu=10)을 표시할 것.
=> 모수를 포함하지 않는 신뢰구간만 다른색으로 표시.
목차
1. 그래프 재현 (5)
2. "surveyH.txt" 자료 불러오기
3. x~normal(10,3) 따를때, 샘플 100개 뽑아서 신뢰구간 구하고,
구한 신뢰구간을 lines를 이용하여 plotting
본문내용
#1. (2)
n2=c();n10=c();n25=c();n100=c()
for (i in 1:500){
n2[i]=mean(runif(2))
n10[i]=mean(runif(10))
n25[i]=mean(runif(25))
n100[i]=mean(runif(100))
}
plot(n2,n100,type="n", xlim=c(0,1),ylim=c(0,14), xlab="Density estimate", ylab="f(x)",main="Simulations related CLT")
lines(density(n2),col=2,lwd=2)
lines(density(n10),col=3,lwd=2)
lines(density(n25),col=4,lwd=2)
lines(density(n100),col=6,lwd=2)
legend(0.7,8,c("n=2","n=10","n=25","n=100"),col=c(2,3,4,6),lty=1,lwd=2)
###2.
surv=read.table("surveyH.txt",header=T)
surv=as.matrix(surv)
gender=surv[,4]
x=which(gender==1)
male=surv[-x,] #gender=0
female=rbind(matrix( NA,61,4),surv[x,]) #gender=1
female2=surv[x,] #gender=1
edu=surv[,2]
#(1) salary 에 대하여 plotting
plot(surv[,1],xlab="번호",ylab="salary",xlim=c(0,100),ylim=c(3500,8300))
#(2-1)남녀 salary를 수평으로 쭉 배열...
plot(male[,1],col=3,xlab="번호",ylab="salary",xlim=c(0,100),ylim=c(3500,8300),pch=2)
points(female[,1],col=6,pch=3)
legend(70,4300,c("male","female"),col=c(3,6), pch=c(2,3))
참고 자료
없음