similar to: replace NA

Dear All   I want to import the vector file (   .shp) to R. I could import the file by rgdal package before, by following:   geology<-readOGR('C:/geology//saga/geo.geom','finalgeology')   but now there is an error:   Error in ogrInfo(dsn = dsn, layer = layer, input_field_name_encoding = input_field_name_encoding) :           GDAL Error 4: .shx file is unreadable, or corrupt.

I applied bootpred for multinomial logistic reg. (with nnet package). I used same as theta.fit and theta.predict of R for my data. but give me error. Can I do this with response vriable;7 levels predictor variables:5 (1 classifier, 4 continuous)?   Thanks alot Azam   [[alternative HTML version deleted]]

Hi   I want to do leave-one-out cross-validation for multinomial logistic regression in R. I did multinomial logistic reg. by package nnet in R. How I do validation? by which function? response variable has 7 levels   please help me   Thanks alot Azam [[alternative HTML version deleted]]

Dear All   Pleas help me to increase the memory in R.   In order increase memory, I read the FAQ and follow the instruction as below     Close R, then right-click on your R program icon (the icon on your desktop or in your programs directory). Select ``Properties'', and then select the ``Shortcut'' tab. Look for the ``Target'' field and after the closing quotes around the

Hi all, This is more a help on ideas that on actuall R code. Those of you which are geologist or work with boreholes woudl understand that i would like to achieve. I'll tr to explain as good as I can...so here it goes.. I have around 1000 geomechanical (geological) borehole logs in a database. The logs contain staial info, lithology and some laboratory data. My initial goal is

Hi   I did validation for prediction by logistic regression according to following:   validationsize <- 23 set.seed(1) random<-runif(123) order(random) nrprofilesinsample<-sort(order(random)[1:100]) profilesample <- data[nrprofilesinsample,] profilevalidation <- data[-nrprofilesinsample,] salich<-profilesample$SALIC.H.1 salic.lr<-glm(salich~wetnessindex, profilesample,

Hi   I want to draw a plot from observed and predicted data and also shows threshold and data before threshold are identified with different color from data after threshold.   Suppose: abserved data are 0 or 1 predicted data= 0 to 1 threshold=0.5   Thanks alot     [[alternative HTML version deleted]]

hello, I am very new to R. My current data set is a mix of values and categories. It is a geoscience data set, with values per rock sample. Case in point, each sample belongs to a lithology class, and each sample has several physical property measurements (density, porosity...). I want to be able to plot these physical properties for all samples in each lithology class. this is how i'm doing

Red Hat Enterprise Linux Server release 6.5 (Santiago) selinux is disabled. The following commands were all run on the RedHat Server on which I am running samba. *The following ports are open* 5 ACCEPT tcp -- 0.0.0.0/0 0.0.0.0/0 state NEW tcp dpt:137 6 ACCEPT tcp -- 0.0.0.0/0 0.0.0.0/0 state NEW tcp dpt:138 7 ACCEPT tcp --

Hello, I am trying to run an ordinal logistic regression (polr) using the package 'MASS'. I have successfully run other regression classes (glm, multinom) without much problem, but with the 'polr' class I get the following error: " Error in svd(X) : infinite or missing values in 'x' " which appears when I run the "summary" command. The data file is

Hi, I am interested in ggplot2 and I found this lattice code very interesting (http://addictedtor.free.fr/graphiques/graphcode.php?graph=48). Code: library(lattice) lattice.options(default.theme = canonical.theme(color = FALSE)) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species = c("ArisDiff",

Dear R users, If I manage to sort out this graph, it is certainly a candidate for the new R graph gallery (http://addictedtor.free.fr/graphiques/displayGallery.php)! I created the following lattice graph: library(lattice) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species =

I am trying to plot a non-linear trend line along with my data. I am sure the solution is simple; any help greatly appreciated. Recent and historic attempts: fit = nls((sd~a+b*exp(-c/time)+d*geology), start=list(a=1, b=1, c=10, d=-1), model=TRUE) plot(time, sd, col=3, main = "Regression", sub=formula(fit)) lines(time, fitted.values(fit), lwd=2) #tt = seq(from=0, to=200, by=10)

Can someone point me to an implementation of the "runs up and runs down" test, or does such beast exist in R? From web searches the "runs up and runs down" test is commonly used for testing pseudo-random number generators and in simulations. John C. Davis describes its use in geology in his book on geological statistics (2002). Having searched the full site, I can locate only

Hi, I have a 3-d geology problem. I have an xyz fault data set. I am able to view the data set in R using plot3-d but cannot create a raster from it. I need to create a raster and export so that it can be read in arcscene. Can anybody help? Thanks -- Shane [[alternative HTML version deleted]]

Dear R graphics gurus, Another question about lattice graphics. This time I would like to plot means and confidence intervals by group factor in a lattice graph. I can not find any working lattice examples. Maybe a custom panel function is the answer, but that is a bit beyond me for now. The individual plots within the lattice graph could look like this: # Example with confidence intervals

Hi, I would like to put the sample number beside each lable in a boxplot. How do I do this? Essentially, I need to count the sample size for each factor, see below: Thanks boxplot(DATA$K_Merge~factor(DATA$UnitName_1),axes=FALSE,col=colours) title(main=list("Tukey Boxplot by Geology:\n K(%)",cex=cexlb)) axis(1, 1:21, labels=FALSE, las=2) text(seq(1, 21, by=1), par("usr")[3],

Dear R users, I face to a nested pattern and despite the numerous examples in the help I am still confused. I sampled bugs in different habitats within sites which were within rivers themselves within different regions. The habitat correspond to different substrata (not systematically present in all sites). For rivers and sites, I have environemental variables (e.g. altitude and slope of

Hi list   I want to define a moving window based on ferquency in a grid file. is there a package for this? I found packages for moving window mean, median or quantile.   thanks  [[alternative HTML version deleted]]

 Hello   I have questions about gbm package.  It seems we have to devide data to two part (training set and test set) for first.   1- trainig set for running of gbm function 2- test set for gbm.perf      is it rigth? I have 123 sample that I devided 100 for trainig and 23 for test.   So, parameter of cv.folds in gbm function is for what?   Thanks alot Azam       [[alternative HTML