那些遇到过的问题

如何在一系列软件图上打印成对比较条和效果大小值？

126*_*7b9 5 loops r boxplot pairwise

我正在努力创建一系列高质量的 ggboxplots，如下所示：

在此输入图像描述

  1. 带有带有 F(df) 测试值、p 值和效应大小的方差分析标签
  2. 通过多对比较，条形图和星形图具有显着差异。

对于上面的示例，事后比较的统计数据已通过您可以在此链接页面找到的方式获得，并且我运行了以下代码

#Compute the post-hocs
postHocs <- df %>%
  tidyr::pivot_longer(., -c(A, C, D),'s')%>%
  mutate(s = fct_relevel(s, 
                               c("E", "F", "G", 
                                 "H", "I", "J", 
                                 "K", "L", "M", 
                                 "N", "O", "P")) %>% 
  arrange(s) %>%
  group_by(s) %>% 
  pairwise_t_test(
    value ~ D, paired = TRUE,
    p.adjust.method = "bonferroni"
  ) %>% 
  #dplyr::select(., -'s')%>% 
  print()
Run Code Online (Sandbox Code Playgroud)

同时得到方差分析统计：

    res.aov <- df %>% 
      tidyr::pivot_longer(., -c(A, C, D),'s')%>%
      mutate(s = fct_relevel(s,c("E", "F", "G", 
                                 "H", "I", "J", 
                                 "K", "L", "M", 
                                 "N", "O", "P") 
                                   )))%>% 
      arrange(s) %>%
      group_by(s) %>% 
      anova_test(., value ~ D, dv = value, wid = A, within = D)%>% 
      print()
Run Code Online (Sandbox Code Playgroud)

我尝试使用以下代码获取感兴趣的统计信息：

unnested_aov <- lapply(res.aov$anova, `[[`, 1)
Run Code Online (Sandbox Code Playgroud)

我编写了这个 for 循环脚本以将统计信息添加到图表中

comparisons <- postHocs %>% add_xy_position(x = "D")
comparisons

for (i in 5:ncol(df)) {
  bxp <- ggboxplot(df,
                   x = 'D', y = colnames(df[i]))
  print(
    bxp +   stat_pvalue_manual(comparisons[,i]) +
      labs(
        subtitle = get_test_label(unnested_aov[[i]], detailed = TRUE),
        caption = get_pwc_label(comparisons[,i])))
  Sys.sleep(1)
}
Run Code Online (Sandbox Code Playgroud)

由于我返回了一些错误，我认为问题是“比较”包含 36 行，并且长度与我应该获得的图表数量（12）不太相符。

我将向您提供下面的代码，并且我将接受您在这方面的最佳建议。

谢谢

