I have a problem with R. I try to compute the confidence interval for my
df. When I want to create the plot I have this problem: Error in
xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differ.
I try this code:
library(dplR)
df.rwi <- detrend(rwl = df, method = "Spline",nyrs=NULL)
write.table(df.rwi,file="rwi.txt",quote=FALSE,row.names=TRUE)
df.crn <- chron(df.rwi, prefix = "RGE",prewhiten = TRUE)
write.table(df.crn,file="crn.txt",quote=FALSE,row.names=TRUE)
##confidence interval
library(boot)
ci.func <- function(y) {
y <- y[!is.na(y)]
mean.fun <- function(y,i)
{
m <- mean(y[i])
n <- length(i)
v <- (n-1)*var(y[i])/n^2
c(m,v)
}
y.boot <- boot (y, mean.fun, R=1000)
y.ci=boot.ci (y.boot, conf=0.99, type = 'basic')
y.ci=y.ci$basic[c(4,5)]
names(y.ci)=c('lower','upper')
y.ci
}
df.trunc<-subset(df.crn,samp.depth>5)
df.rwi.trunc <- subset(df.rwi, df.crn$samp.depth > 5)
df.ci <- data.frame(t(apply(df.rwi.trunc,1,ci.func)))
base<-cbind(df.trunc,df.ci)
op=par(no.readonly=TRUE) # save par to reset
par(mar=c(3,2,2,3)+0.1,mgp=c(1.25,0.25,0),tcl=0.25,tck=0.0125)
yr.vec=as.numeric(rownames(df.trunc))
samp.depth=df.trunc[,2]
crn2=df.trunc[,-2]
plot(yr.vec,crn2,type='n',xlab='Years',ylab='RWT',
main='Gemenele')
axis(1,at=seq(1900,2012,by=4))
polygon(c(yr.vec,rev(yr.vec)), c(df.ci$lower,rev(df.ci#upper)),
col='grey80',border=NA)
lines(yr.vec,crn2,col='black')
abline(h=1)
par(new=TRUE)
plot(yr.vec,samp.depth,type='1',col='black',xlab='',ylab='',axes=FALSE,
lty=2)
axis(4,at=pretty(samp.depth))
mtext('Sample Depth',side=4,line=1.25)
legend("topleft", c('Chronology', '99%
CI', 'Samples'),
col=c('black','grey50','black'),lwd=rep(2,3))
My data is like this:
structure(list(RGE01 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1.47, 3.11,
2.87, 2.69, 1.98, 2.22, 1.93, 1, 0.61, 0.57, 0.55, 1, 1.41, 0.98,
2.32, 1.86, 1.81, 1.48, 2.18, 2.29, 2.31, 2.28, 1.72, 1.8, 1.71,
1.59, 1.11, 1.24, 1.1, 1.83, 0.43, 0.98, 0.71, 1.09, 0.79, 1.11,
0.55, 0.44, 0.36, 0.5, 0.5, 0.51, 0.48, 0.42, 0.44, 0.44), RGE02 = c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
1.97, 0.94, 1.95, 1.95, 1.92, 1.14, 1.84, 1.45, 1.63, 1.81, 1.24,
1.29, 1.31, 1.41, 1.2, 2.04, 2.23, 2.44, 2.11, 2.24, 0.93, 1.89,
1.93, 1.58, 2.36, 2.07, 1.79, 1.4, 1.05, 2.14, 2.14, 2.65, 3.47,
2.93, 1.52), RGE03 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 1.62,
0.83, 0.81, 0.84, 0.83, 1.06, 0.13, 0.72, 0.98, 0.54, 0.8, 0.6,
0.44, 0.59, 0.78, 0.62, 0.7, 0.98, 0.68, 0.58, 0.5, 0.45, 0.37,
0.26, 0.44, 0.6, 0.48, 0.55, 0.88, 0.4, 0.4, 0.45, 0.46, 0.52,
0.46, 0.36, 0.43, 0.36, 0.29, 0.36, 0.64, 0.49, 0.92, 0.9, 0.54,
0.7, 0.84, 0.52, 0.81, 0.95, 0.8, 0.64, 0.84, 0.6, 0.78, 0.76,
0.88, 0.74, 0.9), RGE04 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, 2.29, 2.4, 1.74, 1.56, 2.06,
1.11, 1.63, 1.83, 1.41, 1.37, 1.34, 1.18, 0.92, 1.64, 1.21, 1.5,
1.12, 0.7, 1.05, 1.19, 1.17, 0.68, 0.62, 0.97, 0.71, 0.73, 0.92,
0.64, 0.48, 0.45, 0.96, 0.84, 0.97, 0.33, 0.26, 0.31, 0.6, 0.71,
0.5, 0.46, 0.24, 0.26, 0.65, 0.85, 0.9, 0.27, 0.3, 0.86, 0.82,
0.6), RGE05 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 2.36,
2.65, 1.78, 1.29, 1.08, 1.9, 2.33, 1.72, 1.45, 0.32, 1.42, 1.06,
1.46, 1.61, 1.24, 1.56, 1.58, 1.66, 1.49, 0.55, 1.58, 1.96, 1.26,
0.82, 0.92, 0.84, 0.72, 0.72, 0.68, 0.98, 0.5, 0.84, 0.4, 0.3,
1.22, 0.74, 0.97, 0.52, 0.88, 0.63, 0.58, 0.56, 0.62, 0.47, 0.46,
0.72, 1.2, 1.01, 0.88, 0.87, 1.58, 0.82, 0.67, 1.27, 0.96, 0.54,
0.48, 1.14), RGE06 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, 2.37, 2.32, 1.39, 1.2, 0.97, 1.01, 0.83, 0.86, 0.74,
0.65, 1.04, 0.79, 0.52, 0.54, 0.57, 0.45, 0.44, 0.42, 0.47, 0.75,
0.86, 0.61, 0.79, 0.73, 0.53, 0.61, 0.55, 0.63, 0.7, 0.62, 0.64,
0.64, 0.58, 0.5, 0.46, 0.65, 0.76, 0.44, 0.81, 0.88, 0.48, 0.38,
0.4, 0.57, 0.39, 0.45, 0.41, 0.44, 0.6, 0.37, 0.2, 0.23, 0.49,
0.31, 0.24, 0.38, 0.37, 0.4, 0.48, 0.46, 0.39, 0.26, 0.4, 0.38,
0.52, 0.51, 0.78, 0.68, 0.52, 0.51, 0.5, 0.28, 0.27, 0.42, 0.4,
0.41, 0.46, 0.52, 0.53, 0.33, 0.3, 0.48, 0.48, 0.43, 0.34, 0.36,
0.28, 0.28), RGE07 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, 2.82, 1.47, 1.37, 1.06, 1.01,
0.72, 1.3, 1.16, 0.72, 0.88, 1.02, 0.92, 0.8, 0.89, 0.94, 1,
0.78, 0.89, 1.08, 0.83, 0.99, 0.97, 0.82, 0.9, 0.82, 0.7, 0.87,
0.43, 0.65, 0.55, 0.49, 0.59, 0.63, 0.38, 0.31, 0.64, 0.38, 0.48,
0.47, 0.72, 0.53, 0.46, 0.48, 0.56, 0.56, 0.67, 0.5, 0.7, 0.63
), RGE08 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, 1.03, 1.24, 1.2, 0.34, 0.79, 0.79, 0.72, 1.28, 0.62,
0.94, 0.9, 0.74, 0.96, 1.07, 1.18, 1.28, 1.64, 1.18, 0.98, 0.69,
0.64, 0.74, 0.31, 0.53, 0.26, 0.74, 0.94, 0.64, 0.6, 0.36, 0.31,
0.89, 1.05, 1.21, 1.36, 1.52, 1.94, 0.79, 0.68, 0.98, 0.58, 0.64,
1.22, 1.09, 0.82, 0.76, 1.22, 1.16, 0.58, 0.6, 0.74, 0.67, 0.92,
0.96, 0.9, 0.82, 0.9, 0.76, 1.39, 0.75, 0.93, 0.88, 0.65, 1,
1.01, 0.74, 0.87, 0.86, 0.6, 0.96, 1.04, 1.55, 1.2, 0.96, 0.58,
0.54, 0.48, 1.08, 0.85, 0.73, 0.68, 0.3, 0.57, 0.83, 0.43), RGE09 = c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, 2, 1.3, 0.67, 0.72, 0.88, 0.76, 0.79, 0.78,
0.67, 0.42, 0.35, 0.26, 0.49, 0.68, 1.12, 0.85, 0.96, 0.5, 0.98,
0.66, 0.81, 1, 0.72, 0.78, 1.17, 1.22, 0.78, 1.54, 2.01, 0.99,
0.96, 0.7, 0.62, 0.65, 0.53, 0.61, 1.06, 0.79, 1.29, 0.61, 0.23,
0.32, 0.45, 0.5, 0.23, 0.37, 0.4, 0.45, 0.71, 0.56, 0.44, 0.58,
0.54, 0.56, 0.35, 0.52, 0.44, 0.41, 0.31, 0.37, 0.39, 0.24),
RGE10 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, 1.96, 3.51, 2.77, 3.79,
2.78, 1.81, 1.16, 0.58, 1.16, 2.07, 1.45, 2.08, 2.4, 2.42,
2.07, 1.22, 1.73, 0.82, 1.24, 0.34, 0.4, 0.5, 0.66, 0.53,
0.85, 0.91, 0.86, 0.49, 0.86, 0.8, 0.74, 1.12, 1.74, 1.16
), RGE11 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 3.41, 2.09, 2.38,
3.82, 2.11, 1.86, 1.34, 1.42, 1.23, 0.62, 0.6, 0.96, 1.83,
1.21, 1.1, 1, 0.5, 0.41, 0.64, 0.41, 0.64, 0.6, 0.83, 0.84,
0.52, 0.71, 0.9, 0.97, 0.78, 0.83, 0.48, 0.55, 0.77, 0.65,
0.76, 0.47, 0.67, 0.22, 0.68, 0.74, 0.51, 0.86, 0.96, 0.8,
0.68, 0.57, 0.81, 0.78, 0.62, 0.35, 0.34, 0.3, 0.33, 0.25,
0.39, 0.65, 0.34, 0.46, 0.33, 0.43, 0.74, 0.51, 0.58, 0.37
), RGE12 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, 1.3, 2.76, 1.3, 2.15, 1.69, 2.6, 1.42, 0.94, 1.25,
1.06, 1.05, 0.87, 0.98, 1, 0.68, 0.18, 0.63, 0.44, 0.65,
0.72, 0.51, 0.49, 0.45, 1, 0.96, 0.64, 0.72, 0.92, 0.99,
1.38, 0.72, 1.44, 0.91, 1.06, 0.88, 0.82, 0.37, 0.7, 0.79,
0.86, 0.44, 0.7, 0.5, 0.61, 0.44, 0.53, 0.34, 0.56, 0.37,
0.54, 0.67, 1.1, 0.77, 1.01, 0.83, 0.44, 0.85, 0.76, 0.56,
0.71, 0.53, 0.8, 0.41, 0.42, 0.22, 0.19, 0.66, 0.63, 0.6,
0.45, 0.55, 0.53), RGE13 = c(NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, 3.23, 2.09, 2.71, 2.23, 1.39, 1.28, 1.11, 0.75,
1.32, 0.55, 0.37, 0.53, 0.35, 0.57, 0.97, 1.04, 0.65, 0.5,
0.97, 1.42, 1.18, 0.5, 0.66, 0.51, 0.52, 0.56, 1, 0.82, 1.11,
1.36, 1.26, 0.92, 0.94, 1.25, 0.96, 0.94, 0.76, 1.06, 0.96,
1.06, 0.76, 0.86, 0.99, 0.97, 0.44, 0.41, 0.77, 1.15, 0.77,
0.5, 0.6, 0.93, 0.73, 0.66, 0.42, 0.44, 0.56, 0.49, 0.31,
0.84), RGE14 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 3.97, 1.66,
1.95, 2.3, 1.92, 1.74, 0.75, 0.52, 0.77, 1.02, 1.56, 1.56,
1.27, 0.74, 0.99, 0.68, 0.99, 0.81, 0.43, 0.32, 0.64, 0.15,
0.23, 0.53, 0.2, 0.59, 0.16, 0.34, 0.68, 0.8, 0.7, 0.69,
0.85, 0.61, 0.44, 0.76, 0.82, 0.72, 0.99, 1.09, 0.64, 0.58,
0.6, 0.26, 0.57, 0.59, 0.38, 0.55), RGE15 = c(NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, 3.15, 2.62, 0.92, 1.24, 0.7, 0.93, 0.67, 1.54,
1.14, 0.92, 1.73, 1.16, 1.2, 0.82, 0.78, 1.09, 0.32, 0.94,
1.05, 0.62, 0.44, 0.96, 1.24, 0.86, 1.13, 0.79, 0.45, 0.65,
0.98, 0.72, 1.02, 1.13, 0.79, 0.93, 0.82, 0.78, 0.51, 0.72,
0.57, 0.3, 0.47, 0.64, 0.27, 1.12, 0.68, 0.6, 0.22, 0.58,
0.5, 0.3, 0.36, 0.31, 0.5, 0.56, 0.62, 0.08, 0.67, 0.9, 0.68,
0.51, 0.49, 0.37, 0.68, 0.51), RGE16 = c(NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
1.72, 1.78, 1.14, 0.53, 0.6, 1.45, 1.62, 0.88, 1.94, 1.87,
1.19, 0.77, 0.71, 0.76, 0.69, 0.68, 0.52, 0.64, 0.58, 0.71,
0.64, 0.56, 0.49, 0.48, 0.57, 0.42, 0.86, 0.78, 0.53, 0.54,
0.86, 0.38, 0.59, 0.56, 0.69, 0.49, 0.99, 0.94, 1.22, 0.52,
0.67, 0.43, 0.74, 0.56, 1.01, 0.88, 0.87, 0.7, 0.78, 0.6,
0.86), RGE17 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, 3.64, 2.95, 1.97, 1.98, 3.28, 2.38, 1.72, 1.08, 0.9,
0.13, 1.05, 1.1, 0.87, 1.05, 0.88, 0.72, 0.5, 0.42, 0.34,
0.32, 0.45, 0.54, 0.57, 0.21, 0.2, 0.34, 2.49, 1.41, 0.9,
0.4, 0.28, 0.42, 1.28, 1.56, 1.08, 0.56, 0.3, 0.36, 0.71,
0.96, 0.32, 0.75, 0.58), RGE18 = c(NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, 3.94, 2.52, 2.58, 2.12, 2.04,
1.8, 1.68, 1.87, 1.33, 1.59, 1.31, 1.61, 0.97, 1.46, 1.72,
1.28, 0.42, 1.28, 0.78, 0.83, 0.96, 1.07, 1.2, 1.71, 1.36,
1.64, 1.66, 1.48, 0.34, 1.02, 1.07, 0.84, 1.28, 0.77, 0.5,
1.15, 1.03, 0.95, 0.75, 0.87, 0.56), RGE19 = c(NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, 2.48, 4.02, 3.23, 2.52, 1.67, 2.07, 2.15,
1.83, 1.6, 0.37, 1.39, 1.5, 1.83, 1.22, 1.28, 1.32, 1.3,
1, 0.53, 0.73, 1.33, 1.61, 1.46, 0.56, 0.56, 0.62, 0.2, 0.96,
1.26, 1.44, 1.41, 1.44, 0.91, 1.05, 2.02, 0.46, 0.56, 2.02,
1.82, 1.14, 1.55, 0.64, 0.8, 0.68, 0.82, 0.47, 0.4), RGE20 = c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 3.34,
4.32, 2.89, 2.71, 0.96, 1.12, 1.42, 1.02, 0.97, 0.45, 0.94,
1.05, 0.79, 0.81, 0.91, 0.78, 0.7, 0.5, 0.6, 0.94, 0.98,
0.54, 0.56, 0.7, 0.61, 0.64, 1.08, 1.11, 1.23, 1.14, 0.82,
0.88, 0.46, 0.72, 1, 0.25, 1.1, 0.33, 0.52, 0.54), RGE21 = c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, 2.57, 2.17, 2.32, 2.23, 1.53, 1.82, 1.84, 1.43,
2.32, 3.3, 2.62, 2.85, 3.6, 3.11, 1.1, 1.68, 1.74, 1.59,
1.7, 1.28, 1.08, 1.34, 1.23, 1.26, 1.17, 1.2, 0.91, 0.78,
0.42, 1.06, 1.06, 1, 0.85, 0.78, 0.8, 1.2, 1.3, 1.04, 1.03,
1.14, 0.79, 0.94, 1.11, 0.56, 0.28, 0.45, 0.5, 0.41, 1.11,
0.87, 1.74), RGE22 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, 2.14, 1.75, 1.35, 2.95, 1.51,
1.88, 1.83, 0.78, 0.87, 1.03, 1.05, 1.72, 0.98, 2.18, 2.73,
1.68, 1.82, 2.22, 0.85, 0.71, 1.46, 1.06, 1.3, 1.06, 0.42,
1.31, 1.43, 1.88, 0.97, 0.57, 1.12, 1.07, 1.2, 0.64, 1.44,
0.99, 0.61, 0.71, 0.88, 0.53, 0.51, 0.71, 0.87, 0.53, 1.14,
1.05, 0.85, 0.45, 0.46, 0.37, 0.51, 0.32, 0.65), RGE23 = c(1.18,
1.04, 1.04, 1.18, 1.39, 0.97, 1.32, 1.85, 1.27, 0.99, 0.88,
1.12, 1.2, 0.99, 0.89, 0.88, 0.4, 1.02, 0.8, 1.16, 1.02,
0.75, 0.77, 0.63, 0.72, 0.66, 0.44, 0.62, 0.64, 0.84, 1.16,
0.62, 0.75, 0.64, 0.22, 0.32, 0.82, 0.41, 0.38, 0.7, 0.63,
1.1, 1.01, 1, 0.98, 1.3, 1.02, 0.9, 0.88, 0.46, 0.93, 0.76,
1.16, 0.31, 0.28, 0.72, 0.38, 0.64, 0.71, 1.1, 1.32, 0.76,
0.8, 0.88, 0.93, 0.76, 0.63, 0.7, 0.88, 1.34, 0.77, 0.53,
0.93, 1.04, 1.23, 1, 1.38, 1.08, 0.43, 0.4, 1.16, 1.01, 0.92,
0.57, 0.32, 0.71, 0.52, 0.74, 0.52, 0.69, 0.34, 0.36, 0.4,
0.55, 0.23, 0.28, 0.6, 0.53, 0.64, 0.33), RGE24 = c(NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, 2.88, 1.2, 1.37, 1.35, 2.15, 1.52, 1.42, 1.1, 1.07,
0.77, 0.78, 0.63, 0.88, 0.31, 0.71, 0.52, 0.66, 0.58, 0.79,
0.73, 0.68, 0.34, 0.47, 0.45, 0.36, 0.28, 0.43, 0.63, 0.42,
0.57, 0.53, 0.52, 0.43, 0.51, 0.34, 0.24, 0.24, 0.34, 0.56,
0.25, 0.6, 0.36, 0.45, 0.32, 0.35, 0.36, 0.19, 0.44, 0.39,
0.28, 0.22), RGE25 = c(NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, 2.58, 1.7, 0.54, 0.55,
0.77, 0.54, 0.72, 0.93, 1.42, 0.85, 1.12, 1.12, 1.16, 0.8,
1.41, 0.8, 0.65, 0.97, 0.44, 0.34, 1.03, 0.7, 0.48, 0.88,
0.44, 0.52, 1.39, 1.1, 1.65, 0.43, 0.9, 0.39, 1.17, 1.32,
1.02, 1.15, 1.56, 0.83, 0.98, 1.12, 0.76, 1.9, 2.77, 2.14,
1.24, 2.04, 1.17, 1.98, 0.28, 0.64, 0.96, 1.62), RGE26 = c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, 2.42, 1.39, 2.22, 0.54, 1.76, 1.58, 2.25, 1.2, 1.22,
0.44, 1.27, 1.32, 2.05, 2.01, 2.7, 1.79, 2.25, 3.6, 2.4,
1.77, 1.46, 1.9, 1.95, 1.94, 1.54, 2.2, 2.13, 1.45, 0.92,
0.6, 0.68, 1.41, 1.43, 0.76, 0.23, 0.98, 0.75, 0.48), RGE27 = c(NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,
NA, NA, NA, NA, NA, 2.42, 2.45, 2.27, 2.22, 1.22, 0.7, 2.54,
1.82, 3.04, 3.47, 3.48, 2.75, 2.05, 2.04, 1.86, 2.36, 2.39,
2.32, 1.2)), .Names = c("RGE01", "RGE02", "RGE03", "RGE04",
"RGE05", "RGE06", "RGE07", "RGE08", "RGE09", "RGE10", "RGE11",
"RGE12", "RGE13", "RGE14", "RGE15", "RGE16", "RGE17", "RGE18",
"RGE19", "RGE20", "RGE21", "RGE22", "RGE23", "RGE24", "RGE25",
"RGE26", "RGE27"), row.names = c("1913", "1914", "1915", "1916",
"1917", "1918", "1919", "1920", "1921", "1922", "1923", "1924",
"1925", "1926", "1927", "1928", "1929", "1930", "1931", "1932",
"1933", "1934", "1935", "1936", "1937", "1938", "1939", "1940",
"1941", "1942", "1943", "1944", "1945", "1946", "1947", "1948",
"1949", "1950", "1951", "1952", "1953", "1954", "1955", "1956",
"1957", "1958", "1959", "1960", "1961", "1962", "1963", "1964",
"1965", "1966", "1967", "1968", "1969", "1970", "1971", "1972",
"1973", "1974", "1975", "1976", "1977", "1978", "1979", "1980",
"1981", "1982", "1983", "1984", "1985", "1986", "1987", "1988",
"1989", "1990", "1991", "1992", "1993", "1994", "1995", "1996",
"1997", "1998", "1999", "2000", "2001", "2002", "2003", "2004",
"2005", "2006", "2007", "2008", "2009", "2010", "2011", "2012"
), class = "data.frame")
Thank you for your help!
--
---
Catalin-Constantin ROIBU
Lecturer PhD, Forestry engineer
Forestry Faculty of Suceava
Str. Universitatii no. 13, Suceava, 720229, Romania
office phone +4 0230 52 29 78, ext. 531
mobile phone +4 0745 53 18 01
+4 0766 71 76 58
FAX: +4 0230 52 16 64
silvic.usv.ro
[[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.
