我想知道是否有人可以帮助我理解为什么当我将脚本输入R时我收到错误消息.为了获得背景信息,我正在研究效果6个不同的变量(我认为是63种组合或模型)(X )对于我的环境科学荣誉项目,在不同空间尺度上单独生成初级和净生态系统总生产量(Y).我决定使用详尽的搜索多元回归分析和akaikes信息标准(AIC)来尝试找到一组最适合的模型.(以及用于比较归因于不同X变量的方差的分层分区)我想得到权重,因此我可以对哪些模型"最符合"标准进行排序,看看是否有一个或一组装备其余的,因此是更多可能适合数据.
我最近在关于Cross Validated的hier.part包上发布了一个类似的问题得到了一个很好的答案,如果我将来有任何类似的问题,我被告知来这里.
我用于R的包是glmulti.可以在这里找到
我正在使用的脚本就是这个
require(glmulti)
GPPANDDRIVER<-read.table("C:\\Databases at different scales for R\\River Rhine and Netherlands\\GPP and drivers rhineland (comma delimited).csv",header=T,sep=",")
GPP<-GPPANDDRIVER$GPP
IND_VARS<-subset(GPPANDDRIVER,select=-GPP)
# glmulti S4 generic
glmulti(y=GPP, xr=IND_VARS, data, exclude = c(), name = "glmulti.analysis", intercept = TRUE, marginality = FALSE, bunch=30, chunk = 1, chunks = 1,
level = 2, minsize = 0, maxsize = -1, minK = 0, maxK = -1, method = "h", crit = "aic", confsetsize = 63, popsize = 40, mutrate …Run Code Online (Sandbox Code Playgroud) 我在R中使用ggtern绘制三元图,并希望在ggtern图上具有轴标签和断点,与原始数据相同。对于以下代码中生成的数据,每个轴的最大值为12、10和4。
在上一篇文章之后,我尝试使用中断和标签来执行此操作,但是每个轴仍在0-1比例尺上,缺少标签(由于它们超过1),并且带有标签的轴线不与点相交在情节上。(如何更改ggtern制作的三元图的标签?)
library(ggtern)
labFnc <- function(x,digits=2) format(round(unique(x),digits),digits=digits)
mydata <- data.frame(
x = runif(50, min = 0.25, max = 12),
y = runif(50, min = 0.1, max = 10),
z = runif(50, min = 0.5, max = 4),
value = runif(50, min = 10000, max = 20000))
ggtern(data = mydata,aes(x = x, y = y, z = z,col=value)) +
theme_bw() +
geom_point(alpha = 0.8, size = 3) +
theme_showarrows() +
scale_T_continuous(breaks=unique(mydata$x),labels=labFnc(mydata$x))+
scale_L_continuous(breaks=unique(mydata$y),labels=labFnc(mydata$y))+
scale_R_continuous(breaks=unique(mydata$z),labels=labFnc(mydata$z))
Run Code Online (Sandbox Code Playgroud)
有没有办法做到这一点?任何帮助将不胜感激。
编辑:我也尝试添加tern_limits参数。虽然这看起来会按比例扩大绘图,但数据放置在错误的位置。而且我无法像以前一样添加自己的独特休息时间。
ggtern(data = mydata,aes(x = …Run Code Online (Sandbox Code Playgroud)