similar to: How to code mixed model with nested factors in lmer

I am working with amino acid sequences changing each letter to numbers.I have a data from acf transformation called Zm as shown below. I would like to get Indices D1 to D10 and then create F1 to F10 as indicated below. Is there anyway I can do that in R without typing each of them one by one. For example to get D1 to D3, I have to type D1<-c(Zm[[1]])[1][[1]], D2<-c(Zm[[2]])[1][[1]], and

[mod: Forwarded by Jeff Uphoff. I tried to mangle the headers that it appears as the original post: with an invalid return address. -- REW] A few months back, a very alpha version of bliss got posted. That shouldn''t have happened, but, it was pretty much ignored so I didn''t worry about it. But now it seems there''s a bit of a fuss about this. I''ll post the

Hi list, I was testing out the "doSMP" package from REvolutionR in my regular R2.11.0 installation and I got the following error message.? Well, one obvious thing is that R2.11.0 was built using "i386-pc-mingw32" which is different from what revoIPC used.? I could just use REvolutionR, but all my R peripherals were set up to work with the regular R2.11.0.? So, I really want

Hola, Lo hice de la forma en que propones y no me sale, no sé si algo esté mal con R, porque no me respeta ni el ?key.space?, me lo sigue poniendo en la parte izquierda del gráfico Saludos El 06/05/2014, a las 02:21, Olivier Nuñez <onunez en unex.es> escribió: > El segundo argumento en spplot ha de definir el número/nombre del intervalo al que > pertenece el area a pintar. >

Signed-off-by: Scott Seago <sseago at redhat.com> --- AUTHORS | 17 ++++++ README | 10 +++ conf/ovirt-agent | 12 ++++ conf/ovirt-db-omatic | 12 ++++ conf/ovirt-host-browser | 12 ++++

Hola, El problema con la propuesta de Olivier es que los intervalos son diferentes para cada variable. La forma sencilla y rápida es: spplot(zm["part88"], col.regions=plotclr, at=class$brks) Pero para que quede más elegante hay que dar algunos pasos más: ## Intervalos en forma character op <- options(digits=4) tab <- print(class) options(op) intChar <- names(tab) ## Indice

Hi Aristide and welcome to R-help! Your message was a bit mangled [*]. It's best to compose messages to this mailing list in plain text. Otherwise (when composed in HTML), the mailing list eats the HTML part and we're left with the plain text part automatically generated by your mailer, which isn't always readable. ? Wed, 30 Oct 2024 17:45:29 +0100 Elys?e Aristide <ariel92and at

Hi Ivan, Thank you for your message. Does that mean that I should send a new message? Or is it okay for this time? Best, Aristide On Fri, Nov 1, 2024, 22:29 Ivan Krylov <ikrylov at disroot.org> wrote: > Hi Aristide and welcome to R-help! > > Your message was a bit mangled [*]. It's best to compose messages to > this mailing list in plain text. Otherwise (when composed in

Hi everyone, I am using the gmm function from the gmm package and encountered an unexpected error. No model can be estimated if I load formula.tools?I need to restart R each time. Here is a simple reproducible example: *library(gmm)data(Finance)r <- Finance[1:300, 1:10]rm <- Finance[1:300, "rm"]rf <- Finance[1:300, "rf"]z <- as.matrix(r-rf)zm

We are hiring an R programmer to make biologists better at crop improvement: http://jobs.monsanto.com/missouri/research-and-development/jobid6130734-r-programmer-jobs Please apply if you are passionate about building R culture and infrastructure. . . .. . ... . ... . .. . . ... . ... . . . .. .. .... .. . ... . Barrett Foat, PhD Genome Data Analytics Team Lead

Hola, no me funciona: Te contesto lo mismo que le contesté a otro compañero de la lista, espero ser más claro. Les comparto la carpeta con la capa y los datos que estoy usando. https://www.dropbox.com/s/umy0evt3qm1wr4d/dissolve.zip Esto es lo que estoy haciendo: library(rgdal) library(maptools) library(sp) library(classInt) library(RColorBrewer) zm <- readOGR(?.?, ?zmdis?) data <-

Hi all, I have a question regarding how to properly analyze a data set and then how to perform the analysis in R. First, I have data that I would like to analyze using a mixed GLM (I think this is the most appropriate method, but I am unsure). In a mixed model (y = X*beta+Z*gamma+epsilon), I would like to structure the variance matrices of gamma, G, and the error, R, to take advantage of all my

Hola, Intento representar en un mapa participaciones porcentuales de los sectores económicos y no logro hacerlo con spplot. He intentado con spplot(zm, c(?part88?, ?part93?), cuts=4, col.regions=brewer.pal(4, ?Set3?)) Lo pude hacer utilizando el base graphics de R, definiendo: > plotvar88 <- zm$part88 > nclr <- 8 > plotclr <- brewer.pal(nclr, "PuOr") > plotclr

Hi, I don't know about the lazyeval package or what you are trying to do but to answer the main question "How to create columns dynamically using a function?" I would do something like that: # dataset dem <- structure(list(id = c("L1", "L2", "L3", "M1", "M2", "M3"), TEST_SET_NAME = c("A", "A",

Dear guRus: I am using lmer for a mixed model that includes a random intercept for a set of effects that have the same distribution, Normal(0, sig2b). This set of effects is of variable size, so I am using an as.formula statement to create the formula for lmer. For example, if the set of random effects has dimension 8, then the lmer call is: Zs<-

Your code contains the lines stopifnot(!(any(data1$norm_sd >= 1))) if (!(any(data1$norm_sd >= 1))) { df1 <- dat1 return(df1) } stop() "throws an error", causing the current function and all functions in the call stack to abort and return nothing. It does not mean to stop now and return a result. Does the function give

I am simulating an ANOVA model and get a strange behavior from the summary function. To be more specific: please run the following code and see for yourself: the summary()[["r.squared"]] values of two identical models are quite different!! ## 3 x 3 ANOVA of two factors x and z on outcome y s.size <- 300 # the sample size p.z <- c(0.25, 0.5, 0.25) # the probabilities of factor z ##

Hi, I am writing a function to dynamically create column names and fill those columns with some basic calculations. My function "demo_fn" takes argument "blup_datacut" and I like to use the contents of those arguments to dynamically create new columns in my dataset. Please note that I have another function called "calc_gg" within the function "demo_fn".

hello I make a subroutine that give-me a (mathematical) function in string format. I would like transform this string into function ( R function ). thanks for any tips. cleber #e.g. fun_String = "-100*x1 + 0*x2 + 100*x3" fun <- function(x1,x2,x3){ return( ############ evaluation( fun_String ) ############ ) True String mathematical function :-( :-( > nomes [1]

Inthe code below I was trying to to obtain the GMM estimates for CAPM (REGRESSION) for 36 stocks each have 180 observations,however it only gives me one output rather than 36. In SAS i would just put in a *By statement*. I have a variable TICKER that categorize them into 36 groups. *How can I obtain all 36 output instead of just one.* **