在R中找到两个线性拟合相交的位置
我在R脚本中使用lm函数得到了两个线性拟合。例如... fit1 <- lm(y1 ~ x1) fit2 <- lm(y2 ~ x2) 如果这两条线（fit1 和 fit2）有交点，我希望找到它们相交的 (x,y) 坐标。
如何在高阶参数保持不变的情况下删除模型中的低阶参数？
问题：只要高阶参数（例如交互作用）仍然存在于模型中，我就无法删除低阶参数（例如主效应参数）。即使这样做，模型也会被重构，新模型不嵌套在更高的模型中。 请参见以下示例（因为我来自 ANOVA，所以使用 contr.sum）： d <- data.frame(A = rep(c("a1",...
如何手动设置线性模型中变量的系数？
在R中，如何在lm()函数中为特定变量设置权重而不是观测值？ 背景如下。我正在尝试构建针对特定产品（例如手机）的个人排名系统。我可以基于价格作为因变量，其他特征（如屏幕大小、内存、操作系统等）作为自变量来构建线性模型。然后，我可以使用它来预测电话的真实成本（与公布价格相反），从而找到最佳价...
lm()返回的对象是否有函数或包可以模拟预测？
是否有一个类似于“runif”、“rnorm”等函数的单一函数，可以为线性模型生成模拟预测？我能够自己编写代码，但是代码很丑陋，我认为肯定有人做过这样的事情。 slope = 1.5 intercept = 0 x = as.numeric(1:10) e = rnorm(10, mean=...
`lm` 摘要未显示所有因子水平。
我正在对一些属性进行线性回归，其中包括两个分类属性B和F，但我并没有得到每个因子水平的系数值。 B有9个水平，F有6个水平。当我最初运行模型（带截距）时，我得到了B的8个系数和F的5个系数，这意味着每个属性的第一个水平被包含在截距中。 我想根据它们的系数对B和F中的水平进行排名，因此我在每...
在lm公式中删除变量仍会触发对比错误。
我正在尝试在数据的一个子集上运行lm()，但遇到了问题。 dt = data.table(y = rnorm(100), x1 = rnorm(100), x2 = rnorm(100), x3 = as.factor(c(rep('men',50), rep('women',50)))) ...
lm在dlply内被调用时会抛出“0（非NA）个案例”的错误 [R]
我正在使用dlply()函数和一个自定义函数，对包含一些NA值的数据进行lm()拟合斜率的平均值计算。但是，当我调用具有两个关键变量的dlply时，会出现以下错误："在lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...)中发...
在lm()中强制某些参数具有正系数
我想知道如何在lm()中约束某些参数以具有正系数。有一些包或函数（例如display）可以使所有系数和截距都为正。 例如，在这个例子中，我只想强制x1和x2具有正系数。 x1=c(NA,rnorm(99)*10) x2=c(NA,NA,rnorm(98)*10) x...
从回归估计中提取交互项
这是一个简单的问题，但我找不到任何清晰而令人信服的答案。如果我有一个带有一个或多个交互项的回归模型，例如： mod1 <- lm(mpg ~ factor(cyl) * factor(am), data = mtcars) coef(summary(mod1)) ## ...