similar to: DEV2bitmap: jpeg with res=400 not enough for CORELDRAW poster A0

Hi R helpers, Is there a function (package?) in R available which tests "spatial autocorrelation" between 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

Dear all, Is there a simple trick to read in data with the following format and create a Time Serie object of it? Date CountOfField2 5/10/1998 7 5/11/1998 5 5/12/1998 2 5/14/1998 1 5/15/1998 1 5/19/1998 1 5/20/1998 1 5/21/1998 1 5/24/1998 2 5/25/1998 1 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

Dear R-list, Is there a function or technique by which the probability (or log odds) range of 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

Hi all, I created a graph by plotting 4 time series on top of each other but don't know how to add extra axes on the right hand side of the plot (or leftside). For the first two time series it axis are plotted but the last two time series don't have and axis. How can these be added and where? e.g. par(mfrow=c(1,1)) plot(timeserie1) par(new=T) plot(timeserie2,

Dear R list, From the lrm() binary logistic model we derived the G2 value or the likelihood-ratio chi-squared statistic given as L.R. model, in the output of the lrm(). How can this value be penalized for non-linearity (we used splines in the lrm function)? lrm.iRVI <- lrm(arson ~ rcs(iRVI,5), penalty=list(simple=10,nonlinear=100,nonlinear.interaction=4)) This didn’t work

Dear R-list, *When comparing two logistic regression models with the anova CHi test, I obtain the following error: (there are no NA's in the time series). How can this be solved such that I can compare two models on the same dataset were different explanatory variables are used? l.KBDI <- glm(zna.arson2 ~ zna.KBDI,family = binomial) l.NDWI <- glm(zna.arson2 ~ zna.NDWI,family

Dear r-list, When I try to plot the following 3D lrm fit I obtain only arrows with labels on the three axes of the figure (without values). fit <- lrm(y ~ rcs(x1,knots)+rcs(x2,knots), tol=1e-14,X=T,Y=T) dd <- datadist(x1,x2);options(datadist='dd'); 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

Dear R-helpers, Is it possible to create a type of 'movie' in 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

Dear, Thanks for the previous tips about 'its' for importing the following data. 5/10/1998,7 5/11/1998,5 5/12/1998,2 5/14/1998,1 5/15/1998,1 5/19/1998,1 5/20/1998,1 1. When using the following command; test<-readcsvIts('Fires98.csv',informat=its.format("%m/%d/%Y"),header=FA LSE) the function reads in the data from the csv file as; V2 05/10/1998 7

Hi, When running z <- plot(fit, age=NA, cholesterol=NA, perim=boundaries, method='image') Legend(z, fun=plogis, at=qlogis(c(.01,.05,.1,.2,.3,.4,.5)), zlab='Probability') And after pointing the cursor to the plot() screen in R, I obtain the following message: Using function "locator(2)" to place opposite corners of image.legend Error in

Hi all, How can the label of the x-axes in the plot() of a stl.object be adapted? 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

Hi R-list, I'm working with irregular time series (time series of climate data, daily data 365/6 days a year) and would like to create regular time series from them ( irts e.g. Rain <- irts(as.POSIXct(Climate[,1]),Climate[,5]) to ts e.g. test <- ts(x, start=c(1997,1), frequency=365) ) such that I can find where the gaps (lacking temperature data, ...) are and try

Hi R-help, We are using R 2.2 on Win XP and have a detail question on how the legend can be optimised. We use the following; -> plot(,type="b",...) The lines in the plot do not cross the points. How can we obtain the same effect in the legend? (points without a line through them..) We tried setting the pt.bg to white but this did not help. See script below. thanks, Jan

I have seen screenshots of CorelDRAW X3 running in Wine, but have no idea, how to get this working. How do I convert normal CD version into "portable"? I'd really love to get rid of VirtualBox running WinXP just for CorelDRAW. I would even pay for a complete, step-by-step instructions that will result in CorelDRAW X3 running properly in Wine. My setup: Ubuntu 9.04 Jaunty 64-bit,

Dear R-help, How could a cross-correlation plot be optimized such that the relationship between seasonal time-series can be studied? We are working with strong seasonal time-series and derived a cross-correlation plot to study the relationship between time-series. The seasonal variation however strongly influences the cross-correlation plot and the plot seems to be ?rather? symmetrical (max

Dear List, I am using the bioconductor package Category to do some gene enrichment analysis, and usually save my KEGGmnplot's using a dev2bitmap command. This has worked just fine, until suddenly earlier today I got this error-message: > dev2bitmap("04610_080721.jpg",type="jpeg", height = 10, width = 10, res = 200) Error in dev2bitmap("04610_CSF080721.jpg",

System: Linux kernel 2.6.15 Ubuntu dapper R version 2.5.1 ESS 5.2.11 on Emacs 21.4.1 ------------------------- Colleagues Having read the posts on producing perfect graphs in R, or using inkscape to edit R graphics output, I have a related question. Lately I am publishing in a journal that is very picky about their graphics formats. It is more efficient for me to get the final picky details

Hello Wine experts. I am trying to get CorelDraw 7.0 to work under wine, and I have successfully installed the program, but when I attempt to run it, I get the following dialogue error message: Key: "HKEY_LOCAL_MACHINE\SOFTWARE\Corel\CORELDRAW\7.0" Value: "ConfigDir" not found in the Registry. Along with this error message, comes the exciting output displayed below, that is

Dear, When fitting the following model knots <- 5 lrm.NDWI <- lrm(m.arson ~ rcs(NDWI,knots) I obtain the following result: Logistic Regression Model lrm(formula = m.arson ~ rcs(NDWI, knots)) Frequencies of Responses 0 1 666 35 Obs Max Deriv Model L.R. d.f. P C Dxy Gamma Tau-a R2 Brier 701 5e-07 34.49

Hi, I don't know if I am the first one to report the problem: in 2.8.0 dev2bitmap gained aa support, but "extra" is null if neither taa nor naa is specified. > dev2bitmap("plot.pdf",type="pdfwrite") Error in paste(shQuote(gsexe), " -dNOPAUSE -dBATCH -q -sDEVICE=", type, : object "extra" not found >