[R] how to read this matrix into R

Spencer Graves spencer.graves at pdf.com
Sat Nov 6 02:19:45 CET 2004


      How about the following: 

 > a <- scan("clipboard")
Read 45 items
 > n2 <- length(a)
 > n <- (-1+sqrt(1+8*n2))/2
 > A <- array(NA, dim=c(n,n))
 > A[lower.tri(A,diag=TRUE)] <- a
 > A
        [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8] [,9]
 [1,]  1.000    NA    NA    NA    NA    NA    NA    NA   NA
 [2,]  0.591 1.000    NA    NA    NA    NA    NA    NA   NA
 [3,]  1.000 0.573 0.517    NA    NA    NA    NA    NA   NA
 [4,]  0.356 0.408 0.424 0.320    NA    NA    NA    NA   NA
 [5,]  0.350 0.382 0.303 0.401 0.512    NA    NA    NA   NA
 [6,]  1.000 0.062 1.000 0.479 0.346 1.000    NA    NA   NA
 [7,] -0.098 1.000 0.400 1.000 0.308 0.519 0.167    NA   NA
 [8,]  0.072 0.156 0.414 0.282 0.463 0.484 0.455 0.574   NA
 [9,]  0.380 0.232 0.611 0.375 0.605 0.467 0.311 0.557    1

      hope this helps.  spencer graves

rongguiwong wrote:

>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.
>
>
>  1.000                                                         
>  0.591  1.000                                                 
>  0.356  0.350  1.000                                          
> -0.098  0.072  0.380  1.000                                   
>  0.573  0.408  0.382  0.062  1.000                            
>  0.156  0.232  0.517  0.424  0.303  1.000                     
>  0.400  0.414  0.611  0.320  0.401  0.479  1.000              
>  0.282  0.375  0.512  0.346  0.308  0.463  0.605  1.000       
>  0.519  0.484  0.467  0.167  0.455  0.311  0.574  0.557  1.000
>
>______________________________________________
>R-help at 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
>  
>

-- 
Spencer Graves, PhD, Senior Development Engineer
O:  (408)938-4420;  mobile:  (408)655-4567




More information about the R-help mailing list