[R] Add column to the output of summary(glht).

I am trying to make the leap from an R users to an R aficionado . . . I am trying to understand how add a column to the output of summary (and to understand how summary() works). I have run a glmer fit0 <- glmer(Fall ~ Group+(1|PID),family=poisson(link="log"),data=data[data[,"Group"]!=0,])

Re: [R] Subscripting problem with is.na()

Yes, measurements below detection should be treated differently. I thought about the missing data issue, but there is another context in which spreadsheet data containing count data where 0 entries are deliberately left blank for readability or economy. In that case it is easier to import and

Re: [R] Add column to the output of summary(glht).

Hi you can check structure of any object by str so str(fit0) should give you an insight how fit0 is structured and str(SumTukey) gives you structure of SumTukey object (which is probably list). You can manipulate with the output as you wish according to the rules of R. Cheers Petr >

Re: [R] Add column to the output of summary(glht).

> On Jun 24, 2016, at 8:45 AM, John Sorkin wrote: > > > I am trying to make the leap from an R users to an R aficionado . . . > > I am trying to understand how add a column to the output of summary (and to > understand how summary() works). > > I have run a

Re: [R] Effect size measures for GLM

Gianfranco Lovison unipa.it> writes: > > Is there a library for (friendly) calculation of effect size measures for > Generalized Linear Models? I have found "compute.es", but it seems to be > suitable only for Linear Models. A library taking a glm object and > computing > partial R^2-type

Re: [R] Subscripting problem with is.na()

Hi I do not consider changing NA to 0 as a reasonable approach ( maybe only in some special case, which is not yours). rowSums(is.na(ds_temp1[,2:3])) [1] 1 1 2 2 0 0 2 2 1 2 1 2 1 2 1 2 2 2 2 2 gives you vector of numbers which is equal 2 only if they are both NA. So ds_temp1$open <-

[R] How to use seas()

Dear all, I am trying to run the seas() function. In doing so, I need an object of class "ts". I tried to generate an ts object using the ts() function but it does not work. Does anyone have an idea how to generate an ts object. In addition, I get the error that there are too many observations

Re: [R] Subscripting problem with is.na()

I would tend to agree. But NA is still preferable for both, no? -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Fri, Jun 24, 2016 at 8:42 AM, William

[R] Add column to the output of summary(glht).

I am trying to make the leap from an R users to an R aficionado . . . I am trying to understand how add a column to the output of summary (and to understand how summary() works). I have run a glmer fit0 <- glmer(Fall ~ Group+(1|PID),family=poisson(link="log"),data=data[data[,"Group"]!=0,])

Re: [R] Subscripting problem with is.na()

Is part of the issue that in common parlance "NA" or "N/A" may mean either "not available" or "not applicable" (e.g., isPregnant for a male) but in R NA means only "not available"? Bill Dunlap TIBCO Software wdunlap tibco.com On Fri, Jun 24, 2016 at 8:37 AM, Bert Gunter

Re: [R] Subscripting problem with is.na()

As Petr and Don have shown you, changing NA to 0 is unnecessary to get what you want. However, recoding to 0 may be OK, as NA has a specific meaning in this context, and you are just adding an extra code to a factor for a different level. But it still might cause you trouble later. One of R's

Re: [R] Subscripting problem with is.na()

See insert below. -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 On 6/24/16, 12:14 AM, "R-help on behalf of g.maub...@gmx.de" wrote: >Hi Bert, > >many thanks for all

Re: [R] How to do it in R

myfun <- function(a, b) a/(a + b) Leonardo Ferreira Fontenelle Em Sex 24 jun. 2016, às 13:05, André Luis Neves escreveu: > Dear all, > > I`ve got to calculate the ratio of methanogens to bacteria, but I > wouldn`t > like to divide the total copy numbers of methanogens ( on average 10^8) > by >

Re: [R] How to do it in R

pretty simple: > t_m <- 28e3 > t_b <- 710e3 > ratio <- t_m / (t_m + t_b) * 100 > ratio [1] 3.794038 Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Fri, Jun 24, 2016 at 12:05 PM, André Luis Neves

[R] How to do it in R

Dear all, I`ve got to calculate the ratio of methanogens to bacteria, but I wouldn`t like to divide the total copy numbers of methanogens ( on average 10^8) by bacteria (10^10) because they have different exponents and bases. So, my idea is to standardize both microorganisms counts to 10^3.

Re: [R] How to do it in R

This is very basic. Have you gone through any R tutorials? There are many good ones on the web. e.g., see here: https://www.rstudio.com/online-learning/#R In any case, you should not expect this list to teach you basic R. You *should* expect it to help you learn and improve your own efforts. I

[R] Fwd: Fwd: RE: Heatmap.2 Breaks argument

Dear Mr. or Mrs., Please see the previous messages for information. I had a problem with the heatmap.2 breaks argument and was wondering if someone could help me with that problem. I could not find information for the other authors of heatmap.2 to contact them. Could you provide help or