structure(list(A = 1:20, C = c("Maybe", "Maybe", "Maybe", "Maybe", 
"Maybe", "Maybe", "Maybe", "Maybe", "Maybe", "Maybe", "Maybe", 
"Maybe", "Maybe", "Maybe", "Maybe", "Maybe", "Maybe", "Maybe", 
"Maybe", "Maybe"), D = structure(c(1L, 2L, 3L, 1L, 2L, 3L, 1L, 
2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L, 3L, 1L, 2L), .Label = c("new_value_for_8", 
"new_value_for_6", "new_value_for_4"), class = "factor"), E = c(988.368784828308, 
988.856158671407, 996.004085290553, 999.685844324618, 1000.23888564896, 
1005.03749946898, 999.786378084971, 997.039675082569, 998.028313183065, 
997.168905747014, 1001.09291198164, 993.307008354785, 1004.23849942428, 
1002.98988896299, 1003.55106999008, 1009.57481668809, 1005.41677956183, 
1001.70676077155, 993.869639239065, 997.170442654021), F = c(994.834756009939, 
994.468875098246, 1000.62150212342, 1002.23100741241, 1003.96990710863, 
1007.75899775608, 998.699806256246, 996.401009591011, 998.076594704249, 
1002.19344184533, 1005.87900720863, 994.076210622421, 1002.44958531768, 
1003.10043497883, 1001.65779442628, 1013.71182331817, 1006.86178446511, 
1005.31481098188, 995.867593313022, 1000.16218228559), G = c(1011.88022669726, 
1012.10534266625, 1012.9554415821, 1015.09810043606, 1015.40462298842, 
1016.67103699915, 1003.13771453335, 999.9107434841, 1002.15365554737, 
1013.67789244066, 1014.38627383064, 1006.86762877786, 1007.47946451329, 
1008.93405130319, 1008.45962311068, 1023.4166601996, 1015.18680921429, 
1009.97183712754, 1006.2675210718, 1010.14219845841), H = c(988.221495702721, 
990.850727928741, 992.418094914622, 995.984841639886, 993.398346143465, 
997.971380355398, 1004.4672957051, 1002.54036572775, 1002.2292388993, 
999.116379988893, 997.364309124077, 997.937032776913, 1001.14544537612, 
1002.08056674659, 1000.0422658299, 1013.29862597967, 1005.06669915366, 
1003.93467692475, 1000.02290694207, 1004.31923128858), I = c(994.035709684742, 
994.890815628412, 997.18267770374, 998.564426335124, 996.851278420874, 
1000.16039368502, 1003.52155765272, 1002.1043798945, 1002.7069399281, 
1005.49897156208, 1005.81171180245, 998.62698748611, 999.56563615154, 
1002.87987510596, 998.728473297166, 1017.2093269366, 1007.79412746756, 
1008.11964589961, 1004.9525336386, 1009.50695673265), J = c(1008.23981597718, 
1009.51261484649, 1009.42367409926, 1005.06332653216, 1005.02619159395, 
1009.07903916629, 1007.56089165218, 1005.49719893791, 1004.91476855238, 
1013.03209535721, 1010.84145164945, 1005.86927622259, 1003.25309970443, 
1004.68478802971, 1002.71096740085, 1025.56743956652, 1016.32418136177, 
1013.09901927997, 1011.92002817369, 1014.69013052771), K = c(994.327042030287, 
995.608170991922, 997.033470393412, 1000.15918365269, 998.216388150646, 
1001.97377908784, 1003.17401220482, 1001.60211665164, 1002.27932356239, 
1002.41479226363, 999.832076213262, 1001.37236796086, 1001.17012593697, 
1001.40362599894, 999.964771265342, 1012.75282463779, 1008.65746780516, 
1005.290878105, 999.464067607865, 1005.14963479715), L = c(999.225538268699, 
999.349990537239, 1001.14010250645, 1001.51403741206, 1000.25571835554, 
1003.76051565494, 1002.74763442988, 1001.09116707486, 1003.29833843754, 
1006.55857216695, 1007.06029312947, 1000.60539548502, 999.637387760292, 
1002.72729847885, 998.034039799405, 1016.5065564384, 1009.68783611392, 
1009.47863905986, 1003.56318544047, 1009.23934223585), M = c(1009.99385579756, 
1011.12126521731, 1010.6989716872, 1003.7899021821, 1004.59413830322, 
1008.52123662618, 1006.34418311104, 1004.1077131243, 1004.94124365003, 
1011.89121961563, 1010.13326381032, 1005.33467168056, 1001.04545904874, 
1002.86650202467, 1000.45601490752, 1022.13789464831, 1016.25544969107, 
1012.37379951646, 1008.53744587416, 1012.00856171947), N = c(999.801263745036, 
996.838989582336, 1000.89599227983, 1003.11042068113, 1002.27800090558, 
1003.83846437952, 1000.70169995102, 1001.75290674649, 998.660833714301, 
1006.69246804854, 1004.7636391085, 1005.63873342951, 1001.37744267414, 
1000.97339668679, 1001.98775658049, 1004.70492544978, 1012.11738595707, 
1001.0458886613, 996.725751886115, 1003.17906097432), O = c(1002.96437294923, 
997.870867692911, 1002.94619035116, 1003.44844607015, 1003.02403433836, 
1004.70457675466, 999.880559826981, 1000.66826545719, 999.59436981446, 
1007.32640154038, 1006.7506344557, 1001.59973104217, 1000.71689406196, 
1001.15587576193, 999.988638552344, 1006.4489695839, 1010.51785193511, 
1003.79329103591, 999.118472788132, 1004.99936090838), P = c(1006.28027312932, 
1005.24535230967, 1007.68162285336, 1001.08242973466, 1002.99896314, 
1005.36085942954, 1001.22060069797, 1000.43007709819, 1000.47666761108, 
1008.73650967215, 1007.20593389744, 1004.57722295264, 998.66379615346, 
998.711140983915, 999.452420534917, 1008.11715753014, 1013.30601537204, 
1002.03237948844, 1002.88799699943, 1005.57921718108)), row.names = c(NA, 
20L), class = "data.frame")
Run Code Online (Sandbox Code Playgroud)

