みんな。私は R で DLNM を使用して、病気の有病率に対する気候条件の遅延効果を分析しています。私は他の誰かのプログラムに厳密に従いました。それは avg.temp と max.speed で動作しましたが、avg.ap と avg.hum で「coef/vcov が基底行列と一致しません」というエラーが表示されました。ただし、コードに設定されている変数を変更しただけで、他のコードを変更したことはありません。DLNM は雨天を好まないのではないかという仮説があります。TT どうすればいいのかわからないのですが、助けてもらえますか?
パート 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)
パート 2. エラー コード
attach(cpdlnm)
cb.hum = crossbasis(avg.hum, lag=1 ,
argvar=list(fun="ns",
knots= c(10)),
arglag=list(fun="lin"))
modelhum = glm(pre1 ~ cb.hum +
ns(no,1*1),
family=quasipoisson(), 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)
パート 3. データは次のとおりです。
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 12-OCT 2012 10 990.70 23.60 72.20 6.20 11 4.61 12-NOV 2012 11 991.70 18.00 79.70 6.90 2 993.50 15.30 82.20 6.50 15 3.37 13-MAR 2013 3 989.90 20.20 74.20 8.00 16 2.85 13-APR 2013 4 987.00 21.50 78.50 7.70 17 4.38 13-May 2013 5 9830 25.60 79.20 6.80 6.45 13-JUL 2013 7 981.30 28.00 77.50 7.10 20 6.95 13-AUG 2013 8 980.50 27.90 78.20 78.20 7.90 21 6.51 13-SEP 2013 9 985.90 25.40 77.60 6.00 22 8.16 13-OCT 2013 994.50 18.70 72.30 6.20 24 6.18 13-DEC 2013 12 997.30 11.70 67.20 5.30 25 5.69 14-Jan 2014 1 996.70 12.70 70.30 6.00 26 6.44 14-FEB 2014 2 993.00 12.10 76.90 6.40 27 4.16 14-Mar 2014 2014年4月4 987.60 22.60 82.40 6.70 29 3.96 14-May 2014 5 983.60 25.70 78.80 7.70 30 4.72 14-JUN 2014 6 979.20 27.70 81.40 7.90 31 5.21 14-JUL 2014 7.50 33 6.74 14-SEP 2014 9 984.70 27.10 77.70 8.50 34 4.80 14-OCT 2014 10 991.20 23.90 73.10 5.90 35 4.31 14-NOV 2014 11 993.30 18.60 79.60 6.20 36 4.35 14-DEC14 12 998.70 12.30 67.30 5.95 2015 1 996.70 13.30 76.30 6.20 38 4.63 15-FEB 2015 2 993.50 15.50 78.30 6.50 39 4.00 15-MAR 2015 3 991.70 17.70 83.40 6.30 40 4.16 15-APR 2015 4 988.40 22.80 70.20 7.30 4147 -Jun 2015 6 980.90 28.20 81.00 7.40 43 5.04 15-JUL 2015 7 980.20 27.30 79.40 6.70 44 5.79 15-AUG 2015 8 982.40 27.60 80.10 6.50 45 5.28 15-SEP 2015 9 986.00 84.60 6.60 78.30 6.90 47 4.13 15-Nov 2015 11 993.50 19.40 85.30 6.90 48 3.30 15-dec 2015 12 997.80 13.00 80.90 5.70 49 5.30 16-Jan 2016 1 996.00 11.80 82.30 6.40 50 4.5766 16-MAR 2016 3 991.70 17.00 78.90 7.00 52 4.01 16-APR 2016 4 984.60 23.40 80.90 9.80 53 4.90 16-May 2016 5 983.80 25.50 78.70 8.30 54 3.75 16-JUN 2016 6 981.70 28.20 78.20 78.2070 7.2070 7.2070 7.70 981.10 28.90 77.60 7.60 56 3.25 16-AUG 2016 8 979.00 28.00 79.80 8.70 57 2.93 16-SEP 2016 9 984.30 26.60 75.20 6.40 58 2.93 16-OCT 2016 10 987.90 24.40 16-DEC 2016 12 995.70 15.40 71.70 6.80 61 3.10 17-JAN 2017 1 994.70 14.50 79.20 6.50 62 3.75 17-FEB 2017 2 994.80 14.70 71.50 8.30 63 3.49 17-MAR 2017 21.90 76.70 7.80 65 3.69 17-May 2017 5 985.00 24.80 77.50 10.00 66 3.76 17-JUN 2017 6 980.20 26.90 84.80 8.50 67 2.69 17-JUL 2017 2017 9 984.20 27.60 81.50 7.10 70 2.46 17-OCT 2017 10 990.00 22.80 75.90 7.90 71 2.08 17-NOV 2017 11 993.00 17.80 79.50 79.00 72 2.32 17-dec 2017 12 996.90 13.30 74 3.29 18-FEB 2018 2 992.90 13.40 68.70 7.20 75 3.03 18-MAR 2018 3 988.30 19.20 78.20 9.10 76 2.30 18-APR 2018 4 986.50 21.80 77.30 8.70 77 1.75 6 978.30 26.90 81.60 9.00 79 2.79 18-JUL 2018 7 976.80 27.90 82.10 9.20 80 2.32 18-AUG 2018 8 976.40 27.50 83.40 9.60 81 1.88 18-sep 2018 9 983.50 26.10 80.10 -NOV 2018 11 991.50 18.20 80.30 7.10 84 1.78 18-DEC 2018 12 994.50 13.00 84.00 7.80 85 2.77 19-Jan 2019 1 995.20 11.70 84.50 70 86 4.60 19-FEB 2019 2 990.50 13.70 84.80 8.10 87.7087.7087 2.3. 85.90 9.90 88 2.07 19-APR 2019 4 983.60 23.10 84.80 9.80 89 2.97 19-May 2019 5 981.80 24.30 830 7.70 90 90 2.48 19-Jun 2019 6 977.80 27.50 84.8006 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 7.80.0