我正在尝试做一个练习,以便更熟悉如何在purrr中使用map函数.我正在创建一些随机数据(10列10个数据点),然后我想使用map在数据框中的结果列上执行一系列回归(即lm(y~x,data =)).
如果我只是重复使用第一列为'y',我想要执行10次回归,每列从1到10为'x'.显然结果并不重要 - 它只是方法.我想最终得到10个线性模型对象的列表.
list_of_vecs <- list()
for (i in 1:10){
list_of_vecs[[paste('vec_', i, sep = '')]] <- rnorm(10,0,1)
}
df_of_vecs <- as.data.frame(list_of_vecs)
Run Code Online (Sandbox Code Playgroud)
在这里,我卡住了:
map(df_of_vecs, ~ lm(df_of_vecs[[1]] ~ . ?)
Run Code Online (Sandbox Code Playgroud)
任何提示将不胜感激.
谢谢.
我想要更多地使用 R markdown 来执行分析和生成输出。也许我错过了一些简单的东西,但我只想能够设置小数位数以显示 2 或 3 位数字,具体取决于输出(例如 t 统计量与 p 值)。
我以前使用过r options(digits=2),它一直有效,直到您要包含的最后一个数字为 0。我已经使用 sprintf 函数解决了这个问题,但必须为每个数字指定。
有没有办法设置“全局”sprintf 选项,以便对于后面的所有数字,显示相同的小数位数?
谢谢你,
保罗
我试图通过调整 alpha 值来复制一些纵向数据的轨迹图,并突出显示随机主题轨迹。我已经完成了 95% 的工作,除了当 alpha = 0 时我似乎无法获得线图的完全透明度(它仍然显示为灰色)。有任何想法吗?- 数据在这里:
谢谢你。
# Assign random alpha (0 [90%] or 1 [10%]) values
macs <- ddply(macs, .(id), function(x){
x$alpha = ifelse(runif(n = 1) > 0.9, 1, 0)
x
})
# Loess fit
fit <- loess(CD4 ~ time, data=macs, span = 0.5)
# Plot
ggplot(data = macs, aes(x = fitted(fit), y = resid(fit))) +
geom_point(aes(alpha=1), size = 0.5) +
geom_line(aes(alpha=alpha, group=id)) +
guides(alpha=FALSE)
Run Code Online (Sandbox Code Playgroud)
我希望有人能帮我解决这个问题.数据可以在这里找到.
基本上我想绘制一条直线(平均值)和它测试的4个模型的相关置信区间(下限,上限).我想面对Cat_Auth变量,其中有4个类别(所以4个图).第一个"模型"实际上只是样本数据的平均值,我不想要CI(数据中指定的NA值 - 不确定这是否正确).
我可以通过以下方式获得情节:
newdata <- read.csv("data.csv", header=T)
ggplot(newdata, aes(x = Affil_Max, y = Mean)) +
geom_line(data = newdata, aes(), colour = "blue") +
geom_ribbon(data = newdata, alpha = .5, aes(ymin = Lower, ymax = Upper, group = Model, fill = Model)) +
facet_grid(.~ Cat_Auth)
Run Code Online (Sandbox Code Playgroud)
但我想为每个型号设置不同的彩色线条和阴影色带(例如,红色平均线和红色阴影色带用于模型2,绿色用于模型3等).另外,我无法弄清楚为什么对应于第一组平均值的蓝线原样是脱节的.
真的很感激任何帮助!
我在安装时遇到了问题Flextable,这似乎是由于gdtools包装所致。看起来gdtools依赖于.dylib文件吗?有任何想法吗?
> install.packages("flextable")
There is a binary version available but the source version is later:
binary source needs_compilation
flextable 0.4.0 0.4.2 FALSE
installing the source package ‘flextable’
trying URL 'https://cran.rstudio.com/src/contrib/flextable_0.4.2.tar.gz'
Content type 'application/x-gzip' length 135276 bytes (132 KB)
==================================================
downloaded 132 KB
Warning in strptime(xx, f <- "%Y-%m-%d %H:%M:%OS", tz = tz) :
unknown timezone 'default/Australia/Melbourne'
* installing *source* package ‘flextable’ ...
** package ‘flextable’ successfully unpacked and MD5 sums checked
** R
** …Run Code Online (Sandbox Code Playgroud)