similar to: Accessing elements of plm outputs

Hi, I have the following problem (library Hmisc loaded, iris data loaded, R Version 1.7.0 (2003-04-16), packages updated, running on a linux Debian i386): > summary(Species~Sepal.Length,method="reverse")->a > a Descriptive Statistics by Species +------------+-----------------+-----------------+-----------------+ | |setosa |versicolor |virginica

Dear all Suppose the object below: > require(Hmisc) > require(plyr) > x <- dlply(iris, .(Species), describe) How can I print the object without displaying the attributes? I inspected ?print and ?print.default with no luck. > x $setosa x[, "Sepal.Length"] n missing unique Mean .05 .10 .25 .50 .75 50 0 15 5.006 4.40 4.59

Can someone get me going in the right direction with this test question? We have went as far as t.test in our class, but not to power or anova yet. I have pasted the question and a summary of the dataset "gpa". 5. Given that all the GPA data in this dataset is comprised of means, would you expect the underlying dataset of students? individual GPA to be more variable or less variable

Hello, I have a question about glm: if i have a binary covariate (unit=1,0) the reference group would be 0? (prediction for unit=1) example: dat1<-data.frame(y,unit,x1,x2) log_u <- glm(y~.,family=binomial,data=dat1) summary(log_u) Estimate Std. Error z value Pr(>|z|) (Intercept) -0.54247 0.24658 -2.200 0.0278 * unit1 -0.13052 0.22861 -0.571 0.5680 aps

Hii, I will create boxplots from matrices. I have the following data sets: 5.0 1.78 2.99 2.019 0 10.0 1.79 3.00 1.744 0 15.0 1.78 2.98 1.936 0 20.0 1.78 2.99 1.975 0 25.0 1.73 2.91 3.591 0 30.0 1.79 3.00 1.966 0 35.0 1.79 3.00 2.451 0 40.0 1.79 3.00

I am running 4.8-STABLE. I recently did a cvsup (repeatedly, most recently today). Ever since, portupgrade dies with the following messages: zircon# portupgrade -r pkgconfig /usr/local/lib/ruby/site_ruby/1.6/pkgdb.rb:310:in `deorigin': failed to convert nil into String (PkgDB::DBError) from /usr/local/lib/ruby/site_ruby/1.6/pkgdb.rb:903:in `tsort_build' from

I have a file, each column of which is a separate year, and each row of each column is mean precipitation for that month. Looks like this (except it goes back to 1964). month X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007 X2008 X2009 1 1.600 1.010 4.320 2.110 0.925 3.275 3.460 0.675 1.315 2.920 2 2.960 3.905 3.230 2.380 2.720 1.880 2.430 1.380

Hi, My name is Jos?? Mar??a G??mez, and I am pretty new in R. Thus, I apologize deeply if my questions are extremmely na??ve.I have checked several available books and URL's, without finding any answer. I'm trying to fit Generalized Linear Mixed Models via PQL. Below I provide the structure of my data set. Year and Plot are random variables. Fate is the binomial dependent. I have severe

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

Dear R users, please can you help me understand the behaviour of abline using function lm. I'm trying to learn linearity over groups. So I make three groups with 10 values each: test=data.frame(cbind( l=c(rnorm(10,0,30),rnorm(10,100,30),rnorm(10,200,30)), t = c(rep(0,10), rep(1,10), rep(2,10)) )) when I do: plot(test$l~test$t) abline(lm(test$l~test$t)) the abline is a straight line

Estimados apañeRos: La duda o propuesta que os voy a plantear es a la vez metodológica y relacionada con R. Me encuentro trabajando con tres variables que son el resultado de un computo de porcentajes. Me explico, se toma una muestra de n casos (unos 6.500 aprox) pertenecientes a i individuos (unos 230 aprox) en la que se comprueba si un determinado evento ha ocurrido o no, anotándose 1 en

Dear Professor Lumley I am relatively new to using R and also to logistic regression. We have analysed our Dudley Health Survey using the survey package. I am now trying to look at associations using svyglm but I am unsure of how to interpret the output and present the resulting model or whether there are any other things I should do to check the validity of the model. Below is an example of

Hola Isidro: También había sopesado esa posibilidad pero por una parte me parecía "complicar" el proceso y por otra tengo mis dudas acerca de que en que el análisis cluster pueda prescindir de la hipótesis de Normalidad en las dos variables "raras" que tengo. De cualquier forma muchas gracias por tu idea. Un saludo {In Archive} RE: [R-es] Clasificacion de

En cualquier caso, para nada necesitas normalidad en las variables. De las distribuciones que presentas, y abundando en lo que te indica Olivier, ¿sería interesante considerar sólo 2 grupos (?normal? y ?enhorabuena?) en las variables 2 y 3? Todo esto sin saber de lo que estamos tratando, claro? XD Un saludo, Isidro De: Olivier Nuñez [mailto:onunez en unex.es] Enviado el: martes, 30 de

Dear R-help This is a follow-up to my previous post here: http://groups.google.com/group/r-help-archive/browse_thread/thread/d9b6f87ce06a9fb7/e9be30a4688f239c?lnk=gst&q=dobomode#e9be30a4688f239c I am working on developing an open-source automated system for running batch-regressions on very large datasets. In my previous post, I posed the question of obtaining VIF's from the output of

Hello. We are planning a migration from Sun One Messaging Server to Dovecot+Postfix+Ldap. We are using Dovecot 1.2.11 with Sun One Directory Server 5.2 ldap (we will migrate to Directory Server 6.3.1) soon. In our University we have 65.000 students, 5.500 staff and 6.500 teachers. Our main ldap realm is dc=us,dc=es (us means University of Seville, Spain). We have two e-mail

Hello all! I have a problem to arrange data in another form. My initial data is like this: 'data.frame': 421 obs. of 58 variables: $ 01A: num NA NA NA NA NA NA NA NA NA NA ... $ 01B: num NA NA NA NA NA NA NA NA NA NA ... $ 03A: num NA NA NA NA NA NA NA NA NA NA ... $ 03B: num NA NA NA NA NA NA NA NA NA NA ... $ 05A: num NA NA NA NA NA NA NA NA NA NA ... $ 05B: num NA NA NA NA

As nobody seems to answer my request, I simply post it again. Is there any reason why the r-help-mailinglist should not be converted to a newsgroup? These were the advantages of a newsgroup I mentioned earlier: -) you can easily search the archives -) the discussion is faster (I experience that the R-mailinglist has a lag of about 1 to 2 hours (not for everyone!!!). When I ask a question, I get