ncl多元线性回归函数有问题

dannybear

本人最近在用NCL的多元线性回归的函数，
使用的函数是reg_muiltin, 想得到的方程式y=a1+a2*x+a3*x2
其中x（115）为：，  x2为7.88，就是x2也是一个一维的，115个数组，这115个数都是一样的，和x的大小是一样的。
（0）27.73594
(1)        27.73261
(2)        27.57645
(3)        27.73645
(4)        28.38845
(5)        26.94495
(6)        23.35978
(7)        26.57112
(8)        28.33461
(9)        28.96611
(10)        29.84778
(11)        30.67094
(12)        30.90597
(13)        31.32095
(14)        31.07161
(15)        29.37194
(16)        29.00278
(17)        29.26678
(18)        29.80528
(19)        28.56278
(20)        29.65061
(21)        30.74845
(22)        30.30695
(23)        30.28728
(24)        29.60061
(25)        29.75545
(26)        29.58945
(27)        30.42128
(28)        29.68212
(29)        27.94045
(30)        29.42145
(31)        29.22845
(32)        29.91245
(33)        30.48745
(34)        30.94545
(35)        30.58611
(36)        28.98595
(37)        28.92761
(38)        29.04095
(39)        30.28245
(40)        31.09812
(41)        31.17928
(42)        31.28578
(43)        30.71078
(44)        31.71711
(45)        31.21995
(46)        30.83878
(47)        30.73495
(48)        30.57595
(49)        31.00161
(50)        31.16412
(51)        31.02911
(52)        30.57578
(53)        30.21995
(54)        29.23128
(55)        29.50828
(56)        29.46178
(57)        29.01078
(58)        29.48545
(59)        29.50144
(60)        29.74495
(61)        30.68895
(62)        31.44295
(63)        31.00178
(64)        29.98395
(65)        30.83937
(66)        31.34086
(67)        31.93803
(68)        31.7457
(69)        31.3657
(70)        31.74603
(71)        31.3452
(72)        31.40703
(73)        31.06687
(74)        31.3292
(75)        31.04634
(76)        31.49067
(77)        31.02901
(78)        29.01034
(79)        25.1735
(80)        24.81867
(81)        26.10946
(82)        26.91317
(83)        26.44784
(84)        27.81934
(85)        29.84467
(86)        28.61417
(87)        28.43567
(88)        28.93383
(89)        27.6055
(90)        28.31413
(91)        27.85563
(92)        26.21546
(93)        26.24863
(94)        27.11896
(95)        26.31351
(96)        28.10667
(97)        28.10851
(98)        28.47334
(99)        28.01618
(100)        27.56419
(101)        25.50569
(102)        24.19769
(103)        24.85003
(104)        25.45619
(105)        26.47988
(106)        27.15571
(107)        28.62055
(108)        28.57805
(109)        29.47505
(110)        28.7918
(111)        28.90196
(112)        28.86663
(113)        28.59463
(114)        27.75747
但是最后求解出来的，系数是这样的：
Variable: beta1
Type: float
Total Size: 12 bytes
            3 values
Number of Dimensions: 1
Dimensions and sizes:        [3]
Coordinates:
Number Of Attributes: 2
  _FillValue :        9.96921e+36
  constant :        -4.222121e+12
(0)        2.81475e+14
(1)         0
(2)        -3.518437e+13

beta1是系数，分别代表a1,a2,a3，很好奇为什么计算出来是这样的，请大家帮帮忙，是NCL函数的问题吗

Trinity

你还没有看仔细，用reg_multlin函数,必须单独加入x0 =1 "where the X0 is a dummy variable equal to one (1.0) for every observation."建议你用reg_multlin_stats函数,就没有这个问题了. 另外你的的x2好像有点问题啊,什么叫"x2为7.88,一个一维的，115个数组，这115个数都是一样的，和x的大小是一样的"

dannybear

Trinity 发表于 2016-12-24 11:34
你还没有看仔细，用reg_multlin函数,必须单独加入x0 =1 "where the X0 is a dummy variable equal to one ( ...

意思就是，x2是115个数，每个数都等于7.8，我之前使用了reg_multin_stats这个函数，因为我的资料有缺测，就没有办法使用这个函数，当我把缺测值去掉之后，得到的结果和reg_multin是一样的

zju3060601074

谢谢楼主 学习一下