search for: mass3

I was unable to open this file Bill Venables' excellent "Exegeses on Linear Models" posted at http://www.stats.ox.ac.uk/pub/MASS3/Exegeses.ps.gz I'd be very interested in reading it? Thanks Esther Meenken Biometrician Crop & Food Research Private Bag 4704 Christchurch TEL: (03) 325 9639 FAX: (03) 325 2074 EMAIL:MeenkenE at crop.cri.nz Visit our website at http://www.crop.cri.nz ___________________________...

On p. 48 of "Statistics Complements" to the 3rd MASS edition, http://www.stats.ox.ac.uk/pub/MASS3/VR3stat.pdf I read that the orthogonal rotations of Z Lambda^-1 remain uncorrelated, where Z is the PC and Lambda is the diag matrix of singular values. However, the example below that text is > A <- loadings(ir.pca) %*% diag(ir.pca$sdev) If ir.pca$sdev are the singular values, should...

...message appears: "Error in .Fortran ("vectnc",: "VECTNC" is not a symbol in the load table". May be I did not installed it correctly or (as I suppose) it is incompatible with this version of Splus. I looked in the directory for Splus 6.0 http://www.stats.ox.ac.uk/pub/MASS3/Winlibs but the update of this function is not there. Do you know of some alternative function for the Non-central T in Splus 4.5 or how to solve the problem with noncent.zip? b) I wanted to apply the overall test for coincidental regressions (see Zar, pag. 304), whose statistic is: F = ((SSt-SSp...

I am new to R. I am going to by one of the following book: 1. William N. Venables and Brian D. Ripley. Modern Applied Statistics with S-Plus. Third Edition. Springer, 1999. ISBN 0-387-98825-4. 2. The Fourth Edition of the book from point 1. 3. `S Programming' by W. N. Venables and B. D. Ripley Springer. ISBN 0-387-98966-8, 2000. I can only by one of the above books. Q1. I have found

...s of squares for an ANCOVA for my specific experimental data and hypotheses. I do have a basic understanding of the differences between Type-I, II, and III SSs, have read about the principle of marginality, and read Venable's "Exegeses on Linear Models" (http://www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf). I am pretty new to R and a search of the R-help archive did not answer my question (although I found some good pointers). In brief, leaving my covariates aside, I hypothesize that women (a) generally perform lower then men in a specific task (microworld performance, MWP) and that t...

...andling, and sm uses chull. This will be the last version of the VR5.x software for the second edition, as a new more `R friendly' third edition is due out in July. The software for that is now available for pretest (including for R-0.64 on UNix and Windows) at http://www.stats.ox.ac.uk/pub/MASS3 It is almost completely backwards compatible (just no support for survival4) with eariler versions but has new datasets and functions. Some, such as ppr and lqs, are already in R. Windows versions: ================= I have split the Windows contributed collection into current and rw063x section...

...andling, and sm uses chull. This will be the last version of the VR5.x software for the second edition, as a new more `R friendly' third edition is due out in July. The software for that is now available for pretest (including for R-0.64 on UNix and Windows) at http://www.stats.ox.ac.uk/pub/MASS3 It is almost completely backwards compatible (just no support for survival4) with eariler versions but has new datasets and functions. Some, such as ppr and lqs, are already in R. Windows versions: ================= I have split the Windows contributed collection into current and rw063x section...

Hi to all, how can I perform a multiple correspondence analysys with a block contingence matrix? What's the proper package that must I load into R? Andrea Rossetti. PS Excuse me if my question is stupid and for my bad english!!! -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

...le. Our sincere thanks to all of you who, directly or indirectly, have helped us in this second project. (Some of you may even recognise answers to your questions.) We hope all of you who use it will find it clear and informative. Further details are available at http://www.stats.ox.ac.uk/pub/MASS3/Sprog/ and at the mirror sites for our software. Brian Ripley and Bill Venables -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road,...

Hello, I want to print points having value on 3 variables in a triangle (ternary diagram). How can I di this with R ? -- Cordialement ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~ Emmanuel POIZOT ~ CNAM/INTECHMER ~ B.P. 324 ~ 50103 CHERBOURG CEDEX ~ T?l : (33) 233 887 342 ~ Fax : (33) 233 887 339 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Prof Brian Ripley wrote: > > Even better, use the function fitdistr in package MASS... > I had a look in my MASS library (from the package VR_6.2-6) and couldn't find this function. Is there a newer version available? Thanks for any help. Jim This email message and any accompanying attachments may contain confidential information. If you are not the intended recipient, do not

Hi, Is it any way to convert .Rd files to CHTML Splus 6 files under Windows? Thanks, Philippe Grosjean ...........]<(({?<...............<?}))><............................... ( ( ( ( ( ) ) ) ) ) Philippe Grosjean ( ( ( ( ( ) ) ) ) ) IFREMER Nantes - DEL/AO ( ( ( ( ( rue de l'Ile d'Yeu, BP 21105, 44311 Nantes Cedex 3 ) ) ) ) ) tel: (33) 02.40.37.42.29, fax: (33)

Dear list! I am interested in learning about MLE and I wonder whether it is possible to use the examples for maximum likelihood estimation given in 8.7 in MASSv3 with R? AFAIU R does not have a direct replacement for S-PLUS's ms() which the examples use for the fitting, but optim() may be of help for me. However, I am not sure how I can convert the use of ms() using optim()--is there an

Dear useRs, I need to compute zero of polynomial function fitted by lm. For example if I fit cubic equation by fit=lm(y~x+I(x^2)+i(x^3)) I can do it simply by polyroot(fit$coefficients). But, if I fit polynomial of higher order and optimize it by stepAIC, I get of course some coefficients removed. Then, if i have model y ~ I(x^2) + I(x^4) i cannot call polyroot in such way, because there is

Dear helpeRs, I'm estimating a series of linear models (using lm) in which in every new model variables are added. I want to test to what degree the new variables can explain the effects of the variables already present in the models. In order to do that, I simply observe wether these effects decrease in strength and / or lose their significance. My question is: does any of you know

Dear R Gurus: I vaguely remember reading that if interaction was present in a factorial design, then the main effect results were suspect. However, I was reading a text which now uses the tests for main effects even if interaction is present. Which is correct, please? Thanks, Edna Bell

Dear all,   Does some body have idea on how to extract Type III sum of Square from "lm" or "aov" function in R ? I could not figure out it.   If this is minor and irrelevant to post in this mail, I am sorry for that. Thanks.     Sincerely, Ram Kumar Basnet Wageningen University, Netherlands   [[alternative HTML version deleted]]

...interaction or the quadratic terms: drop1(m) ... Df Sum of Sq RSS AIC <none> 0.000033 -254.306 I(a^2) 1 0.000000 0.000033 -256.306 I(b^2) 1 0.098611 0.098644 -96.239 a:b 1 0.000032 0.000065 -242.674 However, this: http://www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf suggests that marginality rules between powers of variables might not be implemented (although they might have been since 2000). My question is: I am "allowed", according to marginality rules, to remove a^2? I have found plenty of information on how the coefficients corresp...

Hi I have some data (v. large amount) with a (0,1) response where I want to minimise the errors in the betas rather than SS or deviance. So can anyone point me to a ridge regression function or equivalent for such a logistic regression case? John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello all, I've become confused by the output produced by a call to anova(model1,model2). First a brief background. My model used to predict final tree height is summarised here: Df Sum Sq Mean Sq F value Pr(>F) Treatment 2 748.35 374.17 21.3096 7.123e-06 *** HeightInitial 1 0.31 0.31 0.0178 0.89519