similar to: Canonical analysis with class variable?

Dear list, I'm trying to compare two sets of variables using canonical analysis. My X variables include 3 climate indices, all continuous. My Y variables describe a set of 3 environmental moisture measurements. However, I would also like to include a class variable for habitat type. My colleaugue was able to do this in SAS, but we could not find a way to do it in R. Is anyone familiar

Bonjour, we wanted to fit a manova as descripted in Marieta /et al./ 2003, "convergent habitat segregation of /Aedes aegyptii/ and /Aedes albopictus/ (Diptera: /Culicidae/) in Southern Brazil and Florida", /J. Med. Entomol./, *40* (6), 785-794. They did their analysis with SAS software. We have the same kind of data with two insect species and city, habitat and season factors. We

Hello, I have been attempting to run a poisson glmm using lme4 for some time now and have had a lot of trouble. I would say 9 times out of 10 I receive the following warning: CHOLMOD warning: %h Error in mer_finalize(ans) : Cholmod error `not positive definite' at file:../Cholesky/t_cholmod_rowfac.c, line 432 My data are counts of microbe colony forming units (CFUs) collected from

Hello everybody ! I try to obtain in R eigenvectors and canonical analysis on MANOVA results, but I don't find how to process? In particular, I would be interesting to obtain "standardized canonical coefficients" of the canonical variates. There analysis give some information on the correlation between response variates. My data are organised in 2 terms (one is continu, one is a

Hey, all, I'm trying to construct a centroid plot using canonical values from a MANOVA. I know that from the summary.manova object you can get Eigenvalues, and the H and E matrices (from SS$Treatment and SS$Residuals), but I am at a loss to get the loadings for the canonical values, nor values for the centroid centers and radii. Is there a package that does this that I am just missing,

I think i mentioned this in private communication a while back, but in case anyone else cares: I have a version of canonical correspondence analysis (CCA) coded up in R. I based it on Ter Braak and Prentice (1988) Adv. Ecol. Res. 18:271-317, Ter Braak (1986) Ecol. 67(5):1167-1179, and Ter Braak (1995) Section 5.9 in Jongman, Ter Braak, and van Tongeren "Data Analysis in

I am working with a RCB 2x2x3 ANCOVA, and I have noticed a difference in the calculation of sum of squares in a Type III calculation. Anova output is a follows: > Anova(aov(MSOIL~Forest+Burn*Thin*Moisture+ROCK,data=env3l),type=3) Anova Table (Type III tests) Response: MSOIL Sum Sq Df F value Pr(>F) (Intercept) 22.3682 1 53.2141 3.499e-07 *** Forest

Dear list, After searching many old posts, I can't find the solution to a simple problem. can someone tell me how to create a character string with multiple backslashes, as in: file_dir <- c("C:\files\data\") I need to create this string and then paste it to many files names for batch processing in another software program. R won't accept the backslash and removes

Dear List, I'm stuck on what seems like a simple indexing problem, I'd be very grateful to anyone willing to help me out. I queried a dataframe which returns a character vector called "plot". I have another dataframe from which I want to subset or select only those rows that match "plot". I've tried subset, and also the "which" command. plot

Dear list, I'm trying to write my first looping function in R. After many hours of searching help files and previous posts, I'm at wits end. Please forgive my programming ignorance...any help is greatly appreciated. I need to sort through a vector (x) and identify the point at which 2 successive values become smaller than the previous value. I've written a "while"

Dear all, In a greenhouse experiment we tested performance of 4 different species (B,H,P,R) under 3 different water levels in 10 replications. As response variable e.g. the number of emerging sprouts were measured on three dates. A simple Anova considering every measurement date separately shows a higly significant effect of species and moisture (and partly the interaction of both). The

Dear list, I'm trying to download a text file directly from the internet using the RCurl package and the command getURL. Duncan Lang graciously helped me solve the first step in this problem using the following command: ################# txtfile <- getURL('ftp://ftp.wcc.nrcs.usda.gov/data/snow/snow_course/table/history/idaho/13e19.txt', ftp.use.epsv = FALSE) #################

I have a data.frame paint (below-mentioned). In SPlus I used the command raov( MOISTURE ~ BATCH / PROBE ) (you could do raov( MOISTURE ~ BATCH + PROBE%in%BATCH) as well) so that the factors are taken as random. In R this function raov doesn't exist. How can I calculate the same? Maybe with lme, but how? str(paint) ---------- `data.frame': 60 obs. of 5 variables: $ BATCH : Factor

Dear list, I would like to plot 2 series of numbers with very different ranges/scales as lines on the same plot. I assumed this is commonly done and easy, but I have not found any help files (e.g. axis() or matplot() that show how. I've searched many old posts to no avail. I'll be very grateful for any suggestions on how this is done. Best, Zack

Dear list, I am trying to download a text file from an ftp site using download.file(). I used the following code: url <- "ftp://ftp.wcc.nrcs.usda.gov/data/snow/snow_course/table/history/idaho/13e19.txt" dest <- "c:/test/downloadtest.txt" download.file(url, dest) I get this error message, indicating that R is unable to open the connection: trying URL

I am trying to it a particular nonlinear model common in Soil Science to moisture release data from soil. I have written the function as shown below according to the logist example in Ch8 of Pinheiro & Bates. I am getting the following error (R version 2.1.1) *Error in qr(attr(rhs, "gradient")) : NA/NaN/Inf in foreign function call (arg 1)* Below is the function and data. /#

Hello all R listers, I'm struggling to select an appropriate statistical method for my data set. I have collected soil moisture measurements every hour for 2 years. There are 75 sensors taking these automated measurements, spread evenly across 4 treatments and a control. I'm not interested in being able to predict soil future soil moisture trends, but rather in knowing whether the

Dear list, I've recently installed R on a 64 bit machine with 8 GB of RAM. I set this computer up as a dual-boot system, with windows XP 64 and Ubuntu 7.10. I downloaded the Linux 64 bit version of R and installed it. I'm trying to run rather large Random forest models and was running into severe memory limitations with my old windows 32 bit machines. When I run Random Forest models

Greetings R list, I am stuck on a simple syntax problem. I want to list all files in a directory, excluding files of a certain type. I have tried pattern matching as follows: a <- list.files(data, full.name = TRUE, pattern != ".xml") # exclude all .xml files The warning returns that my syntax is incorrect. I have read the regexpr help files and search old posts to no

Dear list, I'd like to use path.analysis in the package agricolae in batch format on many files, retrieving the path coefficients for each run and appending them to a table. I don't see any posts in the help files about this package or the path.analysis package. I've tried creating an object out of the call to path.analysis, but no matter what I try, the function automatically prints