similar to: anova comparisons

Hello list, I'm doing a bootstrap analysis where some models occasionally fail to converge. I'd like to automate the process of restricting AIC to the models that do converge. A contrived example of what I'd like to do is below: resp <- c(1,1,2) pred <- c(1,2,3) m1 <- lm(resp~pred) m2 <- lm(resp~poly(pred,2)) m3 <- lm(resp~poly(pred,3)) # Fails, obviously ## Some

Hi, I am analyzing a repeated-measures Greco-Latin Square with the aov command. I am using aov to calculate the MSs and then picking by hand the appropriate neumerator and denominator terms for the F tests. The data are the following: responseFinger mapping.code Subject.n index middle ring little ---------------------------------------------------------------------------- 1 1

Hi, I have some doubt about how qualitative factors are coded in R. For instance, I consider a response y, a quantitative factor x and a qualitative factor m at 3 levels, generated as follow : y_c(6,4,2.3,5,3.5,4,1.,8.5,4.3,5.6,2.3,4.1,2.5,8.4,7.4) x_c(3,1,3,1,2,1,4,5,1,3,4,2,5,4,3) m_gl(3,5) lm(y~x+m) Coefficients: (Intercept) x m2 m3 3.96364 0.09818

Hi, I have a number of genes (columns) for which I want to examine pairwise associations of genotypes (each row is an individual)...For example (see data below), I would like to compare M1 to M2, M2 to M3, and M1 to M3 (i.e. does ac from M1 tend to be found with bc from M2 more often than expected.) Down stream I will be performing chi square tests for each pair. But I am looking for a way to

I am analyzing from a very simple experiment. I have measured plants of two different colours (yellow and purple) in 9 different populations. So, I have two different factors : a fixed effect (Colour with two levels) and a random one (Population with 9 levels). I first analyzed the data with the aov function LargS is the variable aov(formula = LargS ~ Col + Error(Col/Pop)) Terms:

Kevin Wright wrote: > library(nlme) > m2 <- gnls(conc ~ t1*(1-t2*exp(-k*time)), > data = df.Chloride, > start = list( > t1 = 35, > t2 = 0.91, > k = 0.22)) So my error was to use nls instead that gnls. Thanks a lot, Kevin. > summary(m2) > plot(m2) > lag.plot(resid(m2), do.lines=FALSE) >

I would like to perform partial mantel tests on only within group values, with "between group" values assigned to NA. This is possible in package ncf partial.mantel.test, however this sues a different permutation to that used in ecodist.ecodist will not accept data with NA values, returning a "matrix is not square error. is it possible to perform this test in ecodist? many thanks

Hi, I have a data-matrix: > PID sato hrs fim health 214 3 4.376430 6.582958 5 193 6 4.361825 3.138525 6 8441 6 4.205771 3.835886 7 7525 6 4.284489 3.245139 6 6806 7 4.168926 2.821833 7 5682 7 1.788707 1.212653 7 5225 6 1.651463 1.436980 7 4845 6 1.692710 1.267359 4 4552 5 1.686448 1.220539 6

Dear All, I am new for R and SEM. I try to fit the model with Y (ordinal outcome), X (4 categorical data), M1-M3 (continuous), and 2 covariates (Age&sex) as a diagram. library(polycor) model.ly <-specify.model() 1: x -> m1, gam11, NA 2: x -> m2, gam12, NA 3: x -> m3, gam13, NA 4: age -> m1, gam14, NA 5: age -> m2, gam15, NA 6: age -> m3, gam16, NA 7: sex -> m1,

Hi all, I have a matrix simulating migration in a spatial model. I want to be able to define movement (the values of m1, m2 and m3) as only != 0 between adjacent patches and contingent on certain conditions as defined in the function. Here is the code: WET<-function(t) {everglades$precipitation[t]} #simply reads precipitation data from a csv, value is numeric AB<-function(WET,t)

On Jul 5, 2013, at 9:50 AM, Stéphane Letz <letz at grame.fr> wrote: > > Le 5 juil. 2013 à 04:11, Tobias Grosser <tobias at grosser.es> a écrit : > >> On 07/04/2013 01:39 PM, Stéphane Letz wrote: >>> Hi, >>> >>> Our DSL can generate C or directly generate LLVM IR. With LLVM 3.3, we can vectorize the C produced code using clang with -O3, or

Le 5 juil. 2013 à 04:11, Tobias Grosser <tobias at grosser.es> a écrit : > On 07/04/2013 01:39 PM, Stéphane Letz wrote: >> Hi, >> >> Our DSL can generate C or directly generate LLVM IR. With LLVM 3.3, we can vectorize the C produced code using clang with -O3, or clang with -O1 then opt -O3 -vectorize-loops. But the same program generating LLVM IR version cannot be

Hi I have a small ISDN PABX , (BRI) at home - Siemens Gigaset 4175 which has cordless DECT extensions. I have set up asterisk on FC3 with two HFC cards and I am using the latest bristuff. I am trying to use * between my ISDN line and my PABX. One card is in TE mode and can receive and make calls OK. I can make and receive external calls to a SIP phone. I cannot get the other card to talk to my

> On Jun 6, 2017, at 4:01 AM, Jim Lemon <drjimlemon at gmail.com> wrote: > > Hi Bogdan, > Kinda messy, but: > > N <- data.frame(N=c("n1","n2","n3","n4")) > M <- data.frame(M=c("m1","m2","m3","m4","m5")) > C <-

Thank you David. Using xtabs operation simplifies the code very much, many thanks ;) On Tue, Jun 6, 2017 at 7:44 AM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Jun 6, 2017, at 4:01 AM, Jim Lemon <drjimlemon at gmail.com> wrote: > > > > Hi Bogdan, > > Kinda messy, but: > > > > N <-

Any ideas? Sangoma support is closed for the evening. I have the latest Sangoma drivers and Asterisk 1.4 everything installed. When I fire up asterisk, I keep getting "Primary D-Channel on span 1 up" repeated over and over. The B channels never come up. There are no errors in any of the logs, zttool, or the wanpipe tools. Intense pri debug output: < Unnumbered frame: < SAPI:

Saludos, tengo un problema que no puedo resolver dentro del R Estoy creando una funcion que a partir de un objeto compuesto por diferentes matrices como el que esta a continuacion: [[1]] M1 M2 M3 M4 sp1 2 0 1 8 sp2 4 5 2 4 sp3 0 0 4 0 sp4 5 7 5 0 sp5 0 4 0 0 [[2]] M3 M2 M4 M1 sp1 1 0 8 2 sp2 2 5 4 4 sp3 4 0 0 0 sp4 5 7 0 5 sp5 0 4 0 0 [[3]]

Hello, My name is Ick Hoon Jin and I am Ph. D. student in Texas A & M Univ.. When I run the C embedded in R in the Linux system, I confront the following error after 6,000 iteration. By googling I found this error is from the problem in C. *** caught segfault *** address (nil), cause 'memory not mapped' My C code is following:

Dear Contributors I have a problem with the collection of data from the results of a test. I need to perform a comparative test over groups of data , recall the value of the pvalue and create a table. My problem is in the way to replicate the analysis over and over again over subsets of data according to a condition. I have this database, called y: gg t1 t2 d 40 1 1

Here's another approach: N <- data.frame(N=c("n1","n2","n3","n4")) M <- data.frame(M=c("m1","m2","m3","m4","m5")) C <- data.frame(n=c("n1","n2","n3"), m=c("m1","m1","m3"), I=c(100,300,400)) # Rebuild the factors using M and N C$m <-