PTT
Submit
Submit
选择语言
正體中文
简体中文
PTT
R_Language
[问题] 用for循环算zscore
楼主:
amygm307
(简单点ˇ)
2016-07-05 16:01:18
- 问题: 用for循环算z-score遇到的问题
[问题类型]
程式咨询
[软件熟悉度]:
新手(没写过程式,R 是我的第一次)
[问题叙述]:
不好意思小妹初学R又来打扰
想要把分数依照每位受试者去算z score
每位受试者有8个观察值
于是我先分别计算平均数和标准差,如下
A20是我要标准化的变项,sub是受试者编号
wtpmean=aggregate(dat$A20,by=list(dat$sub),mean)
colnames(wtpmean)=c("sub","mean")
wtpsd=aggregate(dat$A20,by=list(dat$sub),sd)
colnames(wtpsd)=c("sub","sd")
接着我想要用循环来把分数标准化
zscore=0
for(i in 1:length(dat$A20)){
zscore=(subset(dat$A20,dat$sub==i)-wtpmean$mean[i])/wtpsd$sd[i]}
但是这里就出现问题了,我的zscore是空的,我不知道是哪里写错了
还是循环的用法不对呢?
[环境叙述]:
R version 3.3.0 (2016-05-03)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)
locale:
[1] LC_COLLATE=Chinese (Traditional)_Taiwan.950
[2] LC_CTYPE=Chinese (Traditional)_Taiwan.950
[3] LC_MONETARY=Chinese (Traditional)_Taiwan.950
[4] LC_NUMERIC=C
[5] LC_TIME=Chinese (Traditional)_Taiwan.950
attached base packages:
[1] stats graphics grDevices utils datasets methods
[7] base
loaded via a namespace (and not attached):
[1] tools_3.3.0
作者:
clansoda
(小笨)
2016-07-05 16:25:00
我不知道你是不是要这个,
http://pastebin.com/QkmVmXJ5
建议assign上面用<-这个会比较好,=有时候会出错
作者:
swedrf0112
(M)
2016-07-05 16:27:00
zscore[i]=... 要把数值放进去向量里记得要加[]
作者:
clansoda
(小笨)
2016-07-05 16:31:00
对分类数值做function,可以用tapply
作者:
andrew43
(讨厌有好心推文后删文者)
2016-07-05 16:33:00
dat$sub==i 这个判断式怪怪的。如果dat$sub就1,2,3...那就没问题zscore应该要设计成数组向量…实在不明白你的用意。
作者:
swedrf0112
(M)
2016-07-05 22:23:00
你的 wtpmean$mean[1] 是ㄧ个 vector 吧? 所以每次循环算出来是ㄧ个有 8 个元素的 vector?
作者:
andrew43
(讨厌有好心推文后删文者)
2016-07-05 22:24:00
不过我也是建议你不要靠loop,用tapply吧。
作者:
bowin
(尽其在我)
2016-07-06 21:16:00
建议可以vector operation, 运算会比loop, apply等快很多
楼主:
amygm307
(简单点ˇ)
2016-07-08 13:57:00
感谢各位!我后来用tapply解决了!
继续阅读
Re: [问题] 处理资料 效能询问
celestialgod
[问题] 自订网页关键字
Tampa
[问题] 处理资料 效能询问
spiderway
[心得] strptime 与中文
Edster
[问题] R graph with two y-axes
amber219212
[问题] 判断row的重复个数并做某个计算
bedroom0204
[问题] 在mac上plot中文出现方块的永久解决方案
JJabrams
[问题] 请教Rmpi设定哪里出了问题
dp2046
[问题] R的规划求解 optim()、constrOptim()
MrSomebody
[情报] R 语言入门课程(7/21 ~ 7/22,南科)
wanggz
Links
booklink
Contact Us: admin [ a t ] ucptt.com