similar to: Trouble with performing post hoc analysis Tukey for lme model using ghlt

I'm not an experienced user of the 'nlme' package but I ran into a curious problem recently that I couldn't find on the mailing list anywhere. I have simple function which calls 'groupedData' and 'nlme': g <- function(dataset) { unig <- groupedData(logfv ~ Vs+v1+v2+v3+PHA| EQ, data=dataset) modelfv <- nlme(logfv~a + c*log(Vs/a1)+ b1*( v1 +

Dear all, I am trying to get a PERMANOVA with quite large data set. I am reading a lot about this question, but I do not get the answer about it. Although I know that the R function is adonis () (vegan package), it does not work: adonis(Pha.env~SPha, data=Pha, permutations=10) The error message: Error: cannot allocate vector of size 334.2 Mb In addition: Warning messages: 1: In vegdist(lhs,

Dear R-group, I am trying desperately to get this Tukey test working. Its my first time here so I hope my question is not too stupid, but I couldn't find anything helpful in the help or in the forum. I am analysing a dataset of treated grasses. I want to find out, if the grasses from different Provenances react differently. In the aov test I found a significance for the combination Treatment

Hello, i have installed samba alpha3, domain-join works fine. But smbpython setup/adduser is buggy, so i tried to install newest checkout. (I also tried it with a new installed machine) After configure, make, and make install there is no smbpython in source/bin. Everyother seems to bee there. It is a ubuntu with the packets installed mentioned in the tutorial plus python and python-dev and

Dear R-help list, I would like to perform multiple comparisons for lme. Can you report to me if my way to is correct or not? Please, note that I am not nor a statistician nor a mathematician, so, some understandings are sometimes quite hard for me. According to the previous helps on the topic in R-help list May 2003 (please, see Torsten Hothorn advices) and books such as Venables &

Dear list, As a reply to my recent mail: > simint and TukeyHSD work for aov objects. > Can someone point me to similar functions for lme objects? Douglas Bates wrote There aren't multiple comparison methods for lme objects because it is not clear how to do multiple comparisons for these. I don't think the theory of multiple comparisons extends easily to lme models. One could

Hi Folks - I'm experiencing a very strange problem with Server 2008 machines (for all intents and purposes related to Samba, it's Vista) connecting to a Samba Server. The Samba machine is a RHEL4.6 machine running Samba 3.0.25b. I am joined to the mixed mode AD domain via the command "net rpc join -U administrator <domain>". I am not using winbind or kerberos. Or

Dear all, I've discovered the possibility to do the Nemenyi-Damico-Wolfe-Dunn test in the library(coin); oneway_test() With the given example I am unfortunately not able to reproduce the test. What does trafo and contrMat mean? I have a dataframe with 176 elements in 7 classes. It may be a problem, that my model isn't balanced? In Class1 I happen to have 4 elements, while there are 90 in

On Tue, Mar 6, 2018 at 10:22 PM, Paul Anderson <pha at umich.edu> wrote: > Raghavendra, > > I've commited my tests case to https://github.com/powool/gluster.git - > it's grungy, and a work in progress, but I am happy to take change > suggestions, especially if it will save folks significant time. > > For the rest, I'll reply inline below... > > On Mon,

Dear R-Users, I'm trying to set up a repeated measures anova with two within subjects factors. I tried it by 3 different anova functions: aov, Anova (from car package) and lme (from nlme package). I managed to get the same results with aov and Anova, but the results that I get from lme are slightly different and I don't figure out why. I guess I did not set up the error structure

I suppose that once you've got a 6502 working, adding support for a 4510 shouldn't be too difficult.... (http://c65gs.blogspot.com.au/) On Fri, Jul 4, 2014 at 10:38 AM, Bruce Hoult <bruce at hoult.org> wrote: > On Fri, Jul 4, 2014 at 12:02 PM, Edwin Amsler <edwinguy at gmail.com> wrote: > >> Well, the stack pointer be a single byte, so pushing things on there

A couple days ago, Mark Leeds asked about a solution that would basically stagger two lists, a and b, to return a list in the form of a[1], b[1], a[2], b[2], a[3].... In particular, the summary of his question was in reference to lists defined by x <- 5 tempin <- seq(1,1411, by=30) a <- tempin b <- tempin + x I offered the following function everyOther <-

On Tue, Mar 6, 2018 at 10:58 PM, Raghavendra Gowdappa <rgowdapp at redhat.com> wrote: > > > On Tue, Mar 6, 2018 at 10:22 PM, Paul Anderson <pha at umich.edu> wrote: > >> Raghavendra, >> >> I've commited my tests case to https://github.com/powool/gluster.git - >> it's grungy, and a work in progress, but I am happy to take change >>

I just modified mail-filter-out.sh like this: cat &nbsp;mail-filter-out.sh &nbsp; #!/bin/sh &nbsp; USER=$1 &nbsp; cat &gt; /tmp/tempfile &nbsp; cat /tmp/tempfile cat &nbsp;/tmp/tempfile This message is in MIME format. --nsmail-7bjyg3nr5o-7bmibqhut9 Content-Type: text/html; charset="UTF-8" Content-Transfer-Encoding: base64

Hello everyone I am currently trying to conduct analysis of my graduate thesis data using a mixed effects model and I have reached an impass. When I try to conduct a multiple comparison, I get an error (See below): > fm3<- lme(abovegroundbiomass.m.2~medium*amelioration*fertilizer*treatment, random=~1|block/medium/amelioration/fertilizer) > tukeytest<-glht(fm3,

I have run the following logistic regression model: options(contrasts=c("contr.treatment", "contr.poly")) m <- glm(wolf.cross ~ null.cross + feature, family = "binomial") where: wolf.cross = likelihood of wolves crossing a linear feature null.cross = proportion of random paths that crossed a linear feature feature = CATEGORY of linear feature with 5 levels:

I am having a hard time interpreting the phase and coherency plots. x is two timeseries that occur at the same time i.e. a b 1 11.2 12.3 16 11.3 12.4 31 11.4 12.5 46 11.5 12.6 ...etc even though my example is does not show this they are oscillating at more or less the same frequency just shifted by t=x (imagine two sine waves offset with the 2nd sine

Hi, I have a fairly simple repeated measures-type data set I've been attempting to analyze using the lme function in the nlme package. Repeated searches here and other places lead me to believe I have specified my model correctly. However, I am having trouble with post-hoc tests. From what I gather, other people are successfully using the glht function from the multcomp package to

+Csaba. On Tue, Mar 6, 2018 at 2:52 AM, Paul Anderson <pha at umich.edu> wrote: > Raghavendra, > > Thanks very much for your reply. > > I fixed our data corruption problem by disabling the volume > performance.write-behind flag as you suggested, and simultaneously > disabling caching in my client side mount command. > Good to know it worked. Can you give us the

Hi, the following R lines work fine in R 2.4.0, but not in R 2.4.1 or any devel versions of R 2.5.0 (see below for details). library(drc) # to load the dataset 'PestSci' library(nlme) ## Setting starting values sv <- c(0.43355869, 2.49963220, 0.05861799, 1.73290589, 0.38153146, 0.24316978) ## No error m1 <- nlme(SLOPE ~ c + (d-c)/(1+exp(b*(log(DOSE)-log(e)))), fixed =