search for: verbesselt

...? e.g., When plotting: plot(stl(nottem, "per")) In the labels of the x-axes is “time”. How can this be changed to e.g., “Time (dekade) “? It does not work with xlab or others anymore… Thanks, Jan _______________________________________________________________________ Ir. Jan Verbesselt Research Associate Biosystems Department ~ M³-BIORES Vital Decosterstraat 102, 3000 Leuven, Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/ _______________________________________________________________________ Disclaimer: http://www.kuleuven.be/cwis/email_discl...

...n points (e.g. vector layer of weather stations)? (e.g. Moran's I...) Via the archives we found out that there is a package 'spdep' which uses grid data for testing spatial autocorrelation. Thanks a lot, Jan _______________________________________________________________________ Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

...n R based on the output of several figures (e.g., jpegs) via the plot() function. I obtained dynamic results with the plotting function and would like to save these as a movie (e.g., avi or other formats)? Regards, Jan _________________________________________________________________ Ir. Jan Verbesselt Research Associate Group of Geomatics Engineering Department Biosystems ~ M³-BIORES Vital Decosterstraat 102, 3000 Leuven, Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/ _______________________________________________________________________ [[alternative HTML v...

...ge.legend Error in Legend.plot.Design(z, fun = plogis, at = qlogis(c(0.01, 0.05, : couldn't find function "subplot" *How can I position the legend inside the range? (or solve the following error message). Regards, Jan -----Original Message----- From: Jan Verbesselt [mailto:Jan.Verbesselt@biw.kuleuven.be] Sent: Monday, August 22, 2005 4:59 PM To: 'r-help@stat.math.ethz.ch' Subject: How to add legend of plot.Design function ( method=image)? Dear Rlist, How can the Legend of the plot.Design() function can be visualized? Following the docu...

...5/26/1998 2 .... 2002 ... R should recognize that some dates are not available...(NA). You can define start and end date Ok, and frequency= 365 is ok...but is it possible that recognizes the gaps? Thanks a lot, Jan ________________________________________________________________________ __ Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://perswww.kuleuven.ac.be/~u0027178/VCard/mycard.php?name=janv http://gloveg.kuleuven.ac.be/ __________________________________________...

...2 110 113.593 0 9.635 Warning message: NaNs produced in: pchisq(q, df, lower.tail, log.p) *Is this caused by the non-linearity of the data? *How can I obtain a p-value? Thanks, Jan _______________________________________________________________________ ir. Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

...tes=as.POSIXct(x=strptime(dimnames(V1)[[1]] format=its.format("%m/%d/%Y"))) )...Doesn't work! I have difficulties defining the first V1 column as its object. Thanks in advance for helping me out, Jan ________________________________________________________________________ __ Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://perswww.kuleuven.ac.be/~u0027178/VCard/mycard.php?name=janv http://gloveg.kuleuven.ac.be/ __________________________________________...

...ature data, ...) are and try out methods to fill the gaps. The main objective is to detect gaps, how long they are and fill them if possible with average, median values. All ideas & methods are welcome! Regards, Jan _______________________________________________________________________ Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

...icture in CORELDRAW (for a poster A0) format the resolution and colors are not optimal. How can I save pictures (colors/resolution) optimally for import into CorelDraw for an A0 poster? Pdf?/tiff?/bmp?/ Thanks, Jan _______________________________________________________________________ ir. Jan Verbesselt Research Associate Lab of Geomatics Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

...f a logistic model (fit <- lrm(Y~X)) can be derived? The aim is to obtain min & max of the estimated probabilities of Y. Could summary.Design() be used for that or is there another method/trick? Thanks, Jan _______________________________________________________________________ ir. Jan Verbesselt Research Associate Lab of Geomatics Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

...;) par(new=T) plot(timeserie3, type="h", col=2,yaxt="n", ylab="") par(new=T) plot(timeserie4, type="b", col=4,yaxt="n", ylab="") Thanks, Jan _______________________________________________________________________ ir. Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

...rk properly. The aim is to obtain a value that can be used to compare the goodness of fit of the different univariate binary logistic models. (The lower the value, the better the fit) Kind regards, Jan ____________________________________________________________________ Ir. Jan Verbesselt Research Associate Group of Geomatics Engineering Department Biosystems ~ M³-BIORES Vital Decosterstraat 102, 3000 Leuven, Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/ _______________________________________________________________________ [[alternative HTML v...

...ctivated correctly? > acf(ts.Lin, type=c("correlation"), na.action=na.omit) Error in na.omit.ts(as.ts(x)) : time series contains internal NAs ((ts.Lin contains two time series, where one contains internal NAs (-->an NA not a the end/beginning of a time serie))) Thanks a lot! Jan Verbesselt

...9;); par(mfrow=c(1,1)) plot(fit,x1=NA, x2=NA, theta=50,phi=25) How can I add values to the axes of this plot? (axes with the range of values of each of the explanatory variables x1&x2) Thanks, Jan _______________________________________________________________________ Ir. Jan Verbesselt Research Associate Group of Geomatics Engineering Department Biosystems ~ M³-BIORES Vital Decosterstraat 102, 3000 Leuven, Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/ _______________________________________________________________________ [[alternative HTML v...

...("topright", legend.txt,col=1,lty=1, pch=c(19,1),bty="n", pt.bg=1) par(new=T) plot(ts.Y, type="b", lty=1, col=1, pch=1, ylab="", xlab="",yaxt="n",ylim=c(0.1,0.7)) axis(4) mtext(side=4, line=3, "Y", cex=1) -- Ir. Jan Verbesselt Research Associate Lab of Geomatics, K.U.Leuven Vital Decosterstraat 102, 3000 Leuven, Belgium Tel:+32-16-329750 Fax:+32-16-329760 http://gloveg.kuleuven.ac.be/ Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm

...the noise (sensor and atmospheric) with a filter (kalman/ holt-Winters/TsSmooth? Or FFT.) or the BISE method in R? Is the BISE (Best index slope extraction) function already programmed in R I couldn't find it? Much appreciated, Jan ******************************************************* Jan Verbesselt Research Associate Lab of Geomatics and Forest Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel:+32-16-329750 Fax: +32-16-329760 http://perswww.kuleuven.ac.be/~u0027178/VCard/mycard.php?name=janv

Dear all, When fitting an "ols.model", the confidence interval at 95% doesn't cover the plotted data points because it is very narrow. Does this mean that the model is 'overfitted' or is there a specific amount of serial correlation in the residuals? Which R functions can be used to evaluate (diagnostics) major model assumptions (residuals, independence, variance) when

...ical shape? Is ts1 lagged vs. ts2? Jan (*Version R 2.2) Ps. -We tried also ccf() after differencing and decompositioning but seasonality remains in the residuals. -max cross-correlation mostly occurs at lag 0. _______________________________________________________________________ Ir. Jan Verbesselt Research Associate Biosystems Department ~ M??-BIORES Vital Decosterstraat 102, 3000 Leuven, Belgium Tel: +32-16-329750???? Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/ _______________________________________________________________________ Disclaimer: http://www.kuleuven.be/cwis/email_discl...

...Residual Deviance: 31.76 AIC: -167.7 A negative AIC occurs! How can the negative AIC from different models be compared with each other? Is this result logical? Is the lowest AIC still correct? Thanks, Jan _______________________________________________________________________ ir. Jan Verbesselt Research Associate Lab of Geomatics Engineering K.U. Leuven Vital Decosterstraat 102. B-3000 Leuven Belgium Tel: +32-16-329750 Fax: +32-16-329760 http://gloveg.kuleuven.ac.be/

Hi R specialists, When testing the association between two time series the cor.test gives the following message...-> p-values may be incorrect due to tie What does it mean? (it is not described in the help) Thankx, Jan > cor.test(Origi[,1],Origi[,2], alternative = c("two.sided"),method = c("spearman"), conf.level = 0.95) Spearman's rank correlation rho