Warning: file_get_contents(/data/phpspider/zhask/data//catemap/5/objective-c/25.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 为什么我会犯错误;未找到断点";在我的数据子集上的分段回归模型中?_R_Regression - Fatal编程技术网

R 为什么我会犯错误;未找到断点";在我的数据子集上的分段回归模型中?

R 为什么我会犯错误;未找到断点";在我的数据子集上的分段回归模型中?,r,regression,R,Regression,我有一个数据库,里面有在多个地点进行的许多不同实验的数据。我已经为每个实验确定了年龄对谷物产量的影响,现在我的导师希望我对一些选定的实验进行分段回归。我可以为整个模型找到一个断点,但是对于代表单个实验的模型,我得到了错误“找不到断点”。我想知道这是否是因为我在每个单独的实验中只有3年的数据(并且每年都在同一时间进行所有的观察),而在整个模型中,一些实验有4年的数据。如果您有任何见解,我将不胜感激——我已经尝试过搜索,但找不到任何有关此错误的故障排除提示,也找不到任何关于x轴最小范围的文档 dat

我有一个数据库,里面有在多个地点进行的许多不同实验的数据。我已经为每个实验确定了年龄对谷物产量的影响,现在我的导师希望我对一些选定的实验进行分段回归。我可以为整个模型找到一个断点,但是对于代表单个实验的模型,我得到了错误“找不到断点”。我想知道这是否是因为我在每个单独的实验中只有3年的数据(并且每年都在同一时间进行所有的观察),而在整个模型中,一些实验有4年的数据。如果您有任何见解,我将不胜感激——我已经尝试过搜索,但找不到任何有关此错误的故障排除提示,也找不到任何关于x轴最小范围的文档

data <- read.csv("C:/Users/.../Dropbox/WD/Yield decline/data_with_relative_columns.csv")
data2$Grain<- as.numeric(as.character(data2$Grain))
data2$Age <- as.numeric(as.character(data2$Age))
data2 <- subset(data2, Yrs >2) #I'm only interested in experiments with at least 3 years of data

data2$Key4 <- paste(data2$Experiment, data2$Site)

library(segmented)
yielddecline <- lm(Grain~ Age, data2)
seg <- segmented(yielddecline, seg.z= ~Age, psi=2)
summary(seg)

lambleg <- lm(Grain ~ Age, subset(data2, Key4 =="FGI- Legume Lamberton") )
seg.lambleg <- segmented(lambleg, seg.z = ~Age, psi=2)

morclip <- lm(Grain ~ Age, subset(data2, Key4 =="Ceres- Clipping Morris") )
seg.morclip <- segmented(morclip, seg.z = ~Age, psi=2)

mormob <- lm(Grain ~ Age, subset(data2, Key4 =="SARE- Mob Morris") )
seg.mormob <- segmented(mormob, seg.z = ~Age, psi=2)

morstrip <- lm(Grain ~ Age, subset(data2, Key4 =="SARE- Strip Morris") )
seg.morstrip <- segmented(morstrip, seg.z = ~Age, psi=2)

rosrot <- lm(Grain ~ Age, subset(data2, Key4 =="GMI- Rotation Rosemount") )
seg.rosrot <- segmented(rosrot, seg.z = ~Age, psi=2)

stpleg <- lm(Grain ~ Age, subset(data2, Key4 =="FGI- Legume Saint Paul") )
seg.stpleg <- segmented(stpleg, seg.z = ~Age, psi=2)
数据