作者ToastCheng (ToastCheng哥)
看板R_Language
標題[問題] 擷取data.frame中的特定資料
時間Tue Jan 13 18:17:50 2015
[問題類型]:
程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
新手(沒寫過程式，R 是我的第一次)
[問題敘述]:
各位大大好，小弟在做coursera的作業時有一題要求把
一個給定的data.frame（叫做airquality）中某兩個向
量（Ozone和Temp）分別大於31和90時的第三個向量Solar.R
的平均值是多少？
後來我有發現可以用：
attach(data)
sub <- subset(data,Ozone > 31 & Temp > 90)
summary(sub)
的方式求出Solar.R在符合上述兩個條件之下的平均值。
但一開始只是直覺地想要用for loop暴力算，但沒成功，
想問板上是否可以用迴圈的方式算出來？以下是自己嘗試
的程式碼，附帶一提：
data[,1]是Ozone
data[,2]是Solar.R
data[,4]是Temp
[程式範例]:
add <- 0
count <- 0
for(i in 1:153){
if(data[i,1]>31 & data[i,4]>90)
add <- add + data[i,2] #把符合條件的Temp值累加
count <- count + 1
}
mean <- add/count
這裏遇到的問題是data[,1]含有NA值，所以我想把有NA值
的資料刪掉，於是先跑這個for loop：
for(i in 1:153){
if(is.na(hw2[i,1]))
hw2[i,] <- NULL
}
結果：
錯誤在`[<-.data.frame`(`*tmp*`, i, , value = NULL) :
replacement has 0 items, need 6
不知道是類型不一樣，NULL不能隨便套用還是其他原因，
總之先感謝各位大大把這篇看完QQ
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 61.224.80.95
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1421144272.A.DAB.html
1F：推 Edster: 試試!is.na(x) & x > 31; na.omit; 並不能賦值成為 NULL; 01/13 18:56
2F：推 Edster: 若你真要賦值, 則請用 hw2[i,] <- NA 01/13 19:03
3F：推 Carollax: complete.case可以把含有NA的資料刪除 01/13 19:31