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