亲爱的Stack Overflow社区,
目前我正试图在最新版本的R和lme4上重新运行旧数据分析,二项式glmer模型(从2013年初开始),因为我不再使用旧版本的R和lme4了.但是,我通过dmartin和carine(第一个警告消息)以及堆栈溢出之外的其他线程(警告2和3)经历与先前线程类似的警告消息.我使用的早期版本的R和lme4上没有弹出这些警告消息,所以它必须与最新更新有关？
我的数据集的一个子集:
df <- structure(list(SUR.ID = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 1L,
1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L,
3L, 1L, 1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L,
2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L, 3L, 1L,
1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L, 3L, 1L, 1L, 2L, 2L, 3L,
3L, …我通过论坛阅读了如何解决我的问题的线索,但是没有任何相关的线程可用于我,具有有限的编程知识,适用于我的具体问题.
我的问题是:我需要摆脱整个文件中聚集的垃圾行,但是在可用行集群之间.我搜索了有关删除匹配模式的范围的sed手册和其他信息来源,但他们只提到删除UNTIL匹配模式,而不是TILL.
现在我想指定一个范围,sed从第一行开始删除与模式行匹配的行,直到与其他模式匹配的行.此外,sed需要识别行末的模式.
例如:
line 1
blah blah 1
blah blah 2
blah blah 3
blah blah 4
line 2
line 3
结果需要是:
line 1
blah blah 1
line 2
line 3
请注意线和线2之间的多条线.虽然需要保留blah blah 1,但需要删除其他3条线.
谢谢!
我一直在努力将缩放和居中的模型系数从glmer模型转换回未经中断和未缩放的值.
我在lme4(v1.1.7)包中使用GLMM分析了一个数据集.它涉及声学接收器的最大检测范围的计算和环境变量的影响.
样本数据:
dd <- structure(list(SUR.ID = c(10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L, 10186L,
10186L, 10186L, 10186L, 10186L, …我正在尝试在 RStudio 中安装gdalUtils，但它抛出一个错误，指出它与我的 R v 4.2.0 兼容。（最新）不存在。我需要这个包才能使用另一个包。有什么建议如何解决这个问题吗？
Warning in install.packages :\n package \xe2\x80\x98gdalUtils\xe2\x80\x99 is not available for this version of R\n\nA version of this package for your version of R might be available elsewhere,\nsee the ideas at\nhttps://cran.r-project.org/doc/manuals/r-patched/R-admin.html#Installing-packages\n我一直试图解决这个问题超过2天,筛选互联网和教程,但我还没有解决我的问题.我是一个真正的新手,还不知道我在做什么..
我使用的软件:Mac OS X 10.6 Python v3.2.2交互式解释器(IDLE)
问题:IDLE的默认目录是/ Users/ME/Documents /.具有扩展名.py的文件只能在此目录中打开.但是,我创建了一个文件夹,我想保存所有与此软件有关的.py文件等.目前,IDLE无法从我选择的目录中加载.py文件.
我首先做的是我添加到IDLE:import sys.sys.path.append('Users/Mydir /')sys.path
但是,在2010年的现有线程中,我读取sys.path仅用于解释器,如果我要更改它,我需要修改PYTHONPATH环境变量:
PYTHONPATH ="/ Me/Documents/mydir:$ PYTHONPATH"导出PYTHONPATH
但是,我很困惑如何使用它并且无法找到以下问题的答案:1)安装程序时,我的计算机上已经存在PYTHONPATH(.py？)？如果是,那它在哪里？我找不到任何地方.如果不是,我需要创建一个.但是内容应该在哪里以及应该是什么,以便IDLE可以从非默认目录加载文件？它应该只包含粗体字吗？
我希望我能清楚地解决问题.
干杯
我已经搜索了SO,以便在没有运气的情况下实现我的需要,所以在这里.前段时间我发现了包装dplyr及其潜力.我在想这个包可以做我想做的事,我只是不知道怎么做.这是我数据的一小部分,但应该代表我的问题.
dummy<-structure(list(time = structure(1:20, .Label = c("2015-03-25 12:24:00",
"2015-03-25 21:08:00", "2015-03-25 21:13:00", "2015-03-25 21:47:00",
"2015-03-26 03:08:00", "2015-04-01 20:30:00", "2015-04-01 20:34:00",
"2015-04-01 20:42:00", "2015-04-01 20:45:00", "2015-09-29 18:26:00",
"2015-09-29 19:11:00", "2015-09-29 21:21:00", "2015-09-29 22:03:00",
"2015-09-29 22:38:00", "2015-09-30 00:48:00", "2015-09-30 01:38:00",
"2015-09-30 01:41:00", "2015-09-30 01:45:00", "2015-09-30 01:47:00",
"2015-09-30 01:49:00"), class = "factor"), ID = c(1L, 1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
2L), station = c(1L, 1L, 1L, 2L, 3,
4L, 4L, …假设我有以下栅格：
library(raster)
r <- raster(ncols=36, nrows=18)
values(r) <- rnorm(ncell(r))
我可以使用以下方法离散其值cut()：
breaks <- -2:2 * 3
rc <- cut(r, breaks=breaks)
但是，当我尝试指定标签，例如 A、B、C 等时
rc <- cut(r, breaks=breaks, labels = c("A", "B", "C", "D", "E"))
它返回以下错误。
library(raster)
r <- raster(ncols=36, nrows=18)
values(r) <- rnorm(ncell(r))
我尝试用谷歌搜索错误消息，但没有成功。有什么想法出了什么问题吗？
对于出版物,我需要在现有图中添加第二个y轴.我遇到过如何做到这一点的方法(https://rpubs.com/kohske/dual_axis_in_ggplot2).但是,我并不太了解编码.我找不到一种方法来使它成为正确的y轴也显示,而不是顶部边框.我的编码中缺少什么？这是我的虚拟数据:
df1 <- structure(list(month = structure(1:12, .Label = c("Apr", "Aug",
"Dec", "Feb", "Jan", "Jul", "Jun", "Mar", "May", "Nov", "Oct",
"Sep"), class = "factor"), RI = c(0.52, 0.115, 0.636666666666667,
0.807, 0.66625, 0.34, 0.143333333333333, 0.58375, 0.173333333333333,
0.5, 0.13, 0), sd = c(0.327566787083184, 0.162634559672906, 0.299555225848813,
0.172887246493199, 0.293010848165827, 0.480832611206852, 0.222785397486759,
0.381610777775321, 0.219393102292058, 0.3, 0.183847763108502,
0)), .Names = c("month", "RI", "sd"), class = "data.frame", row.names = c(NA,
-12L))
df2<-structure(list(month = structure(c(5L, 4L, 8L, 1L, 9L, 7L, 6L,
2L, 12L, 11L, 10L, 3L), .Label …我想执行引导以获得二项式 GLMM 中固定效应的 95% 顺式：
m <- glmer(cbind(df$Valid.detections, df$Missed.detections) ~ distance +
Habitat + Replicate + transmitter.depth + receiver.depth +
wind.speed + wtc + Transmitter + (1 | Unit) +
(1 | SUR.ID) + distance:Transmitter +
distance:Habitat + distance:transmitter.depth + distance:receiver.depth +
distance:wind.speed, data = df, family = binomial(link=logit),control=glmerControl(calc.derivs=F))
我发现confint()函数可以实现这个，所以我指定了函数：
confint(m, method = "boot", boot.type = "basic", seed = 123, nsim = 1000)
在我决定终止之前，该函数已经运行了 8 个多小时。终止后，我收到以下警告消息 (x10)：
Warning messages:
1: In (function (fn, par, lower = rep.int(-Inf, n), …