similar to: R Commander question

Good morning Does R (or R users) have a formal training workshop/facility? Or do open-source softwares generally put the onus on the user to learn programming techniques? Are the workshops generally offered for the masses or is there one-on-one training available? Does one "go to" R for programming training ... or do R workshops travel around the country? Thanks for your

Good afternoon Has anyone tried using Dr. Elith's BRT script? I cannot seem to run gbm.step from the installed gbm package. Is it something external to gbm? When I run the script itself <- gbm.step(data=model.data, gbm.x = colx:coly, gbm.y = colz, family = "bernoulli", tree.complexity = 5, learning.rate = 0.01, bag.fraction = 0.5) ... I

Dear all, I am currently learning to run logistic regression models with R. Would someone tell me what this exactly means: Estimated scale (compare to 1 ) 1.746724 If the value is higher or lower than 1, what should I do? The complete results of the model were as follows: Generalized linear mixed model fit using Laplace Formula: success ~ diffca + (1 | partner1) + (1 | partner2) + (1 |

Sorry. Let me try again then. I am trying to find "significant" predictors" from a list of about 44 independent variables. So I started with all 44 variables and ran drop1(sep22lr, test="Chisq")... and then dropped the highest p value from the run. Then I reran the drop1. Model: MIN_Mstocked ~ ORG_CODE + BECLBL08 + PEM_SScat + SOIL_MST_1 + SOIL_NUTR + cE + cN +

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 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

Dear all, When saving a plot with the dev2bitmap command: name <- c("test.jpeg") dev2bitmap(name,type="jpeg",height=8,width=13,res=400) Everything seems to be ok... After importing this picture 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

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,

hello r development team, i'm building R 2.3.0 on solaris and when i run the 'make install' i'm getting a syntax error during the "installing etc ..." which causes the installation to fail. i get this error whether i use gnu-make of sun-make, see the error and reasons below. gmake[1]: Entering directory `/export/medusa/darin/build/R-2.3.0/etc' installing etc ...

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 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, 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

Hello All, I'm deploying NUT in a data center and I'm curious to know how others have gone about staging the shutdown of various systems. The systems are broken down into groups of importance, group1 being the most important group4 being the least. In the event of an outage take down group4 after 5 minutes, group3 after 10 minutes, group2 after 15, and finally group4 after 20. I've

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

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

Does anyone know a reason why, in principle, a call to randomForest cannot accept a data frame with missing predictor values? If each individual tree is built using CART, then it seems like this should be possible. (I understand that one may impute missing values using rfImpute or some other method, but I would like to avoid doing that.) If this functionality were available, then when the trees

Good afternoon, This is my first post, so please be gentle ;) Anyway, here's my question....I've recently put together a Linux Server for our company, replacing an old Novel server. I have the box set up with Samba, as a PDC on our network. So far, just about everything is working correctly. One issue I'm coming up against, though, is printing from DOS. I have Samba

I have a set of timestamp data that I have in a text file that I would like to import into R for analysis. The timestamps are formated as follows: DT_1,DT_2 [2006/08/10 21:12:14 ],[2006/08/10 21:54:00 ] [2006/08/10 20:42:00 ],[2006/08/10 22:48:00 ] [2006/08/10 20:58:00 ],[2006/08/10 21:39:00 ] [2006/08/04 12:15:24 ],[2006/08/04 12:20:00 ] [2006/08/04 12:02:00 ],[2006/08/04 14:20:00 ] I can get

Dear R-users, Is there a function in R that classifies data in more than 2 groups using logistic regression/classification? I want to compare the c-indices of earlier research (lrm, binary response variables) with new c-indices obtained from 'multiple' (more response variables) logistic regression. Best regards, Stephanie Delalieux Department Biosystems M?-BIORES Group of Geomatics

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