風中一匹狼

#比較兩品系豬種在10間豬場檢定的平均日增重有沒有差別?
lineA<- c(.70,.80,.65,.78,.62,.67,.78,.61,.76,.64)
lineB<- c(.72,.82,.67,.83,.67,.67,.80,.64,.78,.64)
lineAB<- lineA-lineB #去除各檢定豬場差異
t.test(lineAB)

#由種豬登錄資料逢雞取得品系A與B各十個豬場的檢定資料，得到其平均日增重，兩品系種豬隻平均日增重是否有差別?
lineA<- c(.64,.61,.62,.67,.65,.70,.80,.78,.76,.78)
lineB<- c(.64,.64,.67,.67,.80,.83,.82,.72,.78,.67)
var.test(lineA,lineB)
t.test(lineA,lineB)

#以外觀將豬肉分成正常肉、水漾肉與暗乾肉三類，測定此三類里脊肉血漿某固醇量各120個樣本。所測得水漾肉與暗乾肉平均值差異為-3.3γ/100ml，均方分別為S1^2=128.88與S2^2=136.31。
#水漾肉與暗乾肉之此類固醇平均含量是否有差異?
L<- -3.3 ;S1<- 128.88; S2<-136.31 ;n<- 120
S_pooled<- sqrt(((n-1)*(S1+S2))/(n*2-2))*sqrt(1/n*2)
pt(abs(L/S_pooled), (n*2-2),lower.tail = F)*2

#如果平均值差異在2γ/100ml以上，希望有0.9的機率顯著地檢定出其差異，需要多大的樣本?
#若要估計所需樣本數，我們需要知道變方。假設所得估計值(132.6)為其變方。
L<- 2
n<- ((qnorm(0.05/2,lower.tail = F)+qnorm(0.9))^2)*(((n-1)*(S1+S2))/(n*2-2))*2/L^2
n

#已知餵飼小雞生大豆因含有trpsin inhibitor而造成胰臟腫大，為了解酸處理除去抑制因子的效果，以20隻餵飼未處理的大豆與10隻餵飼酸處理大豆來比較，測得其胰臟佔體重百分比如下，酸處理是否有效?
ctrl<- c(.75,.97,.60,.82,.84,.98,.71,.78,.88,.55,.79,1.03,.91,.63,1.10,.83,.95,.74,.87,.67)
trt<- c(.32,.74,1.05,.53,.27,.40,.97,.65,.89,.68)
t.test(ctrl,trt,var.equal = T)

#產品規格要求標準偏差在0.5以下，若測量36個產品，其標準偏差估計值為4.5，是否有足夠理由認為此產品之規格符合上述標準?
df<- 36-1
pchisq((df*4.5^2)/5^2,df)

#為了解不同劑量四氯化碳驅除鼠鉤蟲的效果，將每隻大鼠灌入500個溝蟲卵，八天後分別給予0.03ml與0.06ml含四氯化碳礦物油。兩天後將之犧牲以計算鼠鉤蟲成蟲數目。以下數據為殺死鉤蟲百分比:
trt1<- c(25.8,17.6,30.0,34.4,27.4)
trt2<- c(58.6,67.6,47.6,85.2,76.8)
var.test(trt1,trt2)
t.test(trt1,trt2,var.equal = T)

#為了解reserpine對血漿serotonin的影響，測量10位健康成人在注射2.5ml reserpine以前與一小時以後血漿serotonin濃度 (μg/10^9血小板)如下:
before<- c(.61,.59,.68,.67,.70,.54,.64,.59,.64,.64)
after<- c(.56,.46,.60,.53,.77,.50,.54,.52,.72,.60)
serotonin<- data.frame(before,after)
t.test(serotonin$before,serotonin$after, paired = T)

#為測定某飼料添加劑對生長豬增重的影響(已知絕不會對增重有負面影響)，使用10對同窩公仔豬來做試驗，每一窩各一頭飼予對照日糧與添加日糧，試驗期間增重(kg)如下:
trt<- c(40,36,36,39,43,30,35,41,48,40)
ctrl<- c(37,35,34,41,39,27,35,40,48,36)
test<-t.test(trt, ctrl, alternative = "greater",var.equal = T,paired = T)
1-0.01466
#我們有多少把握認為添加劑可以增加平均體重?
t.test(trt, ctrl, alternative = "greater",var.equal = T)
#如果添加劑可能也有負面影響，我們有多少把握認為添加劑有利於平均體重?
t.test(trt, ctrl, alternative = "two.sided",var.equal = T,paired = T)

#檢定站上月份檢定了藍瑞斯母豬52頭與杜洛克24頭，其90公斤體重時背脂厚度標準偏差分別為0.24cm與0.36cm，這兩品種背脂厚度的整齊度是否有顯著差異?
n1<- 52;n2<-24
sd1<-0.24;sd2<-0.36
F<- (sd2^2/sd1^2)
F>qf(1-0.05/2,23,51) #邏輯判斷
pf(F,n1-1,n2-1)

附上R script參考