以下是我最初推文时测的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