how to read this matrix into R
rongguiwong <0034058 at fudan.edu.cn> writes:
the following the the lower.tri matrix in a file named luxry.car
and i want to read it in R as a lower.tri matrix.how can i do?
i have try to use help.search("read"),but no result what i want.
Here's one way:
x <- scan()
1: 1.000 2: 0.591 1.000 4: 0.356 0.350 1.000 7: -0.098 0.072 0.380 1.000 11: 0.573 0.408 0.382 0.062 1.000 16: 0.156 0.232 0.517 0.424 0.303 1.000 22: 0.400 0.414 0.611 0.320 0.401 0.479 1.000 29: 0.282 0.375 0.512 0.346 0.308 0.463 0.605 1.000 37: 0.519 0.484 0.467 0.167 0.455 0.311 0.574 0.557 1.000 46: Read 45 items
m <- matrix(0,9,9) m[upper.tri(m,diag=TRUE)] <- x # Yes, upper! m <- t(m) m[upper.tri(m,diag=TRUE)] <- x
(Assuming that you really want a symmetric matrix. If you want just the lower triangle, simply omit the last line.)
O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk) FAX: (+45) 35327907