Displaying 3 results from an estimated 3 matches for "lambda4".
Did you mean:
lambda
2008 Sep 11
0
Loop for the convergence of shape parameter
...a2*(14^a-10^a+14^a-14^a))
-exp(-lambda1*(10^a)-lambda2*(14^a-10^a+18^a-14^a)))
logl<-log(s1_f1*s2_f2)
return(-logl)}
optim(1,L)
alpha<-optim(1,L)$par
}
[II] for 1)-3) (making errors)
alpha[1]<-1
beta0=rep(0,10)
beta1=rep(0,10)
lambda1=rep(0,10)
lambda2=rep(0,10)
lambda3=rep(0,10)
lambda4=rep(0,10)
verpi=rep(0,10)
s1_f1=rep(0,10)
s2_f2=rep(0,10)
a=rep(0,10)
logl=rep(0,10)
beta0[1]=0
beta1[1]=0
lambda1[1]=0
lambda2[1]=0
lambda3[1]=0
lambda4[1]=0
verpi[1]=0
s1_f1[1]=0
s2_f2[1]=0
a[1]=1
logl[1]=0
for(i in 2:11){
verpi[i]<-c(5^alpha[i-1],10^alpha[i-1]-5^alpha[i-1],14^alpha[i-1]-10^al...
2012 Mar 11
1
CRAN (and crantastic) updates this week
CRAN (and crantastic) updates this week
New packages
------------
* EffectStars (1.0)
Maintainer: Unknown
Author(s): Gunther Schauberger
License: GPL-2
http://crantastic.org/packages/EffectStars
The package provides functions to visualize regression models with
categorical response. The effects of the covariates are plotted with
star plots in order to allow for an optical
2011 Feb 04
0
MSBVAR and hc.forecast
....1,1.1,2.1,3.1,100.1,1.1,2.1,3.1,100.1)
b = c(1.1,1.1,1.1,100.1,1.1,1.1,1.1,100.1,1.1,1.1,1.,100.1,1.1,1.1,1.1,100.1,1.1,1.1,1.1,100.1,1.1,1.1,1.,100.1)
K <-ts(cbind(a,b),start=c(1.1,1.1), names=c("a","b"))
fit.bvar <- szbvar(K, p = 1, lambda0=0.6, lambda=0.1, lambda3=2,
lambda4=0.25,lambda5=0,mu5=0, mu6=0, prior=0)
y <- matrix(c(rep(0,8),c(1.10,1.10,1.10,100.1,1.10,1.10,1.0,100.10)),ncol=2)
h <- hc.forecast(fit.bvar,y,nsteps=8,burnin=300, gibbs=500,exog=NULL)