R help - Jan 2007 - for loop problem

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)
}

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 

-- 
View this message in context:
http://www.nabble.com/for-loop-problem-tf3047849.html#a8472217
Sent from the R help mailing list archive at Nabble.com.

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 
>