This example illustrates your problem:
> plot(1:3,1:3)
> plot(1:3,1:4)
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' and 'y' lengths differ
Now all you have to do is figure out which of your plot() commands gives
the error, and check the values you are giving to it to find out which
have different lengths.
Please also read the posting guide.
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
On 5/15/13 11:06 AM, "catalin roibu" <catalinroibu at gmail.com>
wrote:
>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 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.