Re: [R] interpreting repeated measures mixed model

nobody? please - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/interpreting-repeated-measures-mixed-model-tp4647652p4647689.html Sent fr

[R] correct use with fligner test

Hi R users, Iwant to run a model, but want to know previously if my data present homoscedascity: the model is something like this: glm(ax~a*b*c*d) I´ve found four different ways to test for homoscedascity, but i want to know which is correct: 1.- > fligner.test(glm(x~a*b*c*d)) 2.- > fligner

[R] Don´t know what test i have to use

Hello, I´m starting with my PhD and I have to stop because i got a little knowledge in R and statistics. I´ve got a model of this kind: binary response variable: prevalence of infection (0/1) 3 categorical independent variables: sex, month and name of the area I was trying with a full model li

Re: [R] Don´t know what test i have to use

Thanks very much both. I´m starting playing with it, i was a little afaid because it was part of my job, but now i've found it very funny. Josh, I've got just data for 3 representatives months, and it's not a priori rejectable that could be differences in the ratio of changes along the months bet

Re: [R] Don´t know what test i have to use

Thanks for the advice, I will use next years. Till know i´ve just got data for 3 independent months and one of the months it´s the joining for all the summer because of the small sample size, so, I suppose, I can't use it in the way you say. - Mario Garrido Escudero PhD student Dpto. de Biolo

[R] send commands from script to console

Hi, I'm working with R in windows and I wonder if there is any command (of the kind CTRL+ ) to transfer the commands I've worked with (like: > model<-glm(prevalencia~edadysexo*mes*zona,binomial)) to a script automatically, without the results I received from R after execute them. Another question,

Re: [R] send commands from script to console

Thanks very much - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/send-commands-from-script-to-console-tp3216553p3216780.html Sent from t

[R] Problems with TeachingDemos package

Hi, somebody can help me. Idon't understand very well the manual of the this package, I suppose I must enter "wdtxtStart", before I start working, but is supposed that this command should open Word, but nothing happens. I use Rgui, not Rcmdr. Thanks - Mario Garrido Escudero PhD student Dpto.

Re: [R] Problems with TeachingDemos package

R2wd is working but i received an alarm: > wdtxtStart() Error en R2wd::wdGet() : tentativa de aplicar una no-función The translation is "attempt to apply a no-function" - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Univ

[R] Applying a function to a subset

Here´s my problem, i have developed the function kurtosis (using R-book as a guide) with this commands: kurtosis<-function(x) { m4<-sum((x-mean(x))^4)/length(x) s4<-var(x)^2 m4/s4 - 3 } Then create the object fem, which is the difference between the count of a trait in the left side of the body mi

Re: [R] Applying a function to a subset

That's it. thanks very much, I save hours if work qith your reply - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/Applying-a-function-to-a

[R] Problems making subsets with [] or "subset"

Hi, I'm trying to make a model in order to know wich factors got´s influence in the intensity of a infection, but just in the individuals who's got this infection. In my data I've got a variable called "prevalence" with 2 levels: 1.- Infected individual 0.- Non infected So what i'm trying to do i

[R] interpreting interactions in a model

Hi, I´ve got this model > model<-glm(prevalence~agesex+agesex:month,binomial) and the output of anova is like that > anova(model,test="Chisq") Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL524 206.97 agesex

Re: [R] interpreting interactions in a model

Thanks, I´ve found the differencies with a posteriori tests. I don't put the term month because it is not significant. - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context:

[R] duncan test (kruskal post-hoc)

Hi, i´ve read all that i could about the post-hoc test for kruskall-Wallis, now i know that i need to run Duncan test, but i´m not able to find in any package or script. I s anybody knows where can i find it? Thanks - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Paras

[R] Looking for Filliben (correlation test)

Hi R-users, I need to run Filliben correlation test, but I'm not able to find this function in any package. Perhaps it isn't in a package, but sometimes anybody have program this function or create a routine to estimate it. I found it in google and I hope what I found it's in this page, http://gene

Re: [R] Looking for Filliben (correlation test)

Thanks very much, but...the OP, what's that? Sorry - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/Looking-for-Filliben-correlation-test-t

Re: [R] Looking for Filliben (correlation test)

Thanks very much, I got, it runs perfect. - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/Looking-for-Filliben-correlation-test-tp3636069

[R] GLM different results with the same factors

I've read something about this problem, but I don't know how can i avoid this problem. Why the order of the factors give different results? I suppose it's because the order of the factors, i've just changed "lcc" from the first position to the last in the model, and the significance change complete

Re: [R] GLM different results with the same factors

Ok, thanks, I've been reading these days about what you tell me, but i don't understand properly. How could I know, with this tests, which variables are significant? I know my dependent variable depends on the lcc and on the edadysexo. but only one per test seems to be significant. Thanks again -

[R] R in Linux (Ubuntu)

Hi, I´m trying to move to ubuntu at all, but I don´t like R in Ubuntu because it runs in the terminal, so i can´t access directly to some options R gots in Windows. Is there any R can i install in ubuntu which runs in hos own prompt or terminal?? Thanks - Mario Garrido Escudero PhD student Dpt

Re: [R] GLM different results with the same factors

Thanks very much again, I´m reading some papers and articles about this issue and I think i´m starting to understand the problem. And thanks for the link to Professor Fox about the non-sequential Anova. I'll be back with more doubts. I'm sure of that. - Mario Garrido Escudero PhD student Dpto

Re: [R] R in Linux (Ubuntu)

