[R] loop

Luis Silva lm.silva at sapo.pt
Mon Apr 21 12:03:05 CEST 2003


I tried the tips from Sundar but I  still have problems. I want 
both values and vectors, and the latter will be used to project 
my original data. The problem is that with those command lines 
I can't get them as a matrix and when I try as.matrix it 
doesn't work. With only vectors

> b<-sapply(my.list, function(x) eigen(x)$vectors)
> attributes(b)
$dim
[1] 65536    11

$dimnames
$dimnames[[1]]
NULL

$dimnames[[2]]
 
[1] "0"   "0.1" "0.2" "0.3" "0.4" "0.5" "0.6" "0.7" "0.8" "0.9" 
"1"  


> b[1]
[1] 0.02533057+0i

thanks
luis

} 
} 
} Luis Silva wrote:
} > Dear helpers
} > 
} > I have this problem. I want to make a linear
} combination a*A+(1-
} > a)*B where A and B are matrices. I want that a be
} incremented 
} > from 0 to 1 by 0.1 so I made a loop with for. The
} problem is 
} > that I want to keep the result in an object or list
} or 
} > something like that and then apply eigen
} decomposition to all 
} > of the resulting matrices
} > 
} > sapply(my.list,eigen)
} > 
} > the problem is that I don't know how to build that
} list in the 
} > loop. I tried several things but it doesn't work
} (in Matlab i 
} > can do it)
} > 
} 
} a <- seq(0, 1, 0.1)
} my.list <- lapply(a, function(a, A, B) a*A + (a-1)*B,
} A=A, B=B)
} names(my.list) <- as.character(a)
} sapply(my.list, function(x) eigen(x)$value)
} 
} Regards,
} Sundar
} 
} 

--


http://adsl.sapo.pt



More information about the R-help mailing list