Skip to content

matrix indexing in 'for' loop?

7 messages · govindas at msu.edu, Phil Spector, Bert Gunter +1 more

Phil Spector
# 
To make your loop work, you need to learn about the get function.
I'm not going to give you the details because there are better
approaches available.
First, let's make some data that will give values which can 
be verified.  (All the correlations of the data you created
are exactly equal to 1.)  And to make the code readable, I'll
omit the ts.m prefix.
Since cor.test won't automatically do the tests for all columns
of a matrix, I'll write a little helper function:
Now tests is a list of length 2, with a list of the output from
cor.test for the five columns of the each matrix with pe[,2]
(Notice that in your program you made no provision to store 
the results anywhere.)

Suppose you want the correlations:
[,1]       [,2]
cor  0.12723615  0.1342751
cor  0.07067819  0.6228158
cor -0.28761533  0.6218661
cor  0.83731828 -0.9602551
cor -0.36050836  0.1170035

The probabilities for the tests can be found similarly:
[,1]       [,2]
[1,] 0.8727638 0.86572490
[2,] 0.9293218 0.37718416
[3,] 0.7123847 0.37813388
[4,] 0.1626817 0.03974489
[5,] 0.6394916 0.88299648

(Take a look at the "Values" section in the help file for cor.test
to get the names of other quantities of interest.)

The main advantage to this approach is that if you add more matrices
to the allmats list, the other steps automaticall take it into account.

Hope this helps.
 					- Phil Spector
 					 Statistical Computing Facility
 					 Department of Statistics
 					 UC Berkeley
 					 spector at stat.berkeley.edu
On Tue, 21 Dec 2010, govindas at msu.edu wrote:

            

      
      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Liviu Andronic
    

    

      
      #
    

  


  

      
On Wed, Dec 22, 2010 at 2:57 AM, Phil Spector <spector at stat.berkeley.edu> wrote:

            

      
      
      
    

        
Hmisc already provides a function for doing this, rcorr(). Try to
compute correlations using Rcmdr for an example.

Regards
Liviu

            

      
      
      
    

        

      
      
    

        

  
    
  


  


      

    

    

      
      

        
      

    

    

      
      

        


  

        

      Bert Gunter
    

    

      
      #
    

  


  

      
Have you consulted R's extensive documentation? -- in particular, "An
Introduction to R," which would seem like an obvious place for R
newbies to start. If you had done so, you would have found your
question answered there in section 6.1 on lists.

-- Bert Gunter

        
On Wed, Dec 22, 2010 at 9:39 AM,  <govindas at msu.edu> wrote:

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        

      
      
    

        

  
    

      
      
    

  


  


      

    

    

      
      

        


  

        

      Liviu Andronic
    

    

      
      #
    

  


  

      
On Wed, Dec 22, 2010 at 6:39 PM,  <govindas at msu.edu> wrote:

            

      
      
      
    

        
Is this what you need?

            

      
      
      
    

        
[,1]        [,2]        [,3]       [,4]      [,5]
[1,]  0.2523449 -1.26606999  0.27594223 -0.5158524 -0.599899
[2,] -1.9155409  0.01081565  0.44497873  0.1819517  1.187899
[3,]  0.7624774  0.31296468 -0.02187076  0.9389091  0.865527
[4,] -0.8082626 -1.44330148  1.32975075  0.1788399  1.204917

            

      
      
      
    

        
[1]  0.2523449 -1.9155409  0.7624774 -0.8082626

Liviu