Mar*_*łka 6

这需要更多的工作。
\n让我们从数据准备开始

\n
#Loading libraries\nlibrary(tidyverse)\nlibrary(rstatix)\nlibrary(ggpubr)\nlibrary(readxl)\n\n#Upload data\ndf_join <- read_excel("df_join.xlsx")\n\ndf = df_join  %>%\n  mutate_at(vars(ID:COND), factor) %>%\n  pivot_longer(P3FCz:LPP2Pz, names_to = "signals") %>%\n  group_by(signals) %>%\n  nest()\n\n#Let\'s see what we have here\ndf\ndf$data[[1]]\n
Run Code Online (Sandbox Code Playgroud)\n

输出

\n
# A tibble: 12 x 2\n# Groups:   signals [12]\n   signals     data             \n   <chr>       <list>           \n 1 P3FCz       <tibble [75 x 5]>\n 2 P3Cz        <tibble [75 x 5]>\n 3 P3Pz        <tibble [75 x 5]>\n 4 LPPearlyFCz <tibble [75 x 5]>\n 5 LPPearlyCz  <tibble [75 x 5]>\n 6 LPPearlyPz  <tibble [75 x 5]>\n 7 LPP1FCz     <tibble [75 x 5]>\n 8 LPP1Cz      <tibble [75 x 5]>\n 9 LPP1Pz      <tibble [75 x 5]>\n10 LPP2FCz     <tibble [75 x 5]>\n11 LPP2Cz      <tibble [75 x 5]>\n12 LPP2Pz      <tibble [75 x 5]>\n> df$data[[1]]\n# A tibble: 75 x 5\n   ID    GR    SES   COND      value\n   <fct> <fct> <fct> <fct>     <dbl>\n 1 01    RP    V     NEG-CTR -11.6  \n 2 01    RP    V     NEG-NOC -11.1  \n 3 01    RP    V     NEU-NOC  -4.00 \n 4 04    RP    V     NEG-CTR  -0.314\n 5 04    RP    V     NEG-NOC   0.239\n 6 04    RP    V     NEU-NOC   5.04 \n 7 06    RP    V     NEG-CTR  -0.214\n 8 06    RP    V     NEG-NOC  -2.96 \n 9 06    RP    V     NEU-NOC  -1.97 \n10 07    RP    V     NEG-CTR  -2.83 \n# ... with 65 more rows\n
Run Code Online (Sandbox Code Playgroud)\n

现在我们必须准备将返回一些统计数据的函数

\n
#Preparation of functions for statistics\nShapiroTest = function(d, alpha=0.05){\n  st = list(statistic = NA, p.value = NA, test = FALSE)\n  if(length(d)>5000) d=sample(d, 5000)\n  if(length(d)>=3 & length(d)<=5000){\n    sw = shapiro.test(d)\n    st$statistic = sw$statistic\n    st$p.value = sw$p.value\n    st$test = sw$p.value>alpha\n  }\n  return(st)\n}\n\nsum_stats = function(data, group, value, alpha=0.05)data %>%\n  group_by(!!enquo(group)) %>%\n  summarise(\n    n = n(),\n    q1 = quantile(!!enquo(value),1/4,8),\n    min = min(!!enquo(value)),\n    mean = mean(!!enquo(value)),\n    median = median(!!enquo(value)),\n    q3 = quantile(!!enquo(value),3/4,8),\n    max = max(!!enquo(value)),\n    sd = sd(!!enquo(value)),\n    kurtosis = e1071::kurtosis(!!enquo(value)),\n    skewness = e1071::skewness(!!enquo(value)),\n    SW.stat = ShapiroTest(!!enquo(value), alpha)$statistic,\n    SW.p = ShapiroTest(!!enquo(value), alpha)$p.value,\n    SW.test = ShapiroTest(!!enquo(value), alpha)$test,\n    nout = length(boxplot.stats(!!enquo(value))$out)\n  )\n\n#Using the sum stats function\ndf$data[[1]] %>% sum_stats(COND, value)\n\ndf %>% group_by(signals) %>%\n  mutate(stats = map(data, ~sum_stats(.x, COND, value))) %>%\n  unnest(stats)\n
Run Code Online (Sandbox Code Playgroud)\n

