Re: [R] for loop problem

From: Uwe Ligges <ligges_at_statistik.uni-dortmund.de>
Date: Sun 21 Jan 2007 - 14:10:22 GMT

aat wrote:
> Hello R users,
>
> A beginners question which I could not find the answer to in earler posts.
>
> My thought process:
> Here "z" is a 119 x 15 data matrix
> Step 1: start at column one, bind every column with column 1
> Step2: use the new matrix, "test", in the fitCopula package
> Step3: store each result in myfit, bind each result to "answer"
> Step4: return "answer"
>
>
> copula_est <- function(z)
> {
> for(i in 1:length(z[1,]))
> {
> my.cop <- normalCopula(param = 0.5, dim = 2)
> test <- cbind(z[,1],z[,i])
> myfit[i] <- fitCopula(test,my.cop, start=0.3)
> }
> answer <- cbind(myfit[i])
> return(answer)
> }

The example is not reproducible for us, since we do not have z.

I'd try to rewrite it as follows, without having tried anything:

my.cop <- normalCopula(param = 0.5, dim = 2) answer <- apply(z[,-1], 2,

    function(x) fitCopula(cbind(z[,1], x), my.cop, start=0.3),     my.cop = my.cop)

Uwe Ligges

> Errors received:
> Error: object "test" not found
>
> Could my syntax be incorrect, or is it a deeper faulty logic error.
> Thank you for your help, it is much appreciated.
>
> aat
>



R-help@stat.math.ethz.ch 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. Received on Mon Jan 22 01:16:19 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sun 21 Jan 2007 - 15:30:22 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.