similar to: save grid

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

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

Hi   I want to make a grid map using diversity indices. I have a soil grid map with 7 class in r. how can I make a diversity map based on soil map in r?  I found something for calculation diversity in a table (.csv).   Thanks [[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   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, Doing more search i've discovered package RGDAL that can write a geotiff file with projection. I saved a geotiff file in UTM projection and if i read the file back in R and check the projection seems that everything is OK. But if i load the file in ArcGIS (ESRI product) i get the warning that the file is missing spatial reference so it cannot be projected, but it is displayed correctly.

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

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.

Hello all   I have a geology map that has three level, bellow   <-geology lithology    landscape   landform     landform level is used as covariate (with codes=1,2,3,4,5) for training of neural network, but this level has missing data as NA. I want to replace the missing data of landform level with 0 (zero). Finally, landform will have codes  0,1,2,3,4,5.   please help me   Thanks alot.

 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

Hi, I have a matrix of data which i can vizualize as an image - for example. I would like to save this image as a geotiff file or at a tiff file with a world file which holds the projection of my data (ultimately the data represent a map of some sort). I know i can save the data as an ESRI grid, but i am not interested in that. I wonder if anybody knows about any code which will help me do

R-geographers... I'm trying to solve a problem to implement a line-by-line tiled processing using RGDAL (read 1 line of an image, process the one line, write one line of the image to a binary file). Everything except for the final step I'm able to do using a combination of RGDAL and r-base commands. Below is the basic structure, the input file "elev" can be any image file

Dear R users, I have a list of objects of type "im" > mylist$sp1 $sp2 $sp3 and I want to convert them to a list of objects of class "SpatialGridDataFrame" This works for a single object of class "im": a <- mylist$sp1 b <- as(a, "SpatialGridDataFrame") Then I want to write each element in the new list as a TIFF file, with the name of the elements

Hi, I am very new to R and I am trying to write KML files. I am running R-2.13.1 and I have the following packages: - AKIMA - RGDAL - RSAGA - SP I was trying to run the simple example that I found at the following web page: http://spatial-analyst.net/wiki/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid)

Dear All, The following code save my graphs as pdf: pdf("j:/mix.pdf", width = 18, height = 16) grid.newpage() pushViewport(viewport(layout = grid.layout(3,1))) vplayout <- function(x, y) viewport(layout.pos.row = x, layout.pos.col = y) print(Aplot, vp = vplayout(1, 1)) print(Bplot, vp = vplayout(2, 1)) print(Cplot, vp = vplayout(3, 1)) dev.off() How can I save it in

Dear R helpers, I'm new to R and struggling with applying functions to a spatial grid data frame, called data for example.data has a variable named "slope" and "uparea" Basically, I need to calculate the tan(data$slope)Next, I need to divide data$uparea by data$slope and take the log() of that all In the basic examples I found how to make simple multiplications with a

Dear R community I have stored the results of arrays in a list consist of J-components (say 200 components). Each component containing same no of columns but may be different no of rows. e.g [[1]] [,1] [,2] [,3] [,4] [,5] [1,] 4 0 0 0 0 [2,] 4 3 4 0 0 [3,] 4 3 4 0 0 [4,] 4 3 0 0 0 [[2]] [,1] [,2] [,3] [,4] [,5]

Dear friends I have an array consist of r-rows and c-columns e.g. x=c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,4,0,0,0,0,0,0,0,0); x1=array(x, dim=c(4,6)) output is > x1 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1 2 3 4 0 0 [2,] 1 2 3 4 0 0 [3,] 1 2 3 4 0 0 [4,] 1 2 3 4 0 0 How can i ignore columns having zero sums? Help in this regard

¡Buenas tardes! Antes de nada, dar las gracias por permitirme participar en esta lista. Soy un usuario recién llegado a R de la mano del análisis geoespacial. Soy un arqueólogo que se ha topado con un package para R creado por un arqueólogo para la ejecución de procesos para a creación de rutas de menor coste cuyo nombre es leastcostpath y que se encuentra en el repositorio de GitHub. Contiene