R_Language 板


LINE

※ 引述《cog5566 (刃之56)》之銘言: : [問題類型]: : 程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來) : [軟體熟悉度]: : 入門(寫過其他程式,只是對語法不熟悉) : [問題敘述]: : 以 subsetBy 為例好了 (這是doBy package 的一個function) : 這個 function 可以依照組別來篩選 data frame 的資料。 : 假設我想要寫一個 function 會用到 subsetBy : 這個 function 的功能是依組別去除 2 倍 SD 之外的資料 : 問題如下: : subsetBy 的第二個參數是要輸入一個 expression (應該是吧) subsetBy(formula, subset, data = parent.frame(), select, drop=FALSE, join=TRUE, ... ) 第二個參數是 logical 不是 expression,可以看他的說明 subset logical expression indicating elements or rows to keep: missing values taken as false. 用他的 example 看一下就知道了 sbstBy> data(dietox) sbstBy> subsetBy(~Evit, Weight < mean(Weight), data=dietox) data(dietox) dietox$Weight < mean(dietox$Weight) : 如 : subsetBy(~Evit, Weight < mean(Weight), data=dietox) : 這個 expression 是篩選資料的算式,其中的變數會對應到 data frame : 的 column name : 而我想要寫一個 function 是 : sdSubsetFun = function(group, dVar, data) : data 是我要操作的 data frame : dVar 是要去掉的 2SD 的目標 column name : group 是分組的依據 (class 是 formula) : 但是我沒辦法將 dVar 傳進去裏面的 subsetBy 的第二個參數 : 程式碼如下 : ###################################################################### : library('doBy') : sdSubsetFun = function(group, dVar, data) { : - sdVar = substitute(dVar) : # 2SD threshold : - subsetThreshold = substitute( (sdVar > (mean(sdVar) - 2*sd(sdVar))) & : - (sdVar < (mean(sdVar) + 2*sd(sdVar))) ) + subsetThreshold = (dVar > (mean(dVar) - 2*sd(dVar))) & + (dVar < (mean(dVar) + 2*sd(dVar))) : dataOut = subsetBy(group, : - subset = eval(subsetThreshold, group, parent.frame()), + subset = subsetThreshold, : data = data) : } - 是刪除,+ 是新增。這樣應該就可以用了 :p : # Create Data : aa = rep(c('a','b','c'), 50) : bb = 1:150 : tData = data.frame(aa,bb) : subData = sdSubsetFun(~aa, bb, tData) : ###################################################################### : 錯誤訊息為: : > subData = sdSubsetFun(~aa, bb, tData) : Error in eval(subsetThreshold, group, parent.frame()) : : object 'subsetThreshold' not found : 也就是說,裏面的 subsetBy 都不會抓我的 subsetThreshold 的值, : 都還是會抓 subsetThreshold 自己本身的 expression。 : 要怎麼做才會讓他抓到 subsetThreshold 的值呢? --



※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 36.229.4.177 ※ 編輯: Chris7462 來自: 36.229.4.177 (06/02 10:26)







like.gif 您可能會有興趣的文章
icon.png[問題/行為] 貓晚上進房間會不會有憋尿問題
icon.pngRe: [閒聊] 選了錯誤的女孩成為魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一張
icon.png[心得] EMS高領長版毛衣.墨小樓MC1002
icon.png[分享] 丹龍隔熱紙GE55+33+22
icon.png[問題] 清洗洗衣機
icon.png[尋物] 窗台下的空間
icon.png[閒聊] 双極の女神1 木魔爵
icon.png[售車] 新竹 1997 march 1297cc 白色 四門
icon.png[討論] 能從照片感受到攝影者心情嗎
icon.png[狂賀] 賀賀賀賀 賀!島村卯月!總選舉NO.1
icon.png[難過] 羨慕白皮膚的女生
icon.png閱讀文章
icon.png[黑特]
icon.png[問題] SBK S1安裝於安全帽位置
icon.png[分享] 舊woo100絕版開箱!!
icon.pngRe: [無言] 關於小包衛生紙
icon.png[開箱] E5-2683V3 RX480Strix 快睿C1 簡單測試
icon.png[心得] 蒼の海賊龍 地獄 執行者16PT
icon.png[售車] 1999年Virage iO 1.8EXi
icon.png[心得] 挑戰33 LV10 獅子座pt solo
icon.png[閒聊] 手把手教你不被桶之新手主購教學
icon.png[分享] Civic Type R 量產版官方照無預警流出
icon.png[售車] Golf 4 2.0 銀色 自排
icon.png[出售] Graco提籃汽座(有底座)2000元誠可議
icon.png[問題] 請問補牙材質掉了還能再補嗎?(台中半年內
icon.png[問題] 44th 單曲 生寫竟然都給重複的啊啊!
icon.png[心得] 華南紅卡/icash 核卡
icon.png[問題] 拔牙矯正這樣正常嗎
icon.png[贈送] 老莫高業 初業 102年版
icon.png[情報] 三大行動支付 本季掀戰火
icon.png[寶寶] 博客來Amos水蠟筆5/1特價五折
icon.pngRe: [心得] 新鮮人一些面試分享
icon.png[心得] 蒼の海賊龍 地獄 麒麟25PT
icon.pngRe: [閒聊] (君の名は。雷慎入) 君名二創漫畫翻譯
icon.pngRe: [閒聊] OGN中場影片:失蹤人口局 (英文字幕)
icon.png[問題] 台灣大哥大4G訊號差
icon.png[出售] [全國]全新千尋侘草LED燈, 水草

請輸入看板名稱,例如:Boy-Girl站內搜尋

TOP