嘿，大家好，本期小云將通過(guò)兩道題目進(jìn)行R語(yǔ)言的實(shí)戰(zhàn)教學(xué)。跟著小云一起學(xué)起來(lái)吧！

題目1：運(yùn)用三萬(wàn)個(gè)身高數(shù)據(jù)，分別改變樣本容量為10、100、1000隨機(jī)抽取10萬(wàn)次，對(duì)抽取后的結(jié)果取均值，將抽取不同容量樣本的曲線畫(huà)在同一塊畫(huà)布上。

seed<-1:15; ?由于所給的數(shù)據(jù)有15列，可以先隨機(jī)出一列數(shù)據(jù)，從這一列數(shù)據(jù)中抽取1000、100或10個(gè)數(shù)據(jù)，即三次for循環(huán)。使用sample隨機(jī)抽取，mean()函數(shù)取均值。

result=list()

for(i in 1:100000){

??result[i]<-mean(sample(heigh_data[,c(sample(seed,1))],1000))

}

for (i in 1:100000){

??a[i,3]=result[i];

}

將以樣本容量為1000抽取的100000個(gè)數(shù)據(jù)存放到數(shù)據(jù)框a的第三列中，同理，將樣本容量為100、10的100000數(shù)據(jù)分別存放到數(shù)據(jù)框a的第二列和第一列中。

然后使用rbind函數(shù)合并數(shù)據(jù)框成如下形式：

ggplot(all, aes(x=height, ?fill=sample))+

??geom_density(alpha=0.55,bw=1,colour="black",size=0.25)+

??theme(

????text=element_text(size=15,color="black"),

????plot.title=element_text(size=15,family="myfont",face="bold.italic",hjust=.5,color="black"),#,

????legend.position=c(0.8,0.8),

????legend.background = element_blank()

??)

###繪圖即得到圖1的效果圖

題目二：畫(huà)出由1000個(gè)0、1000個(gè)1、1000個(gè)2、……1000個(gè)9醉成的數(shù)據(jù)向量的總體分布圖和隨機(jī)抽取300個(gè)數(shù)據(jù)10000次，每次取平均值的抽樣分布

a<-rep(1:9,1000)

a<-as.data.frame(a) ??#將a轉(zhuǎn)化為數(shù)據(jù)框的形式

hist(a$a) 得到分布圖

使用for循環(huán)，抽取10000次數(shù)據(jù)，for循環(huán)中的是對(duì)每300個(gè)隨機(jī)出的數(shù)據(jù)取均值，即最后得到一萬(wàn)個(gè)均值的數(shù)據(jù)。

Hist(b$V1)得到結(jié)果圖

以上就是本期的實(shí)戰(zhàn)教學(xué)了，如果覺(jué)得對(duì)你有幫助，記得給小云點(diǎn)個(gè)贊哦，更多精彩內(nèi)容就在爾云間，跟著小云每天進(jìn)步一點(diǎn)點(diǎn)哦！今天的介紹就到這里啦，請(qǐng)持續(xù)關(guān)注小云，我們下期見(jiàn)~

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

標(biāo)簽：