産婦人科フクロウ blog 〜PhDからプロの研究者を目指して〜人のまねをせずに、その身に応じ、武器は自分の使いやすいものでなければならぬ

基礎の発生学、細胞物理学について勉強したことを載せていきます。古武道鍛錬中。GitHub;hidem1990

#2 記述統計とグラフ

途中

library(ISwR)

data(juul)

attach(juul)

#factor関数の力

summary(juul)

juul$sex<-factor(juul$sex,labels=c("M","F"))

juul$menarche<-factor(juul$menarche,labels=c("No","Yes"))

juul$tanner<-factor(juul$tanner,labels=c("Ⅰ","Ⅱ","Ⅲ","Ⅳ","Ⅴ"))

attach(juul)

summary(juul)

 

mid.age<-c(2.5,7.5,13,16.5,17.5,19,22.5,44.5,70.5)

acc.count<-c(28,46,58,20,31,64,149,316,103)

age.acc<-rep(mid.age.acc,count)

par(mfrow=c(1,2))

plot(age.acc,type="s")

plot(age.acc,(1:n)/n,type="s")

#縦軸と横軸、indexをどちらに設定するか

 

x<-rnorm(100)

n<-length(x)

plot(sort(x),(1:n)/n,type="s")#累計分布

 

#qq関数

qqnorm(x)

p<-seq(0,1,1/n)

plot(qnorm(p),quantile(x,p))

#こっちだと外れ値はプロット出来ない