[R] binning dates by decade for simulated data
jtgarcia
jtgarcia20 at gmail.com
Thu Mar 8 22:47:39 CET 2012
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
--
View this message in context: http://r.789695.n4.nabble.com/binning-dates-by-decade-for-simulated-data-tp4457766p4457766.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list