2018年4月20日 星期五

R預設的資料集ChickWeight

R語言內建的資料集ChickWeight記錄了雞隻在不同的飼料飼養環境下體重的不同。總共有578筆資料,每筆資料分成四個欄位:
  • weight:體重,單位是公克
  • Time:出生後幾天測量體重
  • Chick:每一隻雞的編號,總共有50隻編號從1到50
  • Diet:四種不同的蛋白質飼料環境,1, 2, 3, 4
藍色是影響體重的解釋變數,綠色的體重是應變數

用指令產生不同的飼料環境的成長圖:
names(ChickWeight)
l <- split(ChickWeight, ChickWeight$Diet)
names(l)
l[["1"]]

par(mfrow=c(2,2))
plot(l[["1"]]$Time, l[["1"]]$weight)
plot(l[["2"]]$Time, l[["2"]]$weight)
plot(l[["3"]]$Time, l[["3"]]$weight)
plot(l[["4"]]$Time, l[["4"]]$weight)


或是根據不同的雞隻畫出體重隨時間的成長圖:
require(graphics)
coplot(weight ~ Time | Chick, data = ChickWeight,
       type = "b", show.given = FALSE)
_EOF_

沒有留言:

張貼留言