Re: [R] Add column to the output of summary(glht).

This won't make you an R aficionado, but depending on your needs library(broom) tidy(SumTukey) might be useful. This converts the output to a familiar data.frame, making it much easier to work with. Best, Ista On Jun 24, 2016 9:48 AM, "John Sorkin" wrote: > > I am

[R] Rattle

hi all, i am getting stuck while using Rattle - specifically, I am unable to get the graphics. i am using R 3.2.5 and some packages do not work in that. can someone pl help ? ** *Deva* ... *in search of knowledge, everyday something is added * *in search of

[R] R package updating

hi all, i notice that the R package available at CRAN is a more recent one compared to what I have. but i failed to update my machine and unfortunately, R does not work any more. can someone kindly suggest what is a way to update the R ver. is there a good way to retain all R files and update

Re: [R] Subscripting problem with is.na()

Hi Bert, many thanks for all your help and your comments. I learn at lot this way. My question was about is.na() at the first sight but the actual task looks like this: I have two variables in my customer data that signal if the customer accout was closed by master data management or by

Re: [R] Reading csv file with missing value

Read yesterday's and today's archives. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Fri, Jun 24, 2016 at 1:55 PM, Steven Yen

Re: [R] Rattle

This is like asking, "My car doesn't work. Can anyone tell me what is wrong?" Please spend some time reading (and paying attention to) the Posting Guide before sending any more emails here. -- Sent from my phone. Please excuse my brevity. On June 24, 2016 11:49:32 AM PDT, deva d

Re: [R] What's box() (exactly) doing?

Hi Jim, Thanks a lot, exactly what I was looking for. Cheers, Marius On Thu, Jun 23, 2016 at 11:06 PM, Jim Lemon wrote: > Hi Marius, > There are a few things that are happening here. First, the plot area > is not going to be the same as your x and y limits unless you

Re: [R] [FORGED] Re: Rattle

On 25/06/16 09:13, Jeff Newmiller wrote: This is like asking, "My car doesn't work. Can anyone tell me what is wrong?" Fortune nomination! cheers, Rolf -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276

Re: [R] What's box() (exactly) doing?

Try this one: myBox <- function (which = c("plot", "figure"), ...) { # draw filled rectangle where box() would draw open rectangle which <- match.arg(which) oldXpd <- par("xpd") on.exit(par(oldXpd)) if (which == "plot") { do.call("rect", c(as.list(par("usr")[c(1, 3, 2,

Re: [R] What's box() (exactly) doing?

Hi Jim, Here is a follow-up question: How would you replicate box("figure") (instead of box() = box("plot"))? I tried to fill the plotted box but there seems to be no argument to box("figure") that does that. If that's indeed the case, one could work again with rect() (thus replicating

[R] Reading csv file with missing value

I read a csv file (with read.csv) containing missing values (as shown below). Is there a convenient way to set these NA into zeros? Better yet, is there an option to assign zeros to these blank cells in reading the csv file? Thank you! NA -1 NA NA NA 1 NA NA NA NA NA NA NA NA NA

Re: [R] Fwd: Fwd: RE: Heatmap.2 Breaks argument

Did you try the maintainer() function? -- Sent from my phone. Please excuse my brevity. On June 24, 2016 10:45:07 AM PDT, fgoetz wrote: >Dear Mr. or Mrs., > >Please see the previous messages for information. I had a problem with >the heatmap.2 breaks argument and was

Re: [R-es] Ayuda sencilla (SQL)

Muchas gracias, Carlos. El 22 de junio de 2016, 21:14, Carlos Ortega escribió: > Hola, > > Hay una forma más sencilla utilizando "table()"... > > > datIn <- read.table("Datos.csv", sep = ";", header = TRUE) > > datIn <- na.omit(datIn) > > Cont_df <-

[R-es] Gráfico de dos Series en un mismo diagrama

Buenas Tardes, Deseo representar dos gráficos de series de tiempo en un mismo gráfico, utilizo plot(x) par(new=T) plot(y) Sin embargo, los ejes verticales manejan diferentes tamaños, mi pregunta es cómo puede hacer utilizando un eje diferente para cada serie pero que me queden en el mismo

Re: [R-es] Fwd: Ayuda ggplot2

Hola, library(reshape2) library(ggplot2) emp <- read.csv("tamano_empresas.csv", header=TRUE, sep=";", comment.char="" , strip.white=FALSE, dec = ",") melted = melt(emp, id.vars="Empresas") # Para lo que quieres esto te valdr� melted$variable<-substring(melted$variable, first = 2, last = 5)

Re: [R-es] Gráfico de dos Series en un mismo diagrama

Hola Elkin. No se si es exactamente lo que buscas, pero te cuento lo que yo hago.primero pinto la primera ajustando el eje 'y' al valor maximo de las dos series, pero no pinto los ejes. Pinto la segunda, y despues  pinto el eje de la primera en el margen izquierdo y el de la segunda en el