输出

\n
# A tibble: 36 x 17\n# Groups:   signals [12]\n   signals     data        COND      n    q1     min   mean median     q3   max    sd kurtosis skewness SW.stat  SW.p SW.test  nout\n   <chr>       <list>      <fct> <int> <dbl>   <dbl>  <dbl>  <dbl>  <dbl> <dbl> <dbl>    <dbl>    <dbl>   <dbl> <dbl> <lgl>   <int>\n 1 P3FCz       <tibble [7~ NEG-~    25 -5.73 -11.6   -1.59  -1.83   1.86   9.57  4.93   -0.612   0.128    0.984 0.951 TRUE        0\n 2 P3FCz       <tibble [7~ NEG-~    25 -4.14 -11.1   -1.39  -0.522  0.659  8.16  4.44   -0.169  -0.211    0.978 0.842 TRUE        1\n 3 P3FCz       <tibble [7~ NEU-~    25 -5.35  -6.97  -1.48  -1.97   1.89   6.24  4.00   -1.28    0.189    0.943 0.174 TRUE        0\n 4 P3Cz        <tibble [7~ NEG-~    25 -2.76  -8.16   0.313  0.906  2.45  13.7   4.58    0.867   0.622    0.952 0.277 TRUE        1\n 5 P3Cz        <tibble [7~ NEG-~    25 -2.22  -9.23   0.739  0.545  3.10  14.5   4.78    1.05    0.527    0.963 0.477 TRUE        1\n 6 P3Cz        <tibble [7~ NEU-~    25 -4.14  -6.07  -0.107  0.622  2.16   7.76  3.88   -1.09    0.0112   0.957 0.359 TRUE        0\n 7 P3Pz        <tibble [7~ NEG-~    25  5.66  -0.856  8.79   7.48  11.9   23.4   5.22    0.502   0.673    0.961 0.436 TRUE        1\n 8 P3Pz        <tibble [7~ NEG-~    25  3.86  -0.888  8.45   7.88  11.2   20.7   5.23   -0.530   0.206    0.982 0.924 TRUE        0\n 9 P3Pz        <tibble [7~ NEU-~    25  2.30  -0.932  6.43   6.82   8.46  16.7   4.60   -0.740   0.373    0.968 0.588 TRUE        0\n10 LPPearlyFCz <tibble [7~ NEG-~    25 -4.02 -11.8   -0.666  0.149  1.90  13.3   5.02    0.916   0.202    0.947 0.215 TRUE        1\n# ... with 26 more rows\n
Run Code Online (Sandbox Code Playgroud)\n

接下来我们需要分析QQ图表。为此，我们准备适当的函数。

\n
#function that creates a QQ plot\nSignalQQPlot = function(data, signal, autor = "G. Anonim") data %>%\n  group_by(COND) %>%\n  mutate(label = paste("\\nSW p-value:" ,\n                       signif(shapiro.test(value)$p.value, 3))) %>%\n  ggqqplot("value", facet.by = "COND") %>%\n  ggpar(title = paste("QQ plots for", signal, "signal"),\n        caption = autor)+\n  geom_label(aes(x=0, y=+Inf, label=label))\n\n#QQ plot for the P3FCz signal\ndf$data[[1]] %>% SignalQQPlot("P3FCz")\n\n#A function that creates a QQ plot and adds it to a data frame\nAddSignalQQPlot = function(df, signal, printPlot=TRUE) {\n  plot = SignalQQPlot(df$data[[1]], signal)\n  if(printPlot) print(plot)\n  df %>% mutate(qqplot = list(plot))\n}\n\n#Added QQ charts\ndf %>% group_by(signals) %>%\n  group_modify(~AddSignalQQPlot(.x, .y))\n
Run Code Online (Sandbox Code Playgroud)\n

