0

我一直在研究R非线性模型，比如我们：

Y = Alpha1*time + Alpha2*sin(2*pi*time/Alpha3) + Alpha4*(-1)^time

我会确定一个伯努利变量是否会影响Alpha1*time和截距。这样的伯努利变量可以是：

varia<-rep("BEFORE","AFTER"),each=30)

有

Y<- -2.5+rnorm(60)+2*sin(2*pi*time/8)+2.5*(-1)^time
time<-seq(1,60)

我发现它nls可以拟合这个模型，但没有这个分类变量的影响，我还发现这个nlme包可以估计变量的影响，而不是非线性模型的项。我的问题是：这些软件包中的哪一个可以帮助我？而且，我怎么能在代码中添加这个伯努利变量？问候和感谢。

4

1 回答 1

1

以下.lin4将是AFTER增加的金额Y：

# test data
set.seed(123)
time<-seq(1,60)
Y<- -2.5+rnorm(60)+2*sin(2*pi*time/8)+2.5*(-1)^time
varia <- rep(c('BEFORE', 'AFTER'), each = 30)

fm <- nls(Y ~ cbind(1, sin(2 * pi * time/Alpha3), (-1)^time, varia == "AFTER"), 
 start = list(Alpha3 = 8), alg = "plinear")

summary(fm)

给予：

Formula: Y ~ cbind(1, sin(2 * pi * time/Alpha3), (-1)^time, (varia == 
    "AFTER"))

Parameters:
       Estimate Std. Error t value Pr(>|t|)    
Alpha3  8.03367    0.02583 310.984  < 2e-16 ***
.lin1  -2.55242    0.16740 -15.247  < 2e-16 ***
.lin2   1.88528    0.16726  11.272 6.45e-16 ***
.lin3   2.53734    0.11828  21.453  < 2e-16 ***
.lin4   0.23992    0.23673   1.013    0.315    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.9159 on 55 degrees of freedom

Number of iterations to convergence: 3 
Achieved convergence tolerance: 6.226e-07

因此，在这种情况下，除了最后一个参数外，所有参数都很重要，正如预期的那样。

REVISED重新阅读问题后修改。

于 2014-04-28T16:09:20.153 回答