在R函式中使用aov執行多個anova

我正在撰寫一個函式來對多個變數運行方差分析。我用一個預測器嘗試了該函式，它作業正常，現在我嘗試使用兩個變數，但該函式正在崩潰。

下面是帶有一個預測器“條件”的代碼

lapply(names(data10)[6:ncol(data10)], function(x) {
  aov.GROUPconint <- lme(reformulate('Condition', x), random = ~1|Subject, data=data10)
  GROUPconint.posthoc <- summary(glht(aov.GROUPconint, linfct = mcp(Condition = "Tukey")))
  GROUPconint.posthoc
}) -> result

使用兩個預測器——“條件”和“測驗”——代碼崩潰了，我不知道如何編輯上面的代碼

lapply(names(data10)[6:ncol(data10)], function(x) {
  aov.GROUPconint <- lme(reformulate(c('Condition', 'test'), x), random = ~1|Subject, data=data10)
  GROUPconint.posthoc <- summary(glht(aov.GROUPconint, linfct = mcp(Condition = "Tukey")))
  GROUPconint.posthoc
}) -> result

data10 <- structure(list(Subject = structure(c(12L, 24L, 12L, 24L, 12L, 
24L, 12L, 24L, 12L, 24L, 12L, 24L), .Label = c("P14", "P15", 
"P16", "P17", "P18", "P19", "P20", "P21", "P22", "P23", "P24", 
"P25", "P26", "P27", "P28", "P29", "P30", "P31", "P32", "P33", 
"P34", "P35", "P37", "P38"), class = "factor"), Condition = structure(c(1L, 
1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), .Label = c("CEN", 
"IPS", "CTL"), class = "factor"), Wave = c(1, 1, 1, 1, 1, 1, 
1, 1, 1, 1, 1, 1), test = structure(c(1L, 1L, 2L, 2L, 1L, 1L, 
2L, 2L, 1L, 1L, 2L, 2L), .Label = c("Block 1", "Block 10"), class = "factor"), 
    ...5 = c("P25-CondA1", "P38-CondA1", "P25-CondA1", "P38-CondA1", 
    "P25-CondB", "P38-CondB", "P25-CondB", "P38-CondB", "P25-CondC", 
    "P38-CondC", "P25-CondC", "P38-CondC"), Fp1 = c(-5.840139, 
    -24.64444, -5.91461, -10.269491, -1.063752, -10.417562, -3.057545, 
    -22.833533, -0.528596, -6.226959, -5.206908, -16.431027), 
    Fz = c(-0.160478, -4.796688, -0.572049, -2.42414, -0.081133, 
    -2.212509, -0.507796, -5.125658, -0.247309, -1.391691, -0.340251, 
    -3.275482), F3 = c(-3.237156, -7.13964, -2.356562, -4.498527, 
    -1.912161, -3.369108, -0.764519, -7.338588, -0.967395, -2.61488, 
    -2.006442, -6.064484), F7 = c(-12.291867, -23.007904, -13.424168, 
    -51.746685, -4.683314, -16.146104, 3.129978, -16.762208, 
    -4.072039, -3.997164, -9.893088, -30.463026), FT9 = c(-5.145666, 
    -19.855898, -7.159363, -79.547642, -2.737132, -24.988066, 
    -0.764389, -8.678315, -2.424962, -3.655473, -16.684891, -40.817396
    ), FC5 = c(-1.912046, -7.299207, -2.782865, -13.433165, -2.717561, 
    -6.404056, -1.020601, -5.731675, -1.449527, -2.455464, -4.816248, 
    -8.451322), FC1 = c(-0.014723, -3.250151, -0.025628, -2.072377, 
    -0.025264, -0.816653, -0.008836, -1.645703, -0.025248, -0.886561, 
    -0.019677, -1.440108), C3 = c(-1.673774, -3.907241, -1.517202, 
    -3.189047, -1.721737, -1.953359, -0.130617, -1.255402, -1.000825, 
    -2.154744, -4.642536, -2.5244), T7 = c(-6.692289, -11.995039, 
    -4.388139, -26.845076, -2.150749, -19.960795, 0.957167, -9.285198, 
    -2.610707, -3.387525, -16.078687, -13.23254), TP9 = c(-6.577069, 
    -24.083291, -5.319204, -17.99928, -2.009148, -22.327639, 
    -0.39423, -11.561658, -11.369999, -1.670421, -27.666286, 
    -8.352513), CP5 = c(-4.256723, -7.168578, -3.69391, -8.595433, 
    -3.172519, -5.686751, -1.439715, -5.684416, -2.23341, -3.819792, 
    -12.817709, -4.332906), CP1 = c(-2.424485, -2.105764, -2.437918, 
    -1.684633, -1.896294, -0.603067, -1.310922, -1.246486, -1.058182, 
    -1.723797, -4.112311, -0.39083), Pz = c(-5.230209, -6.931432, 
    -5.447097, -4.420237, -4.298466, -0.572419, -3.76145, -4.866791, 
    -1.882695, -2.773345, -7.048977, -1.297095), P3 = c(-5.717413, 
    -8.568655, -5.872956, -6.676226, -3.610336, -3.78263, -2.678968, 
    -5.256058, -1.84788, -4.615055, -10.663448, -3.350067), P7 = c(-8.747675, 
    -14.647256, -5.446711, -13.149602, -3.162763, -11.810952, 
    -0.844879, -10.721099, -3.343471, -1.414516, -21.813454, 
    -5.784583), O1 = c(-6.725148, -23.142534, -12.192043, -14.436759, 
    -5.107042, -8.746607, -4.771125, -13.834345, -3.327836, -0.988959, 
    -23.619593, -8.835368), Oz = c(-6.360107, -19.795088, -11.102469, 
    -13.143209, -4.797536, -6.422619, -4.972995, -13.851071, 
    -3.273859, -0.095299, -23.344231, -10.013479), O2 = c(-8.604531, 
    -16.010789, -9.163555, -14.997259, -5.795494, -7.69402, -4.861372, 
    -15.428395, -3.201026, 0.689055, -25.061542, -14.588914), 
    P4 = c(-5.786419, -8.570093, -6.01655, -7.250971, -4.000653, 
    -5.302266, -3.95116, -11.026836, -1.672558, 0.563733, -9.640803, 
    -5.720633), P8 = c(-10.232556, -11.447704, -6.468332, -13.059249, 
    -3.728454, -8.160377, -1.216276, -15.638084, -3.104394, 1.522973, 
    -17.690755, -17.926968), TP10 = c(-8.034667, -11.84704, -6.469727, 
    -17.795867, -3.392685, -14.240206, -0.750774, -14.076242, 
    -2.212408, -0.838559, -24.863607, -13.99878), CP6 = c(-4.075277, 
    -6.763672, -4.167653, -13.016007, -3.034482, -6.815944, -1.626736, 
    -10.962036, -2.496595, -0.236931, -7.798438, -6.962835), 
    CP2 = c(-2.405001, -2.252853, -2.180014, -2.301028, -2.242261, 
    -1.182271, -2.213638, -3.496284, -1.008576, -0.61297, -3.492367, 
    -0.51545), Cz = c(-0.468023, -0.004084, -0.638258, -0.00788, 
    -0.851801, -0.007214, -0.706636, -0.010855, -0.536705, -0.013132, 
    -1.039244, -0.007484), C4 = c(-1.666512, -2.444287, -1.888915, 
    -5.533711, -1.222455, -1.965361, -0.873634, -4.527433, -1.818675, 
    -0.937223, -2.389627, -2.275628), T8 = c(-6.815842, -8.228049, 
    -5.290419, -15.105313, -1.736178, -13.356379, -1.95265, -19.905574, 
    -3.711849, -9.398742, -10.633943, -12.196464), EOG = c(-20.232283, 
    -99.109516, -42.37872, -41.447092, 1.38567, -23.524751, 3.979485, 
    -17.365329, -6.501096, -8.57049, -27.861358, -35.333672), 
    FC6 = c(-3.332464, -3.964186, -3.596093, -15.233898, -0.609353, 
    -6.322399, -0.27383, -7.378963, -2.265338, -1.842532, -3.983366, 
    -3.932509), FC2 = c(-0.026439, -1.089279, -0.01404, -1.967481, 
    -0.021138, -0.912876, -0.016665, -1.514663, -0.016964, -0.101446, 
    -0.017254, -0.920491), F4 = c(-1.289919, -3.210157, -1.225391, 
    -8.599483, -0.13505, -4.758277, -0.16653, -6.277869, -1.030201, 
    -0.780301, -1.717511, -3.839278), F8 = c(-4.269784, -14.505706, 
    -9.1726, -44.947816, 0.269252, -16.018926, -2.184163, -13.581011, 
    -2.464344, -5.565065, -6.70289, -12.713321), Fp2 = c(-10.068992, 
    -12.382243, 0.860865, -23.860863, 0.389788, -13.042772, 0.540714, 
    -21.70098, -2.526254, 2.619264, -3.399223, -59.010281)), row.names = c(NA, 
-12L), class = c("tbl_df", "tbl", "data.frame"))

這是錯誤資訊

Error in na.fail.default(list(Fp1 = c(-4.50016, -1.485162, -2.205266, : missing values in object

uj5u.com熱心網友回復：

library(tidyverse)
library(afex)
library (lsmeans)


lapply(names(data10)[6:ncol(data10)], function(x) {
  a1 <- aov_ez("Subject", x, data10, within = c("Condition",   "test"))
  (ls1 <- lsmeans(a1, c("Condition", "test")))
  update(pairs(ls1), by=NULL, adjust = "none")

}) -> result

result

標籤：r 功能 方差分析

上一篇：創建將4個串列相互疊加的函式？

下一篇：'x'在曲線/dnorm中有什么作用，為什么它在我的函式中起作用？