Re: [問題] 不會傳遞 expression 到sub-function

看板R_Language作者Wush978 (拒看低質媒體)時間11年前 (2013/06/02 22:38)推噓1(1推 0噓 7→)

留言8則, 2人參與討論串6/8 (看更多)

以下是我最初推文時測的code ```r 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))) ) data$temp <- eval(subsetThreshold) dataOut = subsetBy(group, subset = temp, data = data) } # Create Data aa = rep(c('a','b','c'), 50) bb = 1:150 tData = data.frame(aa,bb) subData = sdSubsetFun(~aa, bb, tData) ``` 不過應該和Chris7462的修改類似。 --- 關於物件的是屬於symbol或character，也就是後來cog5566提到關於`rm`函數行為的問題，這也是我過去學R ，到目前覺得很不可思議和不方便的地方。 cog5566可以用以下的方式來探索函數的行為： ```r debug(rm) x <- "y" y <- 10 rm(x) ``` 這之後會進入除錯模式，你可以巨細靡遺的觀察`rm`函數的行為。你會發現有種叫作"name"型態的物件。也可以去追相關函數的說明文件，也許你就可以搞清楚了。 ps. 我最近事情有點多，就沒追下去搞懂他了。我也不太確定Hadley的R advanced裏面有沒有提到這件事情。 ps. 現在github好像掛了，晚點我再補連結。期待你搞清楚後也可以播空發篇文章教教我，先謝謝了。 ※ 引述《cog5566 (刃之56)》之銘言： : ※ 引述《Chris7462 (～烤焦麵包～)》之銘言： : : 你只看到了第二個參數，沒看到後面的 data=dietox : : subsetBy 括號裡面整個的意思是用 dietox 這組資料的 Weight 也就是 dietox$Weight : : 可以把資料 attach 上去就知道了。 : ...嗯... : 你覺得... subsetBy 的範例有用 attach() 嗎？ : 就是因為他沒有用，所以才酷的... : : > data(dietox) : : > attach(dietox) : : > typeof(Weight < mean(Weight)) : : [1] "logical" : : > typeof(dietox$Weight < mean(dietox$Weight)) : : [1] "logical" : : > Weight == dietox$Weight : : 兩個是一樣的 : 當然啊...因為你加了 attach() ... : 你把 dietox 拿到 search path 當然這個時候 Weight < mean(Weight) : 就可以被 evaluate。 : 這不是等於說，你問我有沒有看到一隻貓，我說沒有。 : 然後你放了一隻貓在我前面，再問一次有沒有看到一隻貓... : 然後說我之前錯了，這裡其實是有一隻貓的... : : 因為你開 R 的時候有讀到舊的資料或是 R 的暫存檔，沒有 attach dietox 這個 data : : 的話應該不會 Weight 這個 variable : : 請先把舊的暫存檔或舊資料清掉，不然這樣討論下去你只會一直得到很奇怪的結果 : 這很有趣... : 不過講正經的，這真的不是 attach() 的問題。 : 從頭到尾都沒有必要 attach dietox 到 Search Path : 我很懷疑你自己有沒有真的執行過沒有 attach() 的範例 : 是真的可以執行的。跟 attach() 沒有關係。 : : 請看 ?rm : : rm (..., list = character(), pos = -1, : : envir = as.environment(pos), inherits = FALSE) : : ... : : the objects to be removed, as names (unquoted) or character strings (quoted). : : 裡面放的可以是變數名稱或是字串，所以 rm('x') 跟 rm(x) 一樣 : : 看不懂。 : 那 : x <- 10 : 請問 sum(x,x) 跟 sum(10,10) 會不會一樣呢？ : : 請把值呼叫出來 : : > x # x 是變數，數值為 10 : : [1] 10 : : > y # y 是變數，數值是字元 "x" : : [1] "x" : 所以你有試過 rm(y) 嗎？ : 你覺得答案是？ : 答案是 y 會被刪掉而不是 x : 但是爲什麽這個時候 y 不會被替換成 'x'變成 rm('x') 而刪掉 x？ : 答案是這個完全決定於 rm 的內部寫法。也就是說使用者沒辦法在外部決定輸入值是 : variable name 或是 value。 : 我不知道這有沒有嚇到你啦，但是我第一次看到這個是有被嚇到啦。 : 可能我還太嫩了。 : : 看不懂。 : 這邊我想我們各有一半的責任。 : 我的責任是說錯了一點 data frame 應該沒有包含一個 environment，就只是變數 : 和值的對應關係。 : 你的責任是你小小作弊式的用了 attach() 然後再做 : typeof(Weight < mean(Weight)) : 不然實際情況這個應該會出錯，然後會很驚訝的覺得爲什麽 : subsetBy(~Evit, Weight < mean(Weight), data=dietox) 卻不會有問題 : : 一樣看不懂。 : 但是這是我的最終問題呀。可能是我表達不好，不過你看不懂是如何回答的？ : : 問題是出在這邊，因為 data = data 用的是 data 這個資料的 data$subsetThreshold : : 但是 data 下並沒有 data$subsetThreshold 這個變數，所以就有 error : 沒錯這是個問題，而且是大問題。但是你有看到我貼的 error message 嗎 : 在這個問題發生之前上面就出問題了。 : : 另外這邊也錯，兩個改法。一個是 : : > tData = data.frame(aa <- x, bb <- y) : : 或是 : : > tData = data.frame(aa = x, bb = y) : : > aa <- x : : > bb <- y : : 看不懂上面的幾個錯誤跟 dynamic scope 有什麼關係？ : : 整的問題只是呼叫 function 的時候沒有把該放的參數放對而已。 : 這真的不是這個問題 : : 我把整個改過的程式貼上來如下，如果再不行我也不知道怎麼幫你了.... : 這個有趣... : : library('doBy') : : sdSubsetFun = function(group, dVar, data) { : : # 2SD threshold : : data$subsetThreshold = (dVar > (mean(dVar) - 2*sd(dVar))) & : : (dVar < (mean(dVar) + 2*sd(dVar))) : : dataOut = subsetBy(group, : : subset = subsetThreshold, : : data = data) : : } : : x = rep(c('a','b','c'), 50) : : y = 1:150 : : tData = data.frame(aa = x, bb = y) : : aa <- x : : bb <- y : : subData = sdSubsetFun(~aa, bb, tData) : 你這樣可以跑是因為你用了跟剛剛加了 attach() 的作弊方式 : aa <- x : bb <- y : 這跟 attach(tData) 是類似的意思 : 但是這跟原本的 subsetBy那種寫法意思就差多了。 : 我要的是 bb 不是單純的 value 代進去sdSubsetFun : 而是 bb 本身的 name 會被 sdSubsetFun 在裏面識別出來，進而組成一個新的 : expression 當做 subsetBy 的引數。 : 不過如果你看不懂 : x <- 10 : rm(x) 跟 rm(10) 的這個基本又神秘的地方，那我也不知道該怎麼跟你說明了。 -- 歡迎到ptt R_Language版分享R 的相關知識歡迎加入 Taiwan R User Group : http://www.facebook.com/Tw.R.User 聚會報名 http://www.meetup.com/Taiwan-R/ 聚會影片 https://www.youtube.com/user/TWuseRGroup -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 61.230.51.57