[R] How to calculate the spatial correlation of several files?

Rui Barradas ruipbarradas at sapo.pt
Mon Dec 3 17:23:22 CET 2012


Hello,

Sorry, made a mistake.
Em 03-12-2012 16:12, Rui Barradas escreveu:
> Hello,
>
> Inline.
> Em 03-12-2012 15:15, Jonsson escreveu:
>>      dir1 <- list.files("C:\\Users\\aalyaari\\Desktop\\cor", "*.bin",
>> full.names = TRUE)
>>      dir2 <- list.files("C:\\Users\\aalyaari\\Desktop\\cor2", "*.bin",
>> full.names = TRUE)
>>      results <- list()
>>          for (.files in dir1){   # read in the 365 files as a vector of
>> numbers for dir1
>>          file1 <- do.call(rbind,(lapply(.files, readBin  , integer() 
>> , size =
>> 2 ,
>>                                      n = 360 * 720 , signed = T)))    }
>>          for (.files in dir2){   # read in the 365 files as a vector of
>> numbers for dir2
>>          file2<- do.call(rbind,(lapply(.files, readBin  , integer() , 
>> size =
>> 2 ,
>>                          n = 360 * 720 , signed = T)))   }
>>         # Now each file  in both directories is a vector.
>
> Are you sure? Shouldn't your file reading routines be
>
> do.call(rbind, lapply(dir1, readBin, integer(), size = 2, n = 360 * 
> 720, signed = T))
> do.call(rbind, lapply(dir2, readBin, integer(), size = 2, n = 360 * 
> 720, signed = T))
>
>
> to lapply readBin to each file in dir1/dir2?
>
> Anyway, to correlate the first row in file1 to the first row in file2, 
> etc, try

Here. it should be "column" not "row".

Rui Barradas
>
> for (.f in seq_along(dir1)){
>     results[[.f]]<- cor(file1[, .f] ,file2[, .f])
> }
>
>
> Hope this helps,
>
> Rui Barradas
>
>>   I am not sure how
>> to tell R to correlate the first column in dir1 to the correspond column
>> from dir2. we will finally get only one spatial correlation map.
>> I tried to this:
>>   # calculate the  correlation so we will get a correlation map
>>      for (.files in seq_along(dir1)){
>>          results[[length(results) + 1L]]<- cor(file1 ,file2)
>>          }
>> I got error:Error in cor(file1, file2) : allocMatrix: too many elements
>> specified`
>>
>>
>>
>>
>>
>> -- 
>> View this message in context: 
>> http://r.789695.n4.nabble.com/How-to-calculate-the-spatial-correlation-of-several-files-tp4651888.html
>> Sent from the R help mailing list archive at Nabble.com.
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide 
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.




More information about the R-help mailing list