作者chiushandy (是你制造了天气)
看板Statistics
标题[问题] 用R
时间Tue Jan 10 18:51:21 2012
用R做残差分析的问题
资料结构: (整笔资料叫Adap)
age glucose_consentration .... Result
21 125 1
23 83 0
25 94 0
(以下省略)
....
目标:
我想对Result分别为0和1的资料做t.test 检查其glucose concentraion是否显着差异
在这之前要先诊断资料是否符合常态假设
做法:
resultf <-factor(Adap$Result) #对result 做factor
tapply(Adap$Glucose, resultf, resid) #残差分析
tapply(Adap$Glucose, resultf, qqnorm) #画QQ plot
问题:
原以为使用tapply指令画qq plot会自动跑出两张图 一张是result为0的资料
一张是result为1的资料 结果只有一张@@
而resid函数跑不出来 log显示
"错误在object$na.action : $ operator is invalid for atomic vectors"
请问我的做法哪里错了@@
该如何解决?
先谢谢大家!
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.220.137
1F：→ gsuper:as.data.frame(resultf) 试试看 01/10 20:11
2F：→ gsuper:再来是检查看看资料有没有 NA 或 NaN 01/10 20:20
3F：→ gsuper:即使tapply内建可以处理NA , 但 qqnorm 和 resid 不见得OK 01/10 20:21
4F：→ chiushandy:感谢楼上~我後来把0和1的资料分别存进两个变数 01/10 21:04
5F：→ chiushandy:就解决了:) 01/10 21:04