Trung tâm cá nhân
Bài viết phát hành
từ bỏ
Bài viết phổ biến của tác giả
- html - 出于某种原因，IE8 对我的 Sass 文件中继承的 html5 CSS 不友好？
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景？
25
4
我用对数线性模型拟合了一个 3 维列联表(这里没有提供，但如果它有帮助的话我可以)，包括 loglm 和 glm。我根据系数得到的两个结果是:
> coefficients(nodnox_loglm_model)
$`(Intercept)`
[1] 10.18939
$w
0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45
-1.04596513 -0.41193617 -0.08840858 0.06407334 -0.06862606 0.02999039 0.17084795 0.45838071 0.35307375
0,5
0.53856982
$s
2 3 4 5
0.36697307 0.15164360 -0.48264571 -0.03597096
Và
> coefficients(nodnox_glm_model)
(Intercept) s3 s4 s5 w0.1 w0.15 w0.2 w0.25 w0.3
9.5104005 -0.2153295 -0.8496188 -0.4029440 0.6340290 0.9575566 1.1100385 0.9773391 1.0759555
w0.35 w0.4 w0.45 w0.5
1.2168131 1.5043458 1.3990389 1.5845350
我知道这两种方法有不同的数值程序——我不关心这个——我只想知道如何将 glm 系数与 loglm 系数联系起来？
在使用 stackoverflow 之前，我在互联网上找到的所有内容和我搜索的文档都是这个注释:
The glm coefficient table works just like the summary for ANOVA produced by lm: the level alphabetically first (s2,w0.5) is used as an intercept, and all subsequent levels are tested against the first (thus the remaining coefficients are differences from the mean, not means themselves).
不过，对我来说，这还不足以理解如何以 loglm 的形式从 glm 输出中获取系数。现在，您的问题可能是:“为什么不直接使用 loglm？”Loglm 在我的情况下不起作用(这不是我在这里比较的那个，但它有一个带有一些零的 5 维表。所以如果我在原始表上使用 loglm，它会给我所有系数作为 NaN)。所以我坚持使用 glm，我真的很想获得 loglm 中的系数。
非常感谢!
1 Câu trả lời
您似乎有一个双向交叉表，其中包含 10 个水平的因子 w 和 5 个水平的因子 S，并且模型中没有交互作用。使用 glm()，分类变量的默认编码方案是 treatment coding其中因子中的第一组是引用水平，其余每个组的各自参数是其与该引用的差异。 (Intercept) 估计是针对所有组 = 其因子引用水平的单元格。
hiện hữuloglm()中，参数进行偏差编码，即每组都有自己的参数，一个因子的参数和为零。 (Intercept) 是添加到所有组效应的总均值。
在您的示例中，您可以告诉 glm() 使用偏差编码来获得与 loglm() 相同的参数估计值(请参见下面的示例)，或者您将处理编码的参数估计值转换如下:
w = 0.05 和 S = 2 是引用单元格:glm() 9.5104005 = loglm() 10.18939 + -1.04596513 + 0.36697307w = 0.1 和 S = 2 是 S 的引用水平，但需要与 w = 0.1 到引用 w = 0.05:glm() 9.5104005 + 0.6340290 = loglm() 10.18939 + -0.41193617 + 0.36697307w = 0.1 和 S = 3 但需要从 w = 0.1 到引用 w 的区别= 0.05 以及从 S = 3 到引用 S = 2 的差异:glm() 9.5104005 + 0.6340290 + -0.2153295 = loglm() 10.18939 + -0.41193617 + 0.15164360，等等使用偏差编码的 glm() 示例(UCBAdmissions 是一个具有内置于基 R 中的绝对频率的交叉表):
> library(MASS) # for loglm()
> llmFit <- loglm(~ Admit + Gender + Dept, data=UCBAdmissions)
> coef(llmFit)
$`(Intercept)`
[1] 5.177567
$Admit
Admitted Rejected
-0.2283697 0.2283697
$Gender
Male Female
0.1914342 -0.1914342
$Dept
A B C D E F
0.23047857 -0.23631478 0.21427076 0.06663476 -0.23802565 -0.03704367
> UCBdf <- as.data.frame(UCBAdmissions) # convert to data frame for glm()
> glmFit <- glm(Freq ~ Admit + Gender + Dept, family=poisson(link="log"),
+ contrasts=list(Admit=contr.sum, Gender=contr.sum, Dept=contr.sum),
+ data=UCBdf)
> coef(glmFit)
(Intercept) Admit1 Gender1 Dept1 Dept2 Dept3 Dept4
5.17756677 -0.22836970 0.19143420 0.23047857 -0.23631478 0.21427076 0.06663476
Dept5
-0.23802565
Xin lưu ý,glm() 并未列出那些通过一个因子的参数的总和为零约束完全确定(别名)的参数估计值。
关于r - glm 中的系数与 loglm 中的系数，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/19013793/
25
4
0
我已经训练了一个线性回归模型来以多输出方式进行预测。这是一个时间序列预测问题，根据一组输入来估计 future 12 个月的需求。在过去 - 如果我只预测一个输出值 - 我会简单地调用以下命令来访问模
我的应用程序对原始音频信号执行FFT(所有麦克风读数均为values中的16位整数值，即1024个单元)。它首先根据16位对读数进行归一化。然后提取频率为400Hz的幅度。 int sample_ra
我从包含相同语音内容的两个大约 30 秒的音频文件中提取了两个系列 MFCC 系数。音频文件从不同来源录制在同一位置。应该估计音频是否包含相同的对话或不同的对话。目前我已经测试了两个Mfcc系列的相关
我正在尝试使用两个自变量(批处理和方法)和一个因变量 (conc) 运行线性回归。当我运行回归时，我得到一个系数的 NA 值。当我更改模型中自变量的顺序时，NA 值显示为不同的变量。这是一个可重现的数
我想要的，大图:我想知道如何在 Java 中模仿 Mathematica 的 ExtendedGCD[...] 功能。有关该功能的信息可以在 here 中找到，但为了完整起见，我将对其进行简要描述。
我用 R 编写了这个套索代码，并且得到了一些 beta 值: #Lasso library(MASS) library(glmnet) Boston=na.omit(Boston) x=model.m
我有一个幅频响应点列表。振幅以分贝为单位给出。 任务是将此 AFR 列表导出为某些硬件 DSP 设备的一组系数。我知道此设备的文件格式。该设备的规范表明它需要 12288 个系数，并且 DSP 有 2
是否可以替换 lm 对象中的系数？ 我认为以下会起作用 # sample data set.seed(2157010) x1 <- 1998:2011 x2 <- x1 + rnorm(length(
我正在尝试使用 lm.ridge 方法执行岭回归。我的问题是如何获得拟合模型的系数？我从调用中得到不同的结果: 模型$coef 系数(模型) 哪个是正确的？另外，为什么我从调用中得到不同的结果: co
以下代码片段死于 Eigen 断言: MatrixXd L; VectorXd x, b; ... ASSERT_MATRIX_EQ(L*x, b); 与， template void ASSERT
有没有办法从编码的图像和视频中轻松提取 DCT 系数(和量化参数)？任何解码器软件都必须使用它们来解码 block DCT 编码的图像和视频。所以我很确定解码器知道它们是什么。有没有办法将它们暴露给使
如何找到矢量 w 是多少，即垂直于分离平面？ 最佳答案 我就是这样做的here .如果我没记错的话，这是基于 dual form 的SVM 优化的结果。 model = svmtrain(...);
我是 FFMPEG 的新手。我知道如何提取关键帧， 但现在我想提取那些 I 帧的 DC 系数。 您能否帮助编写 FFMPEG 命令的正确语法以获取 DC 系数？ 问候 最佳答案 FFmpeg 命令不会
任何人都可以共享任何代码或可调用以从 HEVC 比特流中提取变换系数 (DCT/DST) 的函数吗？ 我检查了 HM 16.0 的源代码。文件\source\Lib\TLibCommon\TComTr
我需要从 MPEG4 视频中提取有关运动矢量和 DC 系数的信息。我搜索了相关站点和主题，发现一个好的解决方案是使用 ffmpeg 编解码器的代码。尤其是 libavcodec/mpegvideo.c
我想删除sympy符号表达式中等于1的任何系数，例如:我希望 1.0x**2 成为 x**2 ，有办法做到吗？另外，如果可能的话，可以对整数进行舍入，例如将 2.0x**2 变为 2*x**2 最佳答
我想删除sympy符号表达式中等于1的任何系数，例如:我希望 1.0x**2 成为 x**2 ，有办法做到吗？另外，如果可能的话，可以对整数进行舍入，例如将 2.0x**2 变为 2*x**2 最佳答
我最近的任务是将 PMML 解析回 R 模型。 (我进行了广泛搜索，没有库可以为您进行这种转换。)我正在尝试将包含多项逻辑回归的 PMML 转换回 R 模型，但我不知道如何转换任何PMML 文档中保存
当使用 pandas statsmodels 的 .summary() 函数时，OLS 回归结果包括以下字段。 coef std err t P>|t| [
我正在 STM32F4 微 Controller 上使用 CMSIS 库开发 PID Controller 。我实在无法理解PID参数的归一化。现在我有了 PID Controller 的三个 Kp、
Tôi là một lập trình viên xuất sắc, rất giỏi!