更正作者姓名。\n在此输入图像描述

\n

现在我们将添加 ANOVA 和 tTest 的统计数据

\n
#Add ANOVA test value\nfAddANOVA = function(data) data %>% \n  anova_test(value~COND) %>% as_tibble()\n\n#Method 1\ndf %>% group_by(signals) %>% \n  group_modify(~fAddANOVA(.x$data[[1]]))\n\n#Method 2\ndf %>% group_by(signals) %>% \n  mutate(ANOVA = map(data, ~fAddANOVA(.x))) %>% \n  unnest(ANOVA)\n
Run Code Online (Sandbox Code Playgroud)\n

输出

\n
# A tibble: 12 x 9\n# Groups:   signals [12]\n   signals     data              Effect   DFn   DFd     F     p `p<.05`      ges\n   <chr>       <list>            <chr>  <dbl> <dbl> <dbl> <dbl> <chr>      <dbl>\n 1 P3FCz       <tibble [75 x 5]> COND       2    72 0.012 0.988 ""      0.000338\n 2 P3Cz        <tibble [75 x 5]> COND       2    72 0.228 0.797 ""      0.006   \n 3 P3Pz        <tibble [75 x 5]> COND       2    72 1.61  0.207 ""      0.043   \n 4 LPPearlyFCz <tibble [75 x 5]> COND       2    72 0.885 0.417 ""      0.024   \n 5 LPPearlyCz  <tibble [75 x 5]> COND       2    72 2.65  0.077 ""      0.069   \n 6 LPPearlyPz  <tibble [75 x 5]> COND       2    72 4.62  0.013 "*"     0.114   \n 7 LPP1FCz     <tibble [75 x 5]> COND       2    72 1.08  0.344 ""      0.029   \n 8 LPP1Cz      <tibble [75 x 5]> COND       2    72 2.45  0.094 ""      0.064   \n 9 LPP1Pz      <tibble [75 x 5]> COND       2    72 3.97  0.023 "*"     0.099   \n10 LPP2FCz     <tibble [75 x 5]> COND       2    72 0.103 0.903 ""      0.003   \n11 LPP2Cz      <tibble [75 x 5]> COND       2    72 0.446 0.642 ""      0.012   \n12 LPP2Pz      <tibble [75 x 5]> COND       2    72 1.17  0.316 ""      0.031   \n
Run Code Online (Sandbox Code Playgroud)\n

对于 t 检验

\n
#Add pairwise_t_test\nfAddtTest = function(data) data %>% \n  pairwise_t_test(value~COND, paired = TRUE,\n                  p.adjust.method = "bonferroni")\n\n#Method 1\ndf %>% group_by(signals) %>% \n  group_modify(~fAddtTest(.x$data[[1]]))\n\n#Method 2\ndf %>% group_by(signals) %>% \n  mutate(tTest = map(data, ~fAddtTest(.x))) %>% \n  unnest(tTest)\n
Run Code Online (Sandbox Code Playgroud)\n

输出

