similar to: (no subject)

I generated maps with the function symbols (graphics). These are basic maps generated with : symbols(x,y,circles=myvariable) where x et y are spatial coordinates corresponding to replicates of "myvariable". I would associate legend to this kind of maps, is it possible? Regards, Marion. -- Marion Jacquot Laboratoire de Chrono-environnement UMR UFC/CNRS 6249 USC INRA Universit? de

We want to compute a pseudo R-squared for a model whose parameter estimation was based on maximum likelihood (function likfit, package geoR). I tried to compute the R2 proposed by Maddala (1983) which compare the maximized likelihood for the model without any predictor and the maximized likelihood for the model with all predictors. I got a really low value (0,01%). Did I miss something? Are

Hi I want to sum over one of the dimensions of a n x k1 x k2 array in which each column is the product of the corresponding columns from two matrices with dimensions n x k1 and n x k2. I can see two approaches: a loop on k1 and a loop on k2. But I cannot figure a solution that avoids the loop? Is it possible? (I don't refer to apply or lapply etc either as they are just hidden loops so,

Dear Roger, I am trying to use the write.polylistShape() function of maptools for the first time and realize that it handles list of polygons of class 'polylist'. However, it seems that no as.polylist() function exist in the package. The question behind that is: in your opinion, which would be the best way to convert a list of matrix of polygon nodes coordinates into an object of

Dear Lister, I am working with Windows XP and R 2.1.0 and can check and build home-made packages easily (just *.r, *.rmd, *.rda files, no compiled code). However for some reasons, I cannot check or build the package 'foreign' from the source (I took it as an exercise...). After some lines of sentences OK (here omitted and replaced by ...), I get a message like this: RCMD check

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

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

Exception type exceptions.ValueError Exception Handler Information Traceback (most recent call last): File "/usr/share/rhn/rhn_applet/rhn_applet.py", line 586, in refresh self.model.refresh(force) File "/usr/share/rhn/rhn_applet/rhn_applet_model.py", line 447, in refresh source.refresh(force) File "/usr/share/rhn/rhn_applet/rhn_applet_source.py", line

I have recently configured two servers to be controlers of my domain. The first one is the PDC and is configured like : domain master = yes prefered master = yes local master = yes domain logons = yes security = user os level = 99 And the second, the BDC like : domain master = no prefered master = no local master = no domain logons = yes security = user

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

I am working with the splm package. I use the spgm function: general estimation of a panel data model. Based on this approach, I know it is possible to compute a R2, eg the ratio of variation explained by a given model. My model is : bivmod<-spgm(logIKA~NBLITRE0+NBLITRE1,data=mydatap,listw=comsKnn.nbW,spatial.error=TRUE) I know that we can calculate the R^2 as the variance of the fitted

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

Dear R community, I have currently some problems with non linear regression analysis in R. My data correspond to the degradation kinetic of a pollutant in two different soil A and B, x data are time in day and y data are pollutant concentration in soil. In a first time, I want to fit the data for the soil A by using the Ct = C0*exp(-k*Tpst) with Ct the concentration of pollutant at time t, C0

Dear R community, Please, can you help me with a problem concerning bootstrap. The data table called «RMika», contained times (Tps) and corresponding concentration of a chemical in a soil (SolA). I would like to get, by bootstraping, 10 estimations of the parameters C0 and k from the function: SolA = C0*exp(-k*Tps). # First, I fit the data and all is OK >

Hello, I am using the function mda of the mda library in order to discriminate 4 groups with 8 explanatory variables. I only have 66 observations. I tested all possible combinations of those variable and run for each the Mixture Discriminant Analysis. For some iterations, I got an error message: "error in kmeans(xx, start): initial centers are not distinct". I understood that the

Hello togehter, i have a litte problem to convert a data.frame. My data.frame looks like this one A 1 R5000 2 R4800 3 R4700 4 3500 5 3800 I need now a command, which outputs all the numbers, without the character in front. The solution look like this one: A 1 5000 2 4800 3 4700 4 3500 5 3800 Thanks. Best regards. Mat -- View this message in context:

Hello, It seems that, unlike Splus, the anova.lm function of R does not take the argument ssType=3 into account. How can one get an anova table with the adjusted sum of squares in R? Thanks in advance, 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

Hello, I am trying to use "bwplot" to display whisker boxes according to some conditioning factors ("age" has two values 1/2). I get the following messages: > library(trellis) > bwplot(dvk95~age|site*season,panel=function(x,y){panel.bwplot(x,y)}) Error in tapply(1:0, structure(list(INDICES = numeric(0)), .Names = "INDICES"), : arguments must have

Dear all, I am trying to model variation of distribution of species assemblages according to environmental variables. For that I use a log linear multinomial regression. In order to select variables that mostly discriminate the assemblages, I tried to apply a hierarchical partitioning protocol to my data set. For that I have adapted the all.regs() for multinomial model. The problem is that I

The main limitation of the shapefiles package that I put together is that it does not create shapefiles from R objects - rather it only writes shapefiles that have been read into R and manipulated within the constraints of the existing file structure. By this I mean that for example you can change the coordinates of points and write them back out. Or you can add a bunch of blank columns in the