[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.-

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

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:

[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º

[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:

[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

[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

[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

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] 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

[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

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) 1509

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

[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 edadysexo

[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-16

[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,

[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

[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,

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] 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 2

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:

[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

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

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] 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] 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

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:

[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,

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:

[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

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:

[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 is

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

[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] 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

[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

[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

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