Displaying 20 results from an estimated 5000 matches similar to: "Still needing help with LDA"
2002 Mar 30
0
Wilks' Lambda for lda?!
Hello!
As I remember there were questions concerning wilks' lambda lately.
But I lost those e-mails. So I dare ask again.
Is there a way of calculating wilks' lambda for lda?
Christoph.
--
Christoph Lange clange at epost.de
Verhaltensbiologie, FU Berlin 838-55068
Haderslebener Str. 9, 12163 Berlin
2002 Apr 14
1
Suggestion for implementation
I wonder if it was possible (and desireable) to implement means to
protect data objects (and functions) from overwriting them. So for
instance:
> x <- 5
> protect(x)
> x <- 8
Error: x is read only
Or would that be against the philosophy of R to be compatible to S?!
Just a RFC ;-)
Yours,
Christoph.
ps. Your probably know the situations that let you think of such a
feature
2007 Feb 22
1
MANOVA usage
Hello,
I had a couple questions about manova modeling in R.
I have calculated a manova model, and generated a summary.manova output
using both the Wilks test and Pillai test.
The output is essentially the same, except that the Wilks lambda = 1 -
Pillai. Is this normal? (The output from both is appended below.)
My other question is about the use of MANOVA. If I have one variable which
has a
2007 Mar 16
2
MANOVA permutation testing
Hi,
I've got a dataset with 7 variables for 8 different species. I'd like
to test the null hypothesis of no difference among species for these
variables. MANOVA seems like the appropriate test, but since I'm
unsure of how well the data fit the assumptions of equal
variance/covariance and multivariate normality, I want to use a
permutation test.
I've been through CRAN looking at
2008 Apr 03
3
summary(object, test=c("Roy", "Wilks", "Pillai", ....) AND ellipse(object, center=....)
Dear All,
I would be very appreciative of your help with the following
1). I am running multivariate multiple regression through the manova() function (kindly suggested by Professor Venables) and getting two different answers for test=c("Wilks","Roy","Pillai") and tests=c("Wilks","Roy",'"Pillai") as shown below. In the
2007 Jun 19
2
How to compute Wilk's Lambda
Dear helpeRs,
the following data set comes from Johnson/Wichern: Applied Multivariate
Statistical Analysis, 6th ed, pp. 304-306.
/X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3, 2)))
Y <- structure(c(0, 2, 4, 0), .Dim = as.integer(c(2, 2)))
Z <- structure(c(3, 1, 2, 8, 9, 7), .Dim = as.integer(c(3, 2)))/
I would like to compute Wilk's Lambda in R, which I know is 0.0385.
2008 Aug 10
2
ANOVA help
Hi,
I'm doing anova on a matrix of multivariate data where I want to assess the
effect of each column (element).
My matrix is 86 rows x 31 columns. I've created a grouping factor of length
86 containing group assignments of 6 types.
Then I run:
x<- aov(matrix~grouping.factor)
summary(aov.fit.raw, test="Wilks")
This is working fine enough, but I'm getting different
2008 Aug 13
1
summary.manova rank deficiency error + data
Dear R-users;
Previously I posted a question about the problem of rank deficiency in
summary.manova. As somebody suggested, I'm attaching a small part of
the data set.
#***************************************************
"test" <-
structure(.Data = list(structure(.Data = c(rep(1,3),rep(2,18),rep(3,10)),
levels = c("1", "2", "3"),
class =
2009 Mar 15
1
Bug Report Fwd: MANOVA Data (PR#13595)
Hi.? There appears to be a bug in R function manova.? My friend and I both ran it the same way as shown below (his run) with the shown data set. His results are shown below. we both got the same results.? I was running with R 2.3.1. I'm not sure what version he used.
Thanks very much,
David Booth
Kent State University
-----Original Message-----
From: dvdbooth at cs.com
To: kberk at
2004 Jun 02
1
Manova and contrasts
Hi R-users
I'm trying to do multivariate analysis of variance of a experiment with
3 treatments, 2 variables and 5 replicates.
The procedure adopted in SAS is as follow, but I'm having difficulty in
to implement the contrasts for comparison of all treatments in R.
I have already read manuals and other materials about manova in R, but
nothing about specific contrasts were found in them,
2002 Mar 21
1
Wilks Lamba
Hi all,
I can't figure out how to compute Wilks Lambda in a one way repeated
measure design. My matrix looks like:
> t2.m
Blank ECR ENC UEA UED
1 -0.15 0.14 0.16 0.09 0.14
2 0.30 0.08 0.14 0.14 0.14
[...]
where each row is a case and the columns are levels of one factor (named
trial):
> t2.fit <- manova(t2.m ~ 1)
> summary(t2.fit, intercept=T,
2006 Feb 16
2
MANOVA: how do I read off within and between Sum-of-Squares info from the manova result?
Hi all,
I am experimenting the function "manova" in R.
I tried it on a few data sets, but I did not understand the result:
I used "summary(manova_result)"
and "summary(manova_result, test='Wilks')"
and they gave a bunch of numbers...
But I need the Sum-of-Squares of BETWEEN and WITHIN matrices...
How do I read off from the R's manova results?
Any
2002 Oct 09
1
Large F-value and small P-value
Hi all,
I computed a Wilks Lambda Test with manova:
> df.man <- manova(df.mul ~ 1, na.rm=TRUE)
> summary(df.man, intercep=T, test="Wilks")
Df Wilks approx F num Df den Df Pr(>F)
(Intercept) 1 0.0002824 393.3 9 1 0.03911 *
Residuals 9
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*'
2006 Mar 30
2
Unbalanced Manova
Dear all,
I need to do a Manova but I have an unbalanced design. I have
morphological measurements similar to the iris dataset, but I don't have
the same number of measurements for all species. Does anyone know a
procedure to do Manova with this kind of input in R?
Thank you very much,
Naiara.
--------------------------------------------
Naiara S. Pinto
Ecology, Evolution and Behavior
1
2003 Mar 17
0
Built-In Wilks Lambda for lda?
Hello,
using the lda-method from MASS-package I was wondering whether there is
a built-in method for figuring out Wilks' Lambda?
Searching the Web I found in the r-help archive a thread form june 2002,
but it didn't help me.
I understand I can use manova and its summary-method to get Wilks' Lambda
on the screen, but I don't see the connection to lda from MASS.
Or does
2011 Mar 20
3
manova question
Dear friends,
Sorry for this somewhat generically titled posting but I had a question
with using contrasts in a manova context. So here is my question:
Suppose I am interested in doing inference on \beta in the case of the
model given by:
Y = X %*% \beta + e
where Y is a n x p matrix of observations, X is a n x m design matrix,
\beta is m x p matrix of parameters, and e is a
2018 May 25
2
TukeyHSD for multiple response
Dear all,
I'm testing the effect of species and sex in my sample by using the principal component scores of a PCA analysis.
I have 30 PCs and I tried to see if there is any significant difference from males to females, given that there is a significant effect of phylogeny (factor with several species).
I didi it like this:
Y<-PCA$pc.scores[,1:30]
fit <- manova(Y ~ sp*sex)
2012 Jun 19
1
Stepwise Discriminant Analysis - greedy.wilks
I don't understand which is the problem
Could you help me please?
Thanks in advance
Marta
> str(data_indiciN2)
'data.frame': 200 obs. of 36 variables:
$ gruppo: Factor w/ 2 levels "0","1": 2 2 2 2 2 2 2 2 2 2 ...
$ I001N2: num 19.32 8.22 28.35 7.24 14.7 ...
$ I002N2: num 2.92 2.54 0.11 1.6 7.12 ...
$ I003N2: num -22.362 -0.222 -19.291
2001 Nov 08
2
Manova in R vs. SAS
While I was helping a SAS-using friend with an analysis I noticed some
differences in the multivariate test statistics, approximate F statistics,
and p-values in the manova function using R and proc GLM using SAS. The
univariate coefficients are identical. Is there a reason to expect R and SAS
to give different results?
Thanks,
Bill Kristan.
2010 Jun 15
1
MANOVA proportion of variance explained
Hello everybody
After doing a MANOVA on a bunch of data, I want to be able to make some comment on the amount of variation in the data that is explained by the factor of interest. I want to say this in the following way: XX% of the data is explained by A.
I can acheive something like what I want by doing the following:
X <- structure(c(9, 6, 9, 3, 2, 7), .Dim = as.integer(c(3,