Thanks very much, and...Which do you recommend? - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/R-in-Linux-Ubuntu-tp3703329p3704107.html

Re: [R] R in Linux (Ubuntu)

Thanks very much everybody. I'm going to try all the options and tell you - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/R-in-Linux-Ubunt

[R] Changing labels positions in two graphics

Hi, I´ve got some problems with the labels of the x-axis I,ve got two factors with two categories each: sex (males and females) and area (central and peninsulae), but because of the lenngth in the graphic just appeared two of the four labels. I thought the solution could be rotate it 45º degrees

Re: [R] Changing labels positions in two graphics

Sorry for the late answer. You both put me iin the correct way. Problem is solved - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/Changi

[R] graph of kruskalmc

Hi, I just get interesting results owith non-parametric kruskalmc, and i want to get a graph similar to those obtained for tukey. I know how to do it for tukey ("plot(TukeyHSD(model))"). Anibody knows how to get a graph for kruskalmc Thanks - Mario Garrido Escudero PhD student Dpto. de Biolo

[R] different results from fligner.test

I've made fligner test with the same data, changing the orders of the variables, and this what i get > fligner.test(rojos~edadysexo*zona*ano*estacion) Fligner-Killeen test of homogeneity of variances data: rojos by edadysexo by zona by ano by estacion Fligner-Killeen:med chi-squared =

Re: [R] different results from fligner.test

Yeah! it works, i know it because of the degrees of freedom. This is what i get now > g <- interaction(ano, edadysexo, zona, estacion) > fligner.test(rojos ~ g) Fligner-Killeen test of homogeneity of variances data: rojos by g Fligner-Killeen:med chi-squared = 249.7591, df = 87, p

[R] problems with glht for ancova

I've run an ancova, edadysexo is a factor with 3 levels,and log(lcc) is the covariate (continous variable) I get this results > ancova<-aov(log(peso)~edadysexo*log(lcc)) > summary(ancova) Df Sum Sq Mean Sq F value Pr(>F) edadysexo2 31.859 15.9294 803.9843 <2e-

Re: [R] problems with glht for ancova

I answer myself, perhaps it can be useful to anyone. The problem is that I tried to run a psot hoc test before i get the minimal model by deleting the interaction, if I do it: > ancova2<-aov(log(peso)~edadysexo+log(lcc)) > anova(ancova,ancova2) Res.DfRSS Df Sum of Sq F Pr(>F) 150

Re: [R] kruskal wallis post hoc?

Try to use "kruskalmc" in the package pgirmess I give one of my results > kruskal.test(rojos~mes) Kruskal-Wallis chi-squared = 132.3091, df = 3, p-value < 2.2e-16 > kruskalmc(rojos~mes) Multiple comparison test after Kruskal-Wallis p.value: 0.05 Comparisons obs.dif critical.dif

Re: [R] problems with glht for ancova

Thanks very much - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/problems-with-glht-for-ancova-tp4284908p4288227.html Sent from the R help

[R] glht (multicomparisons) with an interaction factor

Hi, i was working with this model > mq<-glm(rojos~edadysexo*zona*estacion,quasipoisson) and i get this minimal adequate model > anova(mq5,test="F") Df Deviance Resid. Df Resid. DevFPr(>F) NULL518 64799 edadysex

[R] GLHT in multcomp: Two similar models, one doesn't work

i ran this model > model2<-glm(rojos~ageandsex+sector+season+sector:season,quasipoisson) > glht(model2,linfct=mcp(ageandsex="Tukey")) General Linear Hypotheses Multiple Comparisons of Means: Tukey Contrasts Linear Hypotheses: Estimate M - H == 0 0.2898 SUB - H == 0

Re: [R] GLHT in multcomp: Two similar models, one doesn't work

Nobody? - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.789695.n4.nabble.com/GLHT-in-multcomp-Two-similar-models-one-doesn-t-work-tp4291875p4292889.html Sent f

[R] glht (multicomparisons) with a binomial response variable

Hi, I,ve a run a model like this mcrm<-glm(catroj~month,binomial) being catroj a binary response variable with two levels (infected and non infected) > anova(mcrm3,test="Chisq") Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 520 149.81 mes 3

Re: [R] glht (multicomparisons) with a binomial response variable

Thanks very much, Frank, I will try, first I have to adjust your solution to my problem cause i have just for months : april, may, june and july. I will tell you!! Thanks again - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola

[R] pairwise comparisons with multcomp package

Hi, I've got this model and following Hothorn et al advices, I used glht for a post hoc comparison > modezqM<-glm(rojos~estacion*zona3,quasipoisson,subset=(edadysexo=="M")) > anova(modezqM,test="F") Df Deviance Resid. Df Resid. Dev F

[R] problems with lme function

I try to fit a repeated measures model, in which: FA: The difference in the number of scales between the left and right sides for a given trait trait: the each of the measured features are 6 for each individual Condition: is the state of condition of each individual Individual: Each of the specimen

[R] R crash when i'm using lme function

When I try to adjust a mixed model with random effects I can make this order without problem > lm.FA<-lme(absFA~trait*condition,random=~1|individual) But if I try to fit a model in which the response (absFA) is not the same in all individuals at different levels of "trait" factor , but varies ran

[R] Rank a numerical variable

hello, Is there any function in R to transform a numerical continuos variable in a ranked variable? Thanks - Mario Garrido Escudero PhD student Dpto. de Biología Animal, Ecología, Parasitología, Edafología y Qca. Agrícola Universidad de Salamanca -- View this message in context: http://r.78