问题描述
大家。我在 R 中使用 DLNM 来分析气候条件对疾病流行的滞后效应。 我严格按照别人的程序 ，它在 avg.temp 和 max.speed 中工作，但在 avg.ap 和 avg.hum 中显示错误“coef/vcov 与基础矩阵不一致”。但是，我只是更改了代码中设置的变量，并没有更改其他代码。 我有一个假设，也许 DLNM 不喜欢潮湿的天气。 T T 我不知道该怎么办，你能帮我吗？
第 1 部分是成功运行的代码，第 2 部分是显示错误的代码，第 3 部分是我使用的数据。 非常感谢你。希望你能帮我
第 1 部分。成功运行代码
attach(cpdlnm)
cb.temp = crossbasis(avg.temp,lag=1,argvar=list(fun="ns",knots= c(10)),arglag=list(fun="lin"))
modeltemp = glm(pre1 ~ cb.temp +
ns(no,1*1),family=quasipoisson(),cpdlnm)
pred1.temp = crosspred(cb.temp,modelhum,cen=round(median(avg.temp)),bylag=1)
attach(cpdlnm)
cb.hum = crossbasis(avg.hum,arglag=list(fun="lin"))
modelhum = glm(pre1 ~ cb.hum +
ns(no,cpdlnm)
pred1.hum = crosspred(cb.hum,# This step shows "coef/vcov not consistent with basis matrix"
modelhum,cen=round(median(avg.hum)),bylag=0.1)
第三部分数据如下：
no pre1 date year month avg.ap avg.temp avg.hum max.speed
1 3.23 12-Jan 2012 1 996.60 9.00 81.60 5.30
2 6.04 12-Feb 2012 2 993.20 10.90 80.80 6.20
3 5.18 12-Mar 2012 3 991.00 16.40 78.70 7.60
4 4.07 12-Apr 2012 4 985.40 23.50 73.50 7.40
5 4.88 12-May 2012 5 982.60 26.30 77.20 7.00
6 5.11 12-Jun 2012 6 978.10 27.00 81.30 6.20
7 6.18 12-Jul 2012 7 979.50 28.10 77.70 6.40
8 6.17 12-Aug 2012 8 980.40 28.00 75.60 7.90
9 5.18 12-Sep 2012 9 987.60 25.30 73.60 6.30
10 5.16 2012 年 10 月 12 日 10 990.70 23.60 72.20 6.20 11 4.61 2012 年 11 月 12 日 11 991.70 18.00 79.70 6.90 12 5.26 2012 年 12 月 12 日 12 995.00 13.20 74.90 6.50 13 3.79 2013 年 1 月 13 日 1 997.10 11.20 78.40 5.70 14 3.87 2013 年 2 月 13 日 2 993.50 15.30 82.20 6.50 15 3.37 2013 年 3 月 13 日 3 989.90 20.20 74.20 8.00 16 2.85 2013 年 4 月 13 日 4 987.00 21.50 78.50 7.70 17 4.38 2013 年 5 月 13 日 5 983.30 25.60 79.20 6.80 18 5.67 2013 年 6 月 13 日 6 980.60 27.40 76.90 6.60 19 6.45 2013 年 12 月 13 日 7 981.30 28.00 77.50 7.10 20 6.95 2013 年 8 月 13 日 8 980.50 27.90 78.20 7.90 21 6.51 2013 年 9 月 13 日 9 985.90 25.40 77.60 6.00 22 8.16 2013 年 10 月 13 日 10 992.20 22.10 68.80 5.30 23 5.34 2013 年 11 月 13 日 11 994.50 18.70 72.30 6.20 24 6.18 2013 年 12 月 13 日 12 997.30 11.70 67.20 5.30 25 5.69 2014 年 1 月 14 日 1 996.70 12.70 70.30 6.00 26 6.44 2014 年 2 月 14 日 2 993.00 12.10 76.90 6.40 27 4.16 2014 年 3 月 14 日 3 991.60 16.50 83.90 7.30 28 4.13 2014 年 4 月 14 日 4 987.60 22.60 82.40 6.70 29 3.96 2014 年 5 月 14 日 5 983.60 25.70 78.80 7.70 30 4.72 2014 年 6 月 14 日 6 979.20 27.70 81.40 7.90 31 5.21 2014 年 7 月 14 日 7 980.70 28.30 80.20 9.40 32 5.29 2014 年 8 月 14 日 8 982.40 27.50 81.30 7.50 33 6.74 2014 年 9 月 14 日 9 984.70 27.10 77.70 8.50 34 4.80 2014 年 10 月 14 日 10 991.20 23.90 73.10 5.90 35 4.31 2014 年 11 月 14 日 11 993.30 18.60 79.60 6.20 36 4.35 2014 年 12 月 14 日 12 998.70 12.30 67.30 5.90 37 2.95 2015 年 1 月 15 日 1 996.70 13.30 76.30 6.20 38 4.63 2015 年 2 月 15 日 2 993.50 15.50 78.30 6.50 39 4.00 2015 年 3 月 15 日 3 991.70 17.70 83.40 6.30 40 4.16 2015 年 4 月 15 日 4 988.40 22.80 70.20 7.30 41 4.67 2015 年 5 月 15 日 5 982.40 26.70 80.50 8.00 42 5.62 2015 年 6 月 15 日 6 980.90 28.20 81.00 7.40 43 5.04 2015 年 12 月 13 日 7 980.20 27.30 79.40 6.70 44 5.79 2015 年 8 月 15 日 8 982.40 27.60 80.10 6.50 45 5.28 2015 年 9 月 15 日 9 986.30 26.00 84.60 6.50 46 4.39 2015 年 10 月 15 日 10 991.20 23.00 78.30 6.90 47 4.13 2015 年 11 月 15 日 11 993.50 19.40 85.30 6.90 48 3.30 2015 年 12 月 15 日 12 997.80 13.00 80.90 5.70 49 5.30 2016 年 1 月 16 日 1 996.00 11.80 82.30 6.40 50 4.57 2016 年 2 月 16 日 2 997.80 12.20 68.90 7.00 51 4.66 2016 年 3 月 16 日 3 991.70 17.00 78.90 7.00 52 4.01 2016 年 4 月 16 日 4 984.60 23.40 80.90 9.80 53 4.90 2016 年 5 月 16 日 5 983.80 25.50 78.70 8.30 54 3.75 2016 年 6 月 16 日 6 981.70 28.20 78.80 7.70 55 3.13 2016 年 12 月 13 日 7 981.10 28.90 77.60 7.60 56 3.25 2016 年 8 月 16 日 8 979.00 28.00 79.80 8.70 57 2.93 2016 年 9 月 16 日 9 984.30 26.60 75.20 6.40 58 2.93 2016 年 10 月 16 日 10 987.90 24.40 72.90 7.00 59 3.08 2016 年 11 月 16 日 11 993.40 18.10 79.60 6.70 60 2.99 2016 年 12 月 16 日 12 995.70 15.40 71.70 6.80 61 3.10 2017 年 1 月 17 日 1 994.70 14.50 79.20 6.50 62 3.75 2017 年 2 月 17 日 2 994.80 14.70 71.50 8.30 63 3.49 2017 年 3 月 17 日 3 990.20 16.50 83.60 8.50 64 3.36 2017 年 4 月 17 日 4 986.80 21.90 76.70 7.80 65 3.69 2017 年 5 月 17 日 5 985.00 24.80 77.50 10.00 66 3.76 2017 年 6 月 17 日 6 980.20 26.90 84.80 8.50 67 2.69 2017 年 7 月 17 日 7 981.00 27.50 83.60 9.80 68 3.05 2017 年 8 月 17 日 8 980.50 27.70 83.40 9.00 69 3.05 2017 年 9 月 17 日 9 984.20 27.60 81.50 7.10 70 2.46 2017 年 10 月 17 日 10 990.00 22.80 75.90 7.90 71 2.08 2017 年 11 月 17 日 11 993.00 17.80 79.50 7.00 72 2.32 2017 年 12 月 17 日 12 996.90 13.30 69.30 6.90 73 2.53 2018 年 1 月 18 日 1 992.10 12.00 78.40 8.10 74 3.29 2018 年 2 月 18 日 2 992.90 13.40 68.70 7.20 75 3.03 2018 年 3 月 18 日 3 988.30 19.20 78.20 9.10 76 2.30 2018 年 4 月 18 日 4 986.50 21.80 77.30 8.70 77 1.75 2018 年 5 月 18 日 5 982.60 26.70 79.40 8.90 78 2.03 2018 年 6 月 18 日 6 978.30 26.90 81.60 9.00 79 2.79 2018 年 7 月 18 日 7 976.80 27.90 82.10 9.20 80 2.32 2018 年 8 月 18 日 8 976.40 27.50 83.40 9.60 81 1.88 2018 年 9 月 18 日 9 983.50 26.10 80.10 8.90 82 2.76 2018 年 10 月 18 日 10 990.50 21.10 78.70 7.10 83 2.14 2018 年 11 月 18 日 11 991.50 18.20 80.30 7.10 84 1.78 2018 年 12 月 18 日 12 994.50 13.00 84.00 7.80 85 2.77 2019 年 1 月 19 日 1 995.20 11.70 84.50 7.30 86 4.60 2019 年 2 月 19 日 2 990.50 13.70 84.80 8.10 87 2.32 2019 年 3 月 19 日 3 987.70 17.30 85.90 9.90 88 2.07 2019 年 4 月 19 日 4 983.60 23.10 84.80 9.80 89 2.97 2019 年 5 月 19 日 5 981.80 24.30 83.20 7.70 90 2.48 2019 年 6 月 19 日 6 977.80 27.50 84.80 9.00 91 2.32 2019 年 7 月 19 日 7 977.20 27.80 85.00 8.90 92 2.06 2019 年 8 月 19 日 8 977.20 28.30 81.20 10.30 93 2.10 2019 年 9 月 19 日 9 984.60 26.40 72.70 8.20 94 2.89 2019 年 10 月 19 日 10 989.10 22.70 78.00 7.00
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）