@@ -12,7 +12,7 @@ library(stringr)
1212library(cowplot )
1313library(GGally )
1414
15- my_sim = function (parm1 , parm2 ){
15+ my_sim = function (parm1 , parm2 ) {
1616 my_obs = rnorm(n = 10000 , mean = parm1 , sd = parm2 )
1717 return (mean(my_obs ))
1818}
@@ -27,26 +27,22 @@ sim_data1 = expand.grid(mu = settingsA, sigma = settingsB)
2727set.seed(123 )
2828for (i in seq_len(length.out = nrow(sim_data1 ))) {
2929 sim_data1 $ observedmean [i ] = my_sim(
30- parm1 = sim_data1 $ mu [i ],
30+ parm1 = sim_data1 $ mu [i ],
3131 parm2 = sim_data1 $ sigma [i ]
3232 )
3333}
3434
3535# But we should use dplyr
3636set.seed(123 )
37- sim_data2 = sim_data1 | >
37+ sim_data2 = sim_data1 | >
3838 rowwise() | >
3939 mutate(observedmean = my_sim(parm1 = mu , parm2 = sigma ))
4040
41- sim_data1 | >
42- pivot_wider(names_from = mu , values_from = observedmean ) | >
43- gt :: gt() | >
41+ sim_data1 | >
42+ pivot_wider(names_from = mu , values_from = observedmean ) | >
43+ gt :: gt() | >
4444 gt :: tab_options(latex.use_longtable = TRUE )
45- sim_data2 | >
46- pivot_wider(names_from = mu , values_from = observedmean ) | >
47- gt :: gt() | >
45+ sim_data2 | >
46+ pivot_wider(names_from = mu , values_from = observedmean ) | >
47+ gt :: gt() | >
4848 gt :: tab_options(latex.use_longtable = TRUE )
49-
50-
51-
52-
0 commit comments