Displaying 1 result from an estimated 1 matches for "sig_xx2".
Did you mean:
sig_xx
2006 Oct 03
1
do.call with Vectorial Argument
...*myargs[1])-1)/(2*lam1)-2/(lam1+lam2)*(exp(lam1*myargs[1]+lam2*myargs[1])-1)
+ (exp(2*lam2*myargs[1])-1)/(2*lam2) )
}
# Now I use do.call
newtime<-seq(1,5,len=1001)
mypar<-c(q,lam1,lam2)
sig_xx<-do.call("sigma_pos_old",c(list(t=newtime),mypar))
# Now this line does not work; sig_xx2 is not a vector equal to sig_xx
sig_xx2<-do.call("sigma_pos",c(list(myargs =newtime)))
So I am making some mistake with the do.call command, since I get a
single value rather than a vector.
Can anyone tell me what I am doing wrong? I tried several variations
of the syntax, but non...