\n
# A tibble: 36 x 12\n# Groups:   signals [12]\n   signals     data              .y.   group1  group2     n1    n2 statistic    df     p p.adj p.adj.signif\n   <chr>       <list>            <chr> <chr>   <chr>   <int> <int>     <dbl> <dbl> <dbl> <dbl> <chr>       \n 1 P3FCz       <tibble [75 x 5]> value NEG-CTR NEG-NOC    25    25   -0.322     24 0.75  1     ns          \n 2 P3FCz       <tibble [75 x 5]> value NEG-CTR NEU-NOC    25    25   -0.178     24 0.86  1     ns          \n 3 P3FCz       <tibble [75 x 5]> value NEG-NOC NEU-NOC    25    25    0.112     24 0.911 1     ns          \n 4 P3Cz        <tibble [75 x 5]> value NEG-CTR NEG-NOC    25    25   -0.624     24 0.539 1     ns          \n 5 P3Cz        <tibble [75 x 5]> value NEG-CTR NEU-NOC    25    25    0.592     24 0.559 1     ns          \n 6 P3Cz        <tibble [75 x 5]> value NEG-NOC NEU-NOC    25    25    0.925     24 0.364 1     ns          \n 7 P3Pz        <tibble [75 x 5]> value NEG-CTR NEG-NOC    25    25    0.440     24 0.664 1     ns          \n 8 P3Pz        <tibble [75 x 5]> value NEG-CTR NEU-NOC    25    25    3.11      24 0.005 0.014 *           \n 9 P3Pz        <tibble [75 x 5]> value NEG-NOC NEU-NOC    25    25    2.43      24 0.023 0.069 ns          \n10 LPPearlyFCz <tibble [75 x 5]> value NEG-CTR NEG-NOC    25    25   -0.0919    24 0.928 1     ns          \n# ... with 26 more rows\n
Run Code Online (Sandbox Code Playgroud)\n

现在是创建箱线图的时候了

\n
#Function to special boxplot\nSpecBoxplot = function(data, autor = "G. Anonim"){\nbox = data %>% ggboxplot(x = "COND", y = "value", add = "point")\npwc = data %>% \n  pairwise_t_test(value~COND, paired = TRUE,\n                  p.adjust.method = "bonferroni") %>% \n  add_xy_position(x = "COND")\nres.aov = data %>% anova_test(value~COND)\ndata %>% \n  ggboxplot(x = "COND", y = "value", add = "point") + \n  stat_pvalue_manual(pwc) +\n  labs(\n    title = get_test_label(res.aov, detailed = TRUE),\n    subtitle = get_pwc_label(pwc),\n    caption = autor\n  )\n}  \n\n#special boxplot for the P3FCz signal\ndf$data[[1]] %>% SpecBoxplot()\n \n#A function that creates a special boxplot and adds it to a data frame\nAddSignalBoxplot = function(df, signal, printPlot=TRUE) {\n  plot = SpecBoxplot(df$data[[1]], signal)\n  if(printPlot) print(plot)\n  df %>% mutate(boxplot = list(plot))\n}\n\n#Added special boxplot\ndf %>% group_by(signals) %>%\n  group_modify(~AddSignalBoxplot(.x, .y))\n
Run Code Online (Sandbox Code Playgroud)\n

在此输入图像描述

\n

由于修改后的函数与正态分布不一致，我使用Wilcoxon检验进行比较。我还稍微改进了箱线图。

\n
#Function to special boxplot2\nSpecBoxplot2 = function(data, signal, autor = "G. Anonim"){\n  pwc = data %>% pairwise_wilcox_test(value~COND) %>%\n    add_xy_position(x = "COND") %>% \n    mutate(COND="NEG-CTR",\n           lab = paste(p, " - ", p.adj.signif))\n  \n  data %>% ggplot(aes(COND, value, fill=COND))+\n    geom_violin(alpha=0.2)+\n    geom_boxplot(outlier.shape = 23,\n                 outlier.size = 3,\n                 alpha=0.6)+\n    geom_jitter(shape=21, width =0.1)+\n    stat_pvalue_manual(pwc, step.increase=0.05, label = "lab")+\n    labs(title = signal,\n         subtitle = get_pwc_label(pwc),\n         caption = autor)\n}\n\n#special boxplot for the P3FCz signal\ndf$data[[1]] %>% SpecBoxplot2("P3FCz")\n\n#A function that creates a special boxplot2 and adds it to a data frame\nAddSignalBoxplot2 = function(df, signal, printPlot=TRUE) {\n  plot = SpecBoxplot2(df$data[[1]], signal)\n  if(printPlot) print(plot)\n  df %>% mutate(boxplot = list(plot))\n}\n\n#Added special boxplot2\ndf %>% group_by(signals) %>%\n  group_modify(~AddSignalBoxplot2(.x, .y))\n
Run Code Online (Sandbox Code Playgroud)\n

