[R] indexed expression

Uwe Ligges ligges at statistik.tu-dortmund.de
Sat Dec 27 17:17:43 CET 2008



bernardo lagos alvarez wrote:
> Hello expeRts,
> 
> I need generate symbolize the autocovariances matrix of a Gaussian
> ARMA(1,1), for derivate it and evaluate.
> I try this codes, but whitout sucess
> 
> vacv<-NULL
>  vacv[1]<-expression(1-2*phi*theta-theta^2)
>  vacv[2]<-(1-phi*theta)*(phi-theta)
>  vacv[3:n]<-acv[2]*(phi^(1:(n-2)))
>  facv<-list()
> for(i in 1:2)
> facv[[i]]<-deriv(y~vacv[i],c("phi","theta"),function(phi,theta){})
> Erro en deriv.formula(y ~ vacv[i], c("phi", "theta"), function(phi, theta) { :
>         Función '`[`' no está en la tabla de derivadas


1. Please try to provide English error messages.
2. If you want deriv to work on the expressions for vacv[], then please 
do it rather than ttrying to work on already evaluated expressions as 
you do.

One exmaple that works would be to say:
  vacv <- NULL
  vacv[1] <- expression(1-2*phi*theta-theta^2)
  ...
and later on:
   deriv(vacv[i], c("phi", "theta"), function(phi,theta){})

of yourse, you need to evaluate vacv later on if you really want to 
calculate the values....

Uwe Ligges


> 
> dfg<-for(n in 1:2) as.formula(sprintf("y~(1 - phi * theta) * (phi -
> theta)*phi^ %d",n))
>  dfg
> y ~ (1 - phi * theta) * (phi - theta) * phi^2
> deriv(dfg[[1]],phi)
> Erro en deriv.default(expr, namevec, function.arg, tag, hessian) :
>         nombres de variable inválidos
> 
> Give me a help, please.
> 
> Bernardo.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list