小云給大家介紹一個單細(xì)胞分析的小工具：單細(xì)胞分析（http://www.biocloudservice.com/366/366.php），能夠在線輸入數(shù)據(jù)，并快速輸出結(jié)果，大家一起來使用吧。歡迎來和小云交流討論哦。

大家好，本期小云將通過一道R語言的實戰(zhàn)題目帶大家熟悉R語言的內(nèi)容，跟著小云一起來學(xué)習(xí)吧！

1.?有一個外部文件為class.txt，存儲的是某個班級學(xué)生的姓名、年齡、身 高、體重和性別的信息。

1）此文件讀入到R中并把它記為數(shù)據(jù)框class。

> class<-read.table(choose.files())

2）給數(shù)據(jù)框class的列命名為name,age,height,weight,sex

> names(class)<-c("name","age","height","weight","sex")

3）添加一個根據(jù)年齡分組的變量rank,分成兩組： low(小于等于14歲) high（大于等于15歲）

> rank<-class$age

> class<-cbind(class,rank)

> class$rank[rank<=14]<-"low"

> class$rank[rank>=15]<-"high"

4）按性別畫出體重的箱型圖（水平，凹槽，顏色控制等）

>boxplot(weight~sex,class$sex,class$weight,notch=T,col=c("red","blue"),names=c("男","女"),horizontal=T)

5）根據(jù)性別及年齡分組畫出箱型圖

>boxplot(height~sex,data=class,notch=F,col=c("red"),names=c("男","女"),boxwex=0.25,at=1:2-0.2)

>boxplot(height~sex,data=class,add=T,at=1:2+0.2,subset=rank=="low",col=c("green"),boxwex=0.25)

> legend("bottom",c("high","low"),fill=c("red","green"))

6）畫出體重與身高關(guān)系的散點圖

> plot(weight,height,col="blue")

7）通過設(shè)置參數(shù)mfrow將所有圖放到同一頁中并保存為一個png文件

> par(mfrow=c(1,3))

>boxplot(weight~sex,class$sex,class$weight,notch=T,col=c("red","blue"),names=c("男","女"),horizontal=T)

>boxplot(height~sex,data=class,notch=F,col=c("red"),names=c("男","女"),boxwex=0.25,at=1:2-0.2)

>boxplot(height~sex,data=class,add=T,at=1:2+0.2,subset=rank=="low",col=c("green"),boxwex=0.25)

> legend("bottom",c("high","low"),fill=c("red","green"))

> plot(weight,height,col="blue")

> png("boxplot.png")

> dev.off

2.

> a<-t(VADeaths)

> barplot(a,beside = T,col = c("red","blue","green","purple"),legend=rownames(a))

以上就是本期的內(nèi)容啦，跟著小云每天進(jìn)步一點點吧！

請持續(xù)關(guān)注小云，我們下期見~

篩選免疫基因?qū)π」ぞ撸簑ww.biocloudservice.com/589/589.php