→ andrew43: 在R中0/1和F/T是通用的。像是做邏輯式迴歸,0/1即可。 02/02 01:53
→ andrew43: 我認為你原本的資料是洽當的。要再做各別的描述只要 02/02 01:54
→ andrew43: 按二個組切開分別分析就可以了。 02/02 01:55
→ andrew43: 但如果要一鍵完成,大概還是要自己寫。 02/02 01:56
→ andrew43: 我想到的也是用for loop走每一欄,然後用tapply進行。 02/02 02:06
感謝,請問可以直接寫函式來loop,而不用tapply嗎? 還是這樣比較好寫?
※ 編輯: EricTao (140.112.137.56), 02/02/2015 14:12:26
看來 by是傳入一個個list而非list(list()),所以只有單層迴圈
回傳值也是一樣,非單一個資料結構?
用rbind串成frame初步完成,只是as.integer還是會顯示浮點數@@"
df[names(x[i])]<-as.integer(sum(x[i]==1))
※ 編輯: EricTao (140.112.137.56), 02/02/2015 17:02:35
→ andrew43: 是否請你寫出一個雛形,看看你卡在哪裡,再來一起解? 02/03 00:16
→ celestialgod: 給partital data 給你寫的程式碼 再來討論比較好 02/03 13:14
→ celestialgod: 資料不同 怎樣整理比較快也不同 02/03 13:15
抱歉最近在看別的部分,上一次編輯其實是想表達初步解決了的意思。
承蒙如此熱情,就再麻煩大大們,看能不能改進了
資料長這樣 http://ppt.cc/i77w
目前table長這樣 http://ppt.cc/qonZ 有點醜XD 有空再看怎麼調校外觀..
--
output$DFcc<-renderTable({
if(length(input$y)<1)return()
d<-by(cova,eval(parse(text=paste("cova$",input$y,sep=""))),mtab)
#把input$y的值串在cova$上,例如cova$HTN,作為BY的分組欄位。
#用 d<-by(cova,cova[[input$y]],mtab) 好像就可以了?
f<-t(do.call(rbind,d))
#把BY的回傳值rbind再轉秩成我要的table
})
mtab<-function (x){
df<-list()
df["n"]<-nrow(x)
for(i in 1:length(x))
{
if(names(x[i])=="ID") {}#noop
#這邊想要判斷欄的變數型態不過還不知道怎麼寫好
#字串就跳過但是M/F要處理,不知道能不能用boolean?
#數字如果是1/0視為boolean,但年齡應該要show平均值&range
else if(names(x[i])=="GENDER") {
df["Male"]<-as.character(sum(x[i]=="M"))
df["Female"]<-as.character(sum(x[i]=="F"))
}
......
else {
df[names(x[i])]<-as.character(sum(x[i]==1))
}
}
return(df)
}
※ 編輯: EricTao (140.112.137.56), 02/06/2015 19:06:10
→ andrew43: 用欄名(即變數名)來判斷會不會方便一些? 02/07 17:06
→ andrew43: 如果純靠資料內容來判斷,似乎比較困難。 02/07 17:07
→ andrew43: 不過這都視你未來會丟進什麼資料而定。 02/07 17:08
→ EricTao: 對阿 目前是靠欄位名判斷 想順便練習看看就是XD 02/07 22:22
→ andrew43: length(table(x)) 如果等於 2 可以猜測為binary。 02/07 23:14
如果能判斷此欄的數值range,對了,像圖一中age的range應該用什麼方法來算呢?
同樣的方法能判斷1/0的特徵吧?
用max/min好像可以耶? 之前一直腦筋轉不過來
※ 編輯: EricTao (140.112.137.56), 02/09/2015 15:56:16