similar to: conversion of a data.frame of numerics to a data.frame of factors

I am not sure a previous e-mail reached the list (no mail aknowledgement from R-boundle etc.). The question was how to write polygon or segment coordinates into a shapefile set or any other ArcGIS supported format. The library shapefiles seems to do something but the documentation is a bit beyond of my mind.... and I cannot get the meaning of the functions write**** and its application to the case

Dear All, I would like to make partial sums (or means or any other function) of the values in intervals along a sequence (spatial transect) where groups are defined. For instance: habitats<-rep(c("meadow","forest","meadow","pasture"),c(10,5,12,6)) observations<-rpois(length(habitats),2)

Dear all, I'm running into problem in R-1.9.0 that hasn't happened with R-1.8.x If I make a plot with xyplot(), and use the menu to either save to a metafile or copy to clipboard as a metafile to export to eg Powerpoint, I can just copy a blank then. This does not occur with the classical plot(). Thus, I suppose that it may come from lattice or grid (going through the R-archives I

Dear listers, I have written a function to facilitate the drawing of altitude profiles with x (distance), y (altitude) and a z parameter (altitude magnification). profplot<-function(x,y,z=10,...){ op <- par()$mai par(mai=c(0.95625,0.76875,0.76875,0.95625)) plot(x,y*z, type="l",asp=1,las=1,xlab="",ylab="",yaxt="n",...)

Hi, I cannot go through the archives with which() as key-word... so common. Though I am sure to have seen something about this subject in the past could somebody put me on the track. I have a matrix (actually a data.frame) in which I would replace the non-null values by 1. I tried the following: indices<-which(myforetbin > 0,arr.ind=T) myforetbin[indices[,1],indices[,2]]<-1 and get

Dear All, We have a really serious problem. This weekend we upgraded our mixed 1.6.0.1 + 1.6.4.1 cluster based on CentOS 4.4 i386 to CentOS 5.1 x86_64. We reformatted all drbd storages, so recreated the whole cluster and copy back from backup. Windows clients reached the cluster through a samba server, but now it''s just not working at all. This is my smb.conf # testparm Load smb

Is there a fonction for multiple comparison tests (similar to "multicomp" in Splus) in a package of R? Thanks in advance for any hint... Cheers, Patrick Giraudoux University of Franche-Comté Department of Environmental Biology EA3184 af. INRA F-25030 Besançon Cedex tel.: +33 381 665 745 fax.: +33 381 665 797 http://lbe.univ-fcomte.fr [[alternative HTML version deleted]]

Hi, I have just upgraded from the 1.8.1 to the 1.9.0 version of R, and have some trouble to run RWinEdt from the .Rprofile file (in the user folder). The script is: library(MASS) library(lattice) cat("Load editor?(y/n default = y): ") nf <- as.character(readLines(n = 1)) if ((nf=="y")|(nf=="Y")|(nf=="")) {library(RWinEdt)} rm(nf) When run at start,

Hi, Would a specialist of the "point" or "plot" functions try the following: mat<-matrix(c(1.836767,4.025989,0.6396777,0.3764444),ncol=2) plot(mat) points(mat[1,],col="red") ..A lag appears on x for mat [1,1] between the two displays. I wonder if this example may be due to a bug or to the mis-use of a matrix in the plot() points() functions. In case of

Folks, I have got a strange behaviour when testing this: sum(x) != 1 let us set x<-c(70,134,1,5,0) and transform it in a vector of probabilities x<-x/sum(x) One expect sum(x) should be equal to 1, which is apparently the case > sum(x) [1] 1 However, when I try to test it I get: > if(sum(x) !=1) print("lost") else ("OK") [1] "lost" Which means

Dear Aaron, I am really a tool user and not a tool maker (actually an ecologist doing some biostatistics)... so, I take the liberty of sending a copy of this e-mail to the r-help list where capable computer persons and true statisticians may provide more relevant information and also to Paulo Ribeiro and Peter Diggle, the authors of geoR.. I really feel that your huge matrix cannot be

Hi all, I am trying to install add-on in R (rw1070) to work with WinEdit. Libraries Swinregistry and Rwinedt have been installed via "Install Package(s) from local zip files" from R. However, when I run library(Rwinedt) I get the following messages: > library(Rwinedt) Loading required package: SWinRegistry Error in loadNamespace(name) : package `methods' does not have a name

Hello Is anyone else having problems receiving email from the list? No email from R help today! Thanks Marwan ------------------------------------------------------------------- Marwan Khawaja http://departments.aub.edu.lb/~mk36

"in non linear modelling finding appropriate starting values is something like an art"... (maybe from somewhere in Crawley , 2007) Here a colleague and I just want to compare different response models to a null model. This has worked OK for almost all the other data sets except that one (dumped below). Whatever our trials and algorithms, even subsetting data (to check if some singular

Dear listers, I am making a trial to move from S3 to S4... I have created some classes of interest and they work acceptably well for the purpose. I am now wondering how to make them operate in a package. In clear when a package is loaded (eg library(mypackage)) where should I put the class descriptions:

Dear listers, I have developped a set of functions that I would like to package on a Windows XP plateform for some friends (this would be more simple than to deliver them as a source text file without handy help). I am working under Windows XP. Of course I have gone through the manual "Writing R extension" and try to sort out what a most simple "packaging" for beginner

Hi, I am moving from Windows XP to Ubuntu 6.10 and installed R 2.4.0. When I run eg plot.lm (things work fine with plot.default - eg plot(rnorm(30),rnorm(30))) plot(lmobject) I can get the first plot and then this message: Hit <Return> to see next plot: Error in text.default(x, y, labels.id[ind],cex=cex, xpd=TRUE, : could not find any X11 fonts Check that the Font Path is correct

Dear all, After going through the Lattice doc and R-help list and google, I got the feeling that there is no function in lattice or other package to compute a pie chart object of class "trellis". Although pie charts are obviously not considered optimal even in the pie() doc ;-) , pie chart trellis objects would be easy positioned e.g. over a map drawn with the grids package. Can

Dear listers, Here we have a strange result we can hardly cope with. We want to compare a null mixed model with a mixed model with one independent variable. > lmmedt1<-lme(mediane~1, random=~1|site, na.action=na.omit, data=bdd2) > lmmedt9<-lme(mediane~log(0.0001+transat), random=~1|site, na.action=na.omit, data=bdd2) Using the Akaike Criterion and selMod of the package pgirmess

Hi, I am trying to bootstrap the difference between each parameters among two non linear regression (distributed loss model) as following: # data.frame > Raies[1:10,] Tps SolA Solb 1 0 32.97 35.92 2 0 32.01 31.35 3 1 21.73 22.03 4 1 23.73 18.53 5 2 19.68 18.28 6 2 18.56 16.79 7 3 18.79 15.61 8 3 17.60 13.43 9 4 14.83 12.76 10 4 17.33 14.91 etc... # non