我在R中有一个数据集，有两列分别标记为x和y的列，每个列的值都超过1000。我需要找到sum((xi^2-xbar^2)(yi-ybar))/sum((xi-xbar)^4)线性回归问题。我能想到的是：
sum(((data$x)^2-mean(data$x)^2)(data$y-mean(data$y)))/sum((data$x-mean(data$x))^4)
但这给了我Error: attempt to apply non-function。我不知道如何解决这个问题。任何帮助将非常感激。
问题：您如何找出视觉上令人难以置信的表达式中存在的问题？
答：将它拆散。
df <- data.frame(x = rnorm(10), y = rnorm(10))
df$x^2
# works fine
df$x^2 - mean(x)^2
# works fine **SEE NOTE**
sum(df$x^2 - mean(x)^2)
# works fine
# sum(DF$x^2 - mean(x)^2)(data$y-mean.... oh i see
您正在尝试通过将括号彼此相邻来进行乘法。使用*
注意：没有……第二遍，您可能会发现您的值不正确，但是如果您的环境中已有一个x对象（并且该对象没有没有任何NA值）
本文收集自互联网，转载请注明来源。
如有侵权，请联系 [email protected] 删除。
我来说两句