similar to: type III sum of squares - ssType

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

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

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

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

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

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'd like to propose an extension to the function summary.aov. In Splus (2000, I don't know about other versions), summary.aov allows a parameter ssType to be set to 1 or 3 (defaults to 1) to choose the type of Sums of Squares. I know I can get Type III SS in R with drop1(model), but including the functionality into summary.aov would, in my opinion, - yield a more usable table

I have two questions I would appreciate assistance with: (1) I believe "by" is the command used to split a file. In the following example, "mydat" is the dataframe , "group" the variable I want to split the analysis by & "WK1FREQ,WK2FREQ" the variables attach(mydat) by (GROUP) cor.test (WK1FREQ,WK2FREQ) I have also tried: (by (GROUP) cor.test

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

I'm an SPSS user trying to make the transition to R. Can someone help me translate the following SPSS code into R?: GLM Total_tp1 Total_tp2 WITH Age Sex   /WSFACTOR=Time 2 Repeated   /METHOD=SSTYPE(3)   /CRITERIA=ALPHA(.05)   /WSDESIGN= Time   /DESIGN= Age Sex Age*Sex. Also. can anyone recommend any resources to help SPSS users learn to things in R? Thanks, -kristi [[alternative HTML

Dear R people: I wish to switch from SPSS to R, but there is one particular type of ANOVA design that cannot be done in R. Or more likely, it can be done, but it is nowhere documented. The problem is typical for psychologists: You have a repeated measures design with different groups of subjects. Now, this can be done with the aov command, but the number of subjects in both groups must be

Hello, I use lme4 package. library(lme4) fit=lmer(y ~ time+dye+trt+trt:time + (1|rep), data=dataset, na.action='na.omit') anova(fit) The anova gives sequential F-tests and sequential SS. My question is: how I can get partial F-tests and partial SS? For lm (not lmer) anova(lm(y~x+z)) we can use anova(fit, ssType=3) but it is not work for lmer. Natalia.

Hello everybody, This is my firth question in this mailing list. The question is about the how to deal with a model with both random and fixed factors. I know to do the analysis using SPSS. The command line wit be as follows: UNIANOVA alpha.ug.mg BY Line Year Rep /RANDOM=Year Rep /METHOD=SSTYPE(4) /INTERCEPT=EXCLUDE /PRINT=HOMOGENEITY /CRITERIA=ALPHA(0.05) /DESIGN= Year Rep*Year Line

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

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