similar to: Point labels

Dear R What does NaN mean? I recently did a correlation on a batch of data for some reason it didn't like one column cor(sleep,use="complete.obs") BodyWt BrainWt SlowSleep ParaSleep TotalSleep BodyWt 1.00000000 0.95584875 -0.3936373 -0.07488845 -0.3428373 BrainWt 0.95584875 1.00000000 -0.3867947 -0.07427740 -0.3370815 SlowSleep -0.39363729

I am using R for a class which uses S. Often I can not find functions that in S but not in R, example, multicomp(). What is your advise. Thanks. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body",

[Moderator: This was erronously sent to R-announce, and filtered fortunately ] Hi, This is my first time to use R. I'm wondering which function I can use to do multiple comparison. I have an lm object and want to do multiple comparison based on this object. Splus has a function multicomp() but I could not find a similar one in R. Thanks in advance. Really appreciate it. Julie

Hello, I would like to know if this behaviour in R is as expected. I have a data frame 'dat' with column var1 being character (and not factor). Then I create a new column 'var2' by: > dat[,"var2"] <- dat$var1; Column var2 is now a factor. But if I do: > dat$var2 <- dat$var1; Then column var2 is character (and not factor). I don't want to have

I suppose that multcomp in R and multicomp in S-Plus are related and it appears that it is possible to use multicomp with lme in S-Plus given the following correspondence on s-news sally.rodriguez at philips.com 12:57 p.m. 24/04/03 -0400 7 [S] LME summary and multicomp.default() Is it possible to use multicomp with lme in R and if so what is the syntax from a simple readily available

Dear R people, I have heard rumors of a function called (I believe) multicomp in S+ for doing multiple comparisons for linear models. It doesn't seem to be in the version of R I am using, 0.90.0. Is there something equivalent to this in R, or does it simply not exist? I see there is are more recent versions of R. For my platform, SuSE Linux, the most recent rpm available appears to be

Dear list members, First of all thank you for your helpful advices. After your answeres to my firt mail I studied a lot (R-News n?5) and I tried to perform my analysis: First, to fit a GLM with a nested design I decided to use the function "lmer" in package "lme4" as suggested by Spencer Graves and Filippo Piro. I remember to you that my data were: land use classes, 3 levels

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 have a 16 bit image (TIFF) and i want to analyse the pixels distribution. So, i obtain a matrix which values are between 0 and 2^16 -1. Now i would like to represnt this image with the fucntions rgb() and image(). I am not sure , but i think that only 256 colors are available. So is there a solution to represent all the palette of the colors or i have to limit the representations with

Hi all, I was wondering whether it is possible to change the angle of all axis in a scatterplot3d graph. I know that there is one parameter <angle> but I was wondering whether it is possible to control the angle of all different axis individually (eg by using <anglex><angley><anglez>). In addition I was wondering whether there is a way to change the angle of a plot

Hello, I have a data frame whose columns are factors with numeric levels and I want to convert the columns to numeric so that I can treat the data frame as a matrix. I found a way of doing this but would like to know if there is an easier way. My way is: M <- matrix(ncol=ncol(df) , nrow=nrow(df) , as.numeric(as.matrix(df))) where 'df' is the data frame. I am wondering if there is

Hello everyone! If this message was posted twice, my apologies. I'm having issues with my webmail today... I am using the testing version of NUT from SVN at changeset 582. I have successfully created a stub for newhidups driver for my Dynex DX-800U UPS, but I am completely lost on the customization of it. I've followed the hid-subdrivers.txt instructions and I understand I need to

Dear R-users, again two question... # Question 1 Adding two lines of text to a plot, I am using: # ------------------------------- plot(k[,1], k[,2], pch=16, ylim=range((min(k[,2])-0.2):(max(k[,2])+1))) a <- paste("Cor.:" ,cor(k[,2],k[,1])) b <- paste(nrow(k), "Countries") text(90, max((k[,2])-0.51), a) text(90, max((k[,2])-0.83), b) #

Dear R-listers, I am trying to compute simultaneous confidence intervals with simint from the package multcomp. 230 measures (abundance) have been taken in 23 sites (factor) of a data.frame (donnees: a file can be sent on request, saved with save(donnees,file="donnees")). I would like to get all pairwise comparisons with : mc<- simint(ren~ID,type="Tukey",data=donnees) I

Hi to all R lovers, I am thinking of investigating and perhaps improve the anova "aov" function in R. However before I do so I just like to see if anyone knows whether R is capable of multiple comparison tests (eg "Standard Difference of the Means, SED", "Fisher's Least Significant Difference" ...etc)? If it is, how can I do it? Thanks. Ko-Kang Wang

In the model: lm.1 <- lm(variable ~ BLOC + TIL * YEAR , data=selvanera) I found TIL*YEAR interaction significant. Then I am trying to compare means of the different levels of TIL inside every YEAR using: mc.2 <- glht(lm.1, linfct = mcp(TIL*YEAR="Tukey")) summary(mc.2, test = univariate()) but it does not work. There is any way of doing this, like the SLICE option in

I swear I've scoured the help files and several texts before posting what feels like a dumb newbie question. How can I draw two kernel density plots in the same frame ? I have similar variables in two separate data frames, and I would like to show their two histograms/densities in a single picture. Same units, scale, range for both, so I'm simply trying to draw one and then add the

Hello, I am a little out of date and am still using S-Plus instead of R. I haven't been able to find the right place to ask this question, so I thought I would ask it here and hope that someone can help. I am unable to located the TukeyHSD() function in S-Plus. I have been able to use the GUI to run Tukey's test, but the results don't provide a p-value. Does anyone know how I can

Dear List, are there post-hoc tests like Scheffe, LSD, etc. available after ANOVA test is performed with significant F-statistic? I have tried help.search("Scheffe"), but "No documentation found" (and I have most of packages installed). Probably there are such tests in R, and I am just searching badly... My second question is: Which test/method I should use for ANOVA-like

Hello, (maybe a quite basic statistical question, but I'm just struggling with it) I'm doing an anova: > Res1<-aov(H2O~location+topo+horizon+pF+Error(location:topo:horizon)) (water retention of soils) which gives a significant difference at factor "location". Which function should I use to now which locations (there are 3: A, B and C) differ significantly and which do