| 查看: 3430 | 回复: 15 | ||||
monitor2885至尊木虫 (职业作家)
队长
|
[求助]
求广义线性模型GLM的函数表达式 已有1人参与
|
|
y=b0+b1*x1+b2*x2+b3*sin(x3) y=[7.1;4.1;5.2;6.3]; x1=[1;2;3;4]; x2=[100;205;306;150]; x3=[8;10;12;5]; link 函数是 normal，求 glmfit 的表达式，求 b0, b1, b2, b3，谢谢。 |
回复此楼» 猜你喜欢
0856，269分求调剂
已经有18人回复
-
292求调剂
已经有19人回复
-
085601，一志愿厦大334复试被刷求调剂
已经有10人回复
-
生物学308分求调剂（一志愿华东师大）
已经有7人回复
-
学硕机械工程303求调剂
已经有7人回复
-
338求调剂
已经有7人回复
-
285求调剂
已经有5人回复
-
调剂
已经有4人回复
-
材料专硕调剂
已经有18人回复
-
348分环境工程·调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【教程】《实用统计分析及其DPS数据处理系统》
已经有12人回复
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 仿真EPI: 1
- 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
2楼2015-10-29 09:52:11
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39476
- 散金: 147
- 红花: 22
- 帖子: 4994
- 在线: 2174.2小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
3楼2015-10-29 09:55:10
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 仿真EPI: 1
- 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
4楼2015-10-29 09:55:11
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39476
- 散金: 147
- 红花: 22
- 帖子: 4994
- 在线: 2174.2小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
5楼2015-10-29 09:57:40
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 仿真EPI: 1
- 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
monitor2885: 金币+10, ★★★很有帮助 2015-10-29 10:42:13
monitor2885: 金币+10, ★★★很有帮助 2015-10-29 10:42:13
|
实际上b = glmfit([x1 x2 sin(x3)],y)也就可以出结果了，我猜想广义线性模型中的Y=b0+b1*x1+b2*x2+b3*x3，其中用sin(x3)来代替x3这才适用模型，所以你最开始b = glmfit([x1 x2 x3],y,'normal','link','normal')这样直接用x3就有问题，不符合模型假设 b = glmfit([x1 x2 sin(x3)],y,'normal','link','identity')这个格式也是可以的 函数格式为b = glmfit(X,y,distr,param1,val1,param2,val2,...)，distr can be any of the following strings: 'binomial', 'gamma', 'inverse gaussian', 'normal' (the default), and 'poisson'. b = glmfit(X,y,distr,param1,val1,param2,val2,...) additionally allows you to specify optional parameter name/value pairs to control the model fit. Acceptable parameters are as follows. 从matlab函数解释中看出，param参数是来优化distr分布的，你要注意写全了，要么就都缺省让函数内在保持一致，用b = glmfit([x1 x2 sin(x3)],y)，要么就都写全了，用b = glmfit([x1 x2 sin(x3)],y,'normal','link','identity') |
6楼2015-10-29 10:27:09
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 仿真EPI: 1
- 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
7楼2015-10-29 10:30:31
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39476
- 散金: 147
- 红花: 22
- 帖子: 4994
- 在线: 2174.2小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
|
你看stats，哪个是significance? stats is a structure with the following fields: beta — Coefficient estimates b dfe — Degrees of freedom for error s — Theoretical or estimated dispersion parameter sfit — Estimated dispersion parameter se — Vector of standard errors of the coefficient estimates b coeffcorr — Correlation matrix for b covb — Estimated covariance matrix for B t — t statistics for b p — p-values for b resid — Vector of residuals residp — Vector of Pearson residuals residd — Vector of deviance residuals resida — Vector of Anscombe residuals |
8楼2015-10-29 10:44:36
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 仿真EPI: 1
- 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
9楼2015-10-29 10:57:14
monitor2885
至尊木虫 (职业作家)
队长
- 应助: 6 (幼儿园)
- 贵宾: 0.339
- 金币: 39476
- 散金: 147
- 红花: 22
- 帖子: 4994
- 在线: 2174.2小时
- 虫号: 701577
- 注册: 2009-02-15
- 性别: GG
- 专业: 生态系统生态学
10楼2015-10-29 11:01:46
