我在对我的数据运行弗里德曼测试时遇到问题。我正在尝试使用以下命令运行 Friedman 测试：
friedman.test(mean ~ isi | expId, data=monoSum)
在以下数据库 ( https://www.dropbox.com/s/2ox0y1b4gwld0ai/monoSum.csv ) 上：
> monoSum
expId isi N mean
1 m80B1 1 10 100.000000
2 m80B1 2 10 73.999819
3 m80B1 3 10 45.219362
4 m80B1 4 10 116.566174
. . . . .
18 m80L2 2 10 82.945491
19 m80L2 3 10 57.675480
20 m80L2 4 10 207.169277
. . . . . .
25 m80M2 1 10 100.000000
26 m80M2 2 10 49.752687
27 m80M2 3 10 …我对用于计算R中的滚动平均值的所有软件包都比较陌生，希望您能向我展示正确的方向。
我以以下数据为例：
ms <- c(300, 300, 300, 301, 303, 305, 305, 306, 308, 310, 310, 311, 312,
314, 315, 315, 316, 316, 316, 317, 318, 320, 320, 321, 322, 324,
328, 329, 330, 330, 330, 332, 332, 334, 334, 335, 335, 336, 336,
337, 338, 338, 338, 340, 340, 341, 342, 342, 342, 342)
correct <- c(1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0,
1, 1, 0, 0, …我想使用一个分组变量对 data.frame 中的某些数值变量运行 KW 测试。我更喜欢在循环中执行此操作，而不是输入所有测试，因为它们有很多变量（比下面的示例中的变量更多）。
模拟数据：
library(dplyr)
set.seed(123)
Data <- tbl_df(
data.frame(
muttype = as.factor(rep(c("missense", "frameshift", "nonsense"), each = 80)),
ados.tsc = runif(240, 0, 10),
ados.sa = runif(240, 0, 10),
ados.rrb = runif(240, 0, 10))
) %>%
group_by(muttype)
ados.sim <- as.data.frame(Data)
下面的代码在循环之外工作得很好。
kruskal.test(formula (paste((colnames(ados.sim)[2]), "~ muttype")), data =
ados.sim)
但它不在循环内：
for(i in names(ados.sim[,2:4])){
ados.mtp <- kruskal.test(formula (paste((colnames(ados.sim)[i]), "~ muttype")),
data = ados.sim)
}
我收到错误：
terms.formula(formula, data = data) 中的错误：模型公式中的项无效
有谁知道如何解决这个问题吗？
我有以下代码：
ggplot(iris, aes(x=Species, y=Sepal.Length)) +
stat_summary(fun.y=mean, geom='point', size=2, fill='white')
我知道可能还有其他方法可以使用iris数据绘制此平均值。但是，对于我自己的数据，这是唯一的方法。
问题：上面的代码没有给出白色填充点，而是纯黑色点。使用stat_summary参数时是否有设置填充颜色的方法？
谢谢！