在此输入图像描述

\n

来自我朋友的特别更新

\n

嘿@ma\xc5\x82y_statystyczny！你会很棒的！

\n

这是为您准备的特别更新。下面是一个使用参数和非参数统计数据准备图表的函数。

\n
#Function to special boxplot3\nSpecBoxplot3 = function(data, signal, parametric = FALSE, autor = "G. Anonim"){\n  if(parametric) {\n    pwc = data %>%\n      pairwise_t_test(value~COND, paired = TRUE,\n                      p.adjust.method = "bonferroni") %>%\n      add_xy_position(x = "COND") %>%\n      mutate(COND="NEG-CTR",\n             lab = paste(p, " - ", p.adj.signif))\n    res.test = data %>% anova_test(value~COND)\n  } else {\n    pwc = data %>% pairwise_wilcox_test(value~COND) %>%\n      add_xy_position(x = "COND") %>%\n      mutate(COND="NEG-CTR",\n             lab = paste(p, " - ", p.adj.signif))\n    res.test = data %>% kruskal_test(value~COND)\n  }\n  \n  data %>% ggplot(aes(COND, value, fill=COND))+\n    geom_violin(alpha=0.2)+\n    geom_boxplot(outlier.shape = 23,\n                 outlier.size = 3,\n                 alpha=0.6)+\n    geom_jitter(shape=21, width =0.1)+\n    stat_pvalue_manual(pwc, step.increase=0.05, label = "lab")+\n    ylab(signal)+\n    labs(title = get_test_label(res.test, detailed = TRUE),\n         subtitle = get_pwc_label(pwc),\n         caption = autor)\n}\n\n\n#special boxplot for the P3FCz signal\ndf$data[[1]] %>% SpecBoxplot3("P3FCz", TRUE)\ndf$data[[1]] %>% SpecBoxplot3("P3FCz", FALSE)\n\n#A function that creates a special boxplot3 and adds it to a data frame\nAddSignalBoxplot3 = function(df, signal, printPlot=TRUE) {\n  plot1 = SpecBoxplot3(df$data[[1]], signal, TRUE)\n  plot2 = SpecBoxplot3(df$data[[1]], signal, FALSE)\n  if(printPlot) print(plot1)\n  if(printPlot) print(plot2)\n  df %>% mutate(boxplot1 = list(plot1),\n                boxplot2 = list(plot2),\n  )\n}\n\n#Added special boxplot3\ndf %>% group_by(signals) %>%\n  group_modify(~AddSignalBoxplot3(.x, .y))\n
Run Code Online (Sandbox Code Playgroud)\n

以下是其运行结果\n在此输入图像描述

\n

在此输入图像描述

\n

归档时间：

查看次数：

627 次

最近记录：

2 年，8 月 前
相关归档
在特定的行索引处添加新行到dataframe,而不是追加？ 159
汇总/汇总每组的多个变量(例如总和,平均值) 143
推荐用于R中非常大的数据集处理和机器学习的软件包 41
<my code>出错:赋值目标扩展为非语言对象 28
使用R中的hist()函数获得百分比而不是原始频率 26
Roxygen很好的参考？ 20
R中的函数可以返回两个向量吗？ 19
扭转累积和的功能方法？ 18
迭代一次和迭代两次之间的性能差异？ 5
以ReSharper方式循环使用嵌套条件有什么好处？ 1
难疑归档
樱桃挑选Git意味着什么？ 2117
如何在JavaScript中将数字格式化为美元货币字符串？ 1711
确定整数平方根是否为整数的最快方法 1403
如何在HTTP POST请求中发送参数？ 1396
如何在Git中有选择地合并或选择来自另一个分支的更改？ 1374
何时使用struct？ 1347
什么是右值,左值,x值,glvalues和prvalues？ 1291
wait()和sleep()之间的区别 1158
我怎么知道分支是否已经合并为主分支？ 1077
#Hibernate hbm2ddl.auto配置的可能值是什么,它们做了什么 1046