similar to: Sums of sq in car package Anova function

Hello all, As a general programming question I can't seem to figure out how to make a list of lists in R. As matrices won't work as they have to be rectangular. I am sure that there is an easy solution but... the specific situation is this: - I have created a Tukey confidence interval table and have listed the means that are not significantly different - then using these not

Hi, I need construct box plot graph, but I want keep Groups order karla = data.frame( Groups = factor(rep(c('CPre','SPre','C7','S7','C14','S14','C21','S21'), 11)), Time = rep(c(0,7,14,21), 11), Resp = valor ) boxplot(Resp~Groups, order=T) doesn't work. How do this? -------------------------------------- Silvano Cesar da

I have been trying to install the new, patched R 2.1.0a along with the Rcmdr on my mac. I got the basic R package installed, but in Project Manager, the tcltk package will not install. When I downloaded the Rcmdr binary package from the Berkeley site, it did not show up on the package installer list until a second attempt, but it still does not show up on the package manager list. On

Buenas tardes, Quiero aplicar la función rbind y necesito tener los mismos nombres de columnas. Como tengo unas 195 variables en cada dataframe, necesito hacerlo de una forma rápida. Tengo 9 bases de datos y tengo que fusionar todas. ¿Como puedo comprobar que los nombres de las variables son los mismos? Y de lo contrario, ¿como detecto las diferencias? He probado con

------------------------------------------------------------------------------------------------ CALL FOR PAPERS 15th International Conference on Managed Languages & Runtimes (ManLang'18) September 10-14, 2018, Linz, Austria http://ssw.jku.at/manlang18/ ------------------------------------------------------------------------------------------------ ManLang (formerly PPPJ) is a

Hi there, I have 2 questions which I cannot find answers for in Dalgaard or the helpfiles currently available. (1) I wish to plot 3 ecdf plots on one graph. Is there any way of holding a figure to plot all on the same plot. I can't get an ecdf for any plot command other than "plot" or "ecdf.plot" so using lower level command is no good, have also tried x11(). (2) If

(2) If I have three curves on a plot how can I label each curve so the > label sits directly above the curve? [A]As to your second question: I don't think there will be a simple answer because on has to decide where "above" is. Could you please be more specific? Dietrich I have a picture of S+ output which has the label sitting on top of each curve and curving with it

Hi, I would like to plot a graph which sits in the background as a watermark with other plots in the foreground - on top. I have looked through the threads on the r-project website but they seem to concern background colours rather than actual background plots. I have also searched through the demo's and pars but can't find any eg's any ideas - I only want to plot a simple graph

Hi, I just wanted to plot a boxplot with a nice curve going through it, I thought this would be a simple task but for some reason I can't get the two graphs on the same page accurately. Enclosed is the code showing the two plots seperately and together. I would have thought it should work if I could use boxplot() then plot() overlayed but it won't allow the command add=TRUE (which

on hp-ux 11 i see: $ date;ssh jenny Fri Oct 12 14:44:13 PDT 2001 Last successful login for stevesk: Fri Oct 12 10:45:42 PST8PDT 2001 on pts/2 Last unsuccessful login for stevesk: Mon Sep 24 22:55:53 PST8PDT 2001 Last login: Fri Oct 12 10:45:43 2001 from 172.31.1.53 You have mail. so solaris PAM is different. can other solaris+PAM users confirm this? On Fri, 12 Oct 2001, Benn Oshrin wrote:

Dear Person That Posted This Email, Please be aware that calls for papers are only on-topic for this email list if they are specifically about LLVM. For example, a call for a workshop on LLVM in Scientific Computing would be on topic whereas a general compiler call for papers (like the one below) is not. That said, I hope you receive many strong submissions. Your conference looks

Hello dear R members. I have been learning the Anova syntax in order to perform an SS type III Anova with repeated measures designs (thank you Prof. John Fox!) And another question came up: where/what are the (between/within) residuals for my model? ############ Play code: phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),

On the Squirrelmail website I setup in the FAQ a quick install guid on how to get squirrelmail to work with Dovecot. Maybe someone can verify my setup procedure really quick. Maybe someone thinks this is wrong. Also if you have something to add please do so. No one has bothered to do this before so I took it upon my self. Here is the link in the FAQ its about 8 lines down

Hi, I have been trying to do a simple GLM with count data using a poisson distribution. The results show evidence of underdispersion. I have only ever encountered overdispersion. Am I still able to use family=quasipoisson to correct for underdispersion? Thank you, Karla -- View this message in context:

Hi, I am trying to test the random labeling hypothesis with my data. I have a point pattern of bird nests and each nest has a mark that it is either predated or not-predated. I want to test the null hypothesis that the K distribution of predation events is equal to the K distribution of all events (nests). This is the code I am using: Kdif <- function(X,...,i) { Kidot <-

I've been using the nlme library to construct lme models and I'd like to extract an r-sq value for the model - how do I do this? With glm I just used the reported (null deviance-deviance)/(null deviance). Deviance is not reported directly for lme nor can it be extracted using the deviance() function. How do I work out the deviance or can I work out R-sq in some other way? Thanks, Yvonne

Dear R-Users, I have a trivial problem, but extensive googling and ??'ing did not solve it: I want to obtain the sums of squares from a summary.aov() object for later use. Example: > DV <- rnorm(100) > IV1 <- as.factor(rep(c("male", "female"), times = 50)) > IV2 <- as.factor(rep(c("young", "old"), times = 50)) > >

Hi, On Tue, 2023-03-07 at 09:49 +0800, Xuan Zhuo wrote: > On Mon, 6 Mar 2023 12:58:22 -0500, "Michael S. Tsirkin" <mst at redhat.com> wrote: > > On Mon, Mar 06, 2023 at 12:15:33PM +0800, Xuan Zhuo wrote: > > > If the queue of xdp xmit is not an independent queue, then when the xdp > > > xmit used all the desc, the xmit from the __dev_queue_xmit() may

On Wed, Mar 08, 2023 at 04:13:12PM +0800, Yunsheng Lin wrote: > On 2023/3/8 15:14, Xuan Zhuo wrote: > > On Wed, 8 Mar 2023 14:59:36 +0800, Yunsheng Lin <linyunsheng at huawei.com> wrote: > >> On 2023/3/8 10:49, Xuan Zhuo wrote: > >>> If the queue of xdp xmit is not an independent queue, then when the xdp > >>> xmit used all the desc, the xmit from

Hi I spent hours looking over my formula. Somehow I cant find the reason why it gives me different answer. help appreciated. x = as.matrix(read.table("http://www.niehs.nih.gov/research/atniehs/core/microarrays/docs/heinloth.txt",1)) x = t(x) #now rows are subjects, cols are genes x = x[order(rownames(x)),] #order by treatment group oxygen, ultra-violet, gamma radiation y =