[R] Non-linearity with Parametric data

2008-04-14 Thread Charles Willis
Hello, I am trying to test for non-linearity in a set of non-parametric data. Furthermore, I would like to do so in a multi-variate model. Frankly, I don't even know if this is possible, and if it is possible, I don't know if it is implemented in an R package. Any advice would greatly appreciated.

[R] Writing GLM/GEE Ouput to Text File

2008-03-17 Thread Charles Willis
Hello, How can one write the output of a GLM/GEE model to a text file, such that the results appear in text file in the same (or similar) format they appear in the R console? For instance, I have the following model: Traitresult <- compar.gee(Y~X, data = data, family = "binomial", phy = tree, sc

[R] Convert a List of Distances to a Distance Matrix

2008-03-12 Thread Charles Willis
Hello, Is there an easy function for switching list to matrix. My list is of genetic distances between species pairs: A A 0 A B .5 A C .25 B C .5 and I want a distance matrix such as: A B C A 0 .5 .25 B .5 0 .5 C .25 .5 0 for use in a mantel test. Thank you for the help! cheers, charl

[R] Matrix Row name Compatibility

2008-03-09 Thread Charles Willis
Hello, I am trying to preform a comparative analysis using the partial mantel method. However, i am having trouble creating compatible matrices. I can create matrices from my numerical data using the 'distance(x, method)' function from ecodist, but the row/column names are out put as numbers. Whil

[R] COMPAR.GEE Output

2008-03-01 Thread Charles Willis
Hello, I am running the program COMPAR.GEE within the package APE. My dependent variable is binomial, while my independent variable is a multi-state categorical variable. The output reports an estimate for each state of the independent variable except the first one. For example, for the variable X

[R] COMPAR.GEE error with logistic model

2007-12-29 Thread Charles Willis
Hello, I am trying to run the APE program COMPAR.GEE with a model containing a categorical response variable and a mixture of continuous and categorical independent variables. The model runs when I have categorical (binary) response and two continuous independent variables (VAR1 and VAR2), but whe