[R] Changing order of x factor in interaction plot

2008-12-05 Thread Kota Hattori
Hi all, I would like to ask one simple question. I wanna know how you can ask R to change the order of x factor when you draw interaction plot. I find that R puts factors in an alphabetic order (e.g., dog, lamb, and monkey etc). But I simply want to control the order. How can I do this? Aloha,

[R] Specifying data point locations on x-axis using boxplot() and points()

2011-08-17 Thread Kota Hattori
Hi all, I would like to ask a question with regard to R graphics. I am currently trying to add data points in boxplots using points(). In my script, I add up boxplots three times using add=TRUE. I use points() every time I add boxplots, and I see data points. However, I see all data points are

[R] Post hoc power analysis for mixed-effects models

2013-06-04 Thread Kota Hattori
Dear all, I have been searching ways to run power analysis for mixed-effects models. However, I have not been successful in the research. Today I would like to ask your help. As long as I see from my search, Martin Julien wrote a package called pamm for the power analysis. One of the

Re: [R] Post hoc power analysis for mixed-effects models

2013-06-06 Thread Kota Hattori
Hi David, Thank you very much for your reply. I do understand your point. I was requested to do power analyses for null findings by reviewers. I think simulations would be an alternative choice given the validity of post hoc power analysis is questionable. The package, pamm does simulation.

[R] Permutation test for GLM with proportional data

2013-06-12 Thread Kota Hattori
Dear all, I am trying to run permutation tests for GLM using the glmperm package. When I created my model, I followed The R book and Mixed effects models and extensions in ecology with R. In both books, response variables were specified by using cbind(). That is, response variable is in a

[R] Kendall package tau-a, b, and c

2012-08-19 Thread Kota Hattori
Hi all, I would like to ask a question related to Kendall package. I ran Kendall (x,y) and saw the results. But I am not sure which tau values R reported. I have ties in my data set, so I want tau-b. Can anybody tell how Kendall package is calculating tau values? I have looked at the