关于非线性哑变量询问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

在做非线性哑变量模型时，遇到一点问题：
代码如下（数据在附件中）

1. rd<-read.csv("data1.csv")
   
2. attach(rd)
   
3. head(rd)
   
4. rd$region<- as.factor(rd$region)
   
5. total <- nls(TOTAL~a*DBH^b+factor(region),
   
6. data=rd,start=list(a=0.1,b=1))

复制代码

出现：
Error in numericDeriv(form[[3L]], names(ind), env) :
  Missing value or an infinity produced when evaluating the model
In addition: Warning messages:
1: In Ops.factor(a * DBH^b, factor(region)) : + not meaningful for factors
2: In Ops.factor(a * DBH^b, factor(region)) : + not meaningful for factors

哪位高手指点一下，应该是稍微修改一下就可以，但是我不知道怎么修改

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：非线性 哑变量 模型

nlm是做非线性优化的函数，跟回归没关系啊。你再仔细读读帮助文件吧。

hugebear 发表于 2014-4-8 12:38
nlm是做非线性优化的函数，跟回归没关系啊。你再仔细读读帮助文件吧。

是我自己写错了，谢谢，现在改成nls,还是做不出来，你能帮一下吗？

1. total <- nls(TOTAL~a*DBH^b+as.numeric(region),
   
2. data=rd,start=list(a=0.1,b=1))

复制代码

yywan0913 发表于 2014-4-8 16:36

谢谢指点
结果中每个区域的参数值还是出不来呀？
结果是这样：
Formula: TOTAL ~ a * DBH^b + as.numeric(region)

Parameters:
  Estimate Std. Error t value Pr(>|t|)   
a  0.12513    0.01834   6.823 6.83e-11 ***
b  2.40360    0.04759  50.502  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 18.87 on 247 degrees of freedom

Number of iterations to convergence: 12
Achieved convergence tolerance: 9.828e-08

1. total <-sapply(unique(as.numeric(rd$region)),function(x) nls(TOTAL~a*DBH^b+x,data=rd,start=list(a=0.1,b=1)))
   
2. colnames(total)=unique(region)

复制代码

难道要这样？
region也当做变量？