I have a simulated matrix of dates that I generated from a probability 
function. Each column represents a single iteration. 

I would like to bin each run _separately_ by decades and dump them into a new 
matrix where each column is the length of all decades a single run with the 
number dates binned by decade. 

I have successfully done this for a single vector of dates, but not for a 
matrix: 

"dates" is a vector of observed data representing when certain trees 
established in a population 

#-----find min and max decade -----# 

        minDecade <- min(dates) 
        maxDecade <- max(dates) 

#-----create vector of decades -----# 
        
        allDecades <- seq(minDecade, 2001, by=10) 

#-----make empty vector of same length as decade vector-----# 
        
        bin.vec <- rep(0,length(allDecades)) 

#-----populate bin.vec (empty vector) with the number of trees in each 
decade-----# 

        for (i in 1:length(allDecades)) { 
                        
                        bin.vec[i] <- length(which(dates==allDecades[i])) 
                } 
                        

bin.vec : [1]    0  0  0  0  0  0  0  0  0  0  1  1  1  0  1  2  0  1 
               [19]  3  0  1  3  8  5  9  8  5  5  4 10  3  6  9 17 32 37 
               [37] 35 25 31 41 41 44 45 40 50 43 59 42 46 28 16 18 20 16 
               [55] 11  4  7  1 


My matrix looks like this (it actually had 835 rows, I used head (x) just to 
demonstrate). 

head(bin.mat) 
     [,1]     [,2]     [,3]   [,4]   [,5]    [,6]    [,7]   [,8]    [,9]     
[,10] 
[1,] 1831 1811 1841 1881 1851 1871 1921 1821 1781  1561 
[2,] 1851 1931 1821 1701 1841 1961 1941 1931 1891  1841 
[3,] 1751 1861 1861 1751 1841 1841 1771 1971 1811  1871 
[4,] 1831 1871 1741 1881 1871 1771 1821 1901 1901  1851 
[5,] 1681 1861 1871 1811 1711 1931 1891 1771 1811  1821 
[6,] 1931 1841 1841 1861 1831 1881 1601 1861 1891  1891 

Each column is a separate run (runs <- 10 ) .  How can I bin each column into 
decades separately? 

I'll bet this is super easy, but my R-skills are seriously limited!!! 

Thanks for any help! 
~Jeff
        [[alternative HTML version deleted]]

______________________________________________
R-help@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.

Reply via email to