similar to: Package candisc

Hello, I am having trouble loading the package candisc onto my R distribution. I am using 2.7.1-2. I do a "> install.packages("candisc" and get the following output. Warning in install.packages("candisc") : argument 'lib' is missing: using '/usr/local/lib/R/site-library' --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk

Hi, I am doing canonical discriminant analysis using candisc function from the candisc package. My input is a table of species distribution (columns = abundance of each species in each sample) in samples that are split by categories (rows), and I want to know whether each category is associated with a particular set of species and their abundances. I have 20 rows (samples) split into 6

Hello, I have a question regarding the candisc package. My data are: species three five 1 2.95 6.63 1 2.53 7.79 1 3.57 5.65 1 3.16 5.47 2 2.58 4.46 2 2.16 6.22 2 3.27 3.52 I put these in a table and then a linear model >newdata <- lm(cbind(three, five) ~ species, data=rawdata) and then do a candisc on them >candata<-candisc(newdata)

Hello, I have a file with two dependent variables (three and five) and one independent variable. I do i.mod <- lm(cbind(three, five) ~ species, data=i.txt) and get the following output: Coefficients: three five (Intercept) 9.949 9.586 species -1.166 -1.156 I do a" i.can<-candisc(i.mod,data=i): and get the following output: Canonical Discriminant Analysis

Hi listers, I have the following code... data<-as.matrix(c(1,2,3,4,5,6,7,8,9,10),nrow=10,ncol=1,byrow=TRUE) n<-nrow(data) m<-n-1 boot<-data[sample(m,replace=T),] So, I need to count the number of times each observation was selected at the sample with replacement... Suppose I sampled... 4 5 1 3 7 5 1 4 7 So, I would count... x count 1 2 2 0 3 1 4 2 5 2 6 0 7 2 8 0 9 0 10 0

Hi listers, I have a simple doubt... I need to manipulate the data base as: var1 var2 7 0.1 7 0.1 8 0.12 10 0.15 12 0.18 20 0.31 I would like to aggregate the equal cases as var1 making the sum of var2... I would get... var1 var2 7 0.2 8 0.12 10 0.15 12 0.18 20 0.31 Thanks in advance, Marcio -- View this message in context:

Hi listers, I have a data frame with 5 variables... And I want reorder the data frame according a specific variable X... I checked the sort and order functions... In my case I think that the ORDER function is applicable... But I did not understand at all... The idea is simple... Data order by X... I know that this is very simple... Anyways... How I would do that... Thanks in advance, Marcio --

I've just released the latest R-Forge versions of heplots 0.9-3 and candisc 0.5-18 to CRAN. They should appear there within a day or two. == heplots The heplots package provides functions for visualizing hypothesis tests in multivariate linear models (MANOVA, multivariate multiple regression, MANCOVA, etc.). They represent sums-of-squares-and-products matrices for linear hypotheses and for

I'm happy to announce the candisc package, v 0.5-9, now on CRAN. Generalized Canonical Discriminant Analysis Description This package includes functions for computing and visualizing generalized canonical discriminant analyses for a multivariate linear model. They are designed to provide low-rank visualizations of terms in a mlm via the plot method and the heplots package. The methods

I'm happy to announce the candisc package, v 0.5-9, now on CRAN. Generalized Canonical Discriminant Analysis Description This package includes functions for computing and visualizing generalized canonical discriminant analyses for a multivariate linear model. They are designed to provide low-rank visualizations of terms in a mlm via the plot method and the heplots package. The methods

I'm trying to write a formula method for canonical correlation analysis, that could be called similarly to lm() for a multivariate response: cancor(cbind(y1,y2,y3) ~ x1+x2+x3+x4, data=, ...) or perhaps more naturally, cancor(cbind(y1,y2,y3) ~ cbind(x1,x2,x3,x4), data=, ...) I've adapted the code from lm() to my case, but in this situation, it doesn't make sense to include an

Hi listers, I am working on a simulation... But I am having some troubles... Suppose I have a function A which produces two results (mean and variance)... Then I would like to simulate this function A with a function B many times using the results from function A For example: #Function A boot<-function(a,b,c){ mean_boot<-(a+b)/2 var_boot<-c #list(a=a,b=b,c=c) return(a) } Then I would

Hi listers, I am trying to arrange my data and I didn't find any information how to do it! I have a data with 3 variables: X Y Z 1-I would like to multiplicate de information of X according to the number I have for my Y variable... 2-Then I want to identify with a dicotomic variable by the number according my variable Z from X... I can do the first part by... z<-rep(x,y) But I don't

Hi Listers, Is there a way that I can transpose an array... Suppose I have the following array... x<-array(c(1,2,3,4),dim=c(1,2,2)) , , 1 [,1] [,2] [1,] 1 2 , , 2 [,1] [,2] [1,] 3 4 And I would like to get the following result... , , 1 [,1] [1,] 1 [,2] 2 , , 2 [,1] [,2] [1,] 3 [,2] 4 Thanks in advance, Marcio -- View this message in context:

[Using R 2.7.2 on Windows XP] After re-building our heplots package, I've begun to get the following error from sessionInfo(), even though it passes R CMD check and builds without errors: > sessionInfo() Error in x$Priority : $ operator is invalid for atomic vectors In addition: Warning message: In FUN(c("MASS", "heplots", "car", "rgl",

Hi Listers, I am producing some graphics that the commands are in a FUNCTION... The problem is that I end up viewing just last graphic and in my FUNCTION there are 4 graphics with the PAR command function... Like those below... How do I view/get the other 3 graphics? Any help? Thanks in advance... histogram<-par(mfrow=c(1,2)) hist(rw_mean_app,main='Bootstrap Method

I'm working on some functions for generalized canonical discriminant analysis in conjunction with the heplots package. I've written a candisc.mlm function that takes an mlm object and computes a candisc object containing canonical scores, coeficients, etc. But I'm stumped on how to construct a mlm for the canonical scores, in a function using the *same* right-hand-side of the model

Hi listers, I've been looking for a procedure, but I am not succeding... I have a function that give multiple results... Then, I would like to simulate this function n times, so I need to save/keep the n multiple results, in order to calculate my desired statistics... I have already tried with the RETURN and LIST FUNCTION, but I am not getting it right... An example of what I am looking for

Hi Sergio Doing those tests 30 times is going to give you a huge Type I error rate, even if there was a function that did that. There is a reason why TukeyHSD doesn't make it easy. In general, if there are useful comparisons among the species, you are better off setting up and testing contrasts than doing all-pairwise Tukey tests. Also, the PCA scores are ordered in terms of variance

Dear R list, I'm trying to calculate Mahalanobis distances for 'Species' of 'iris' data as obtained below: Squared Distance to Species From Species: Setosa Versicolor Virginica Setosa 0 89.86419 179.38471 Versicolor 89.86419 0 17.20107 Virginica 179.38471 17.20107 0 This distances above were obtained with proc