search for: mandeville

...rrelation Analysis for the Behavioral Sciences, Third Edition) on page 273 state the covariance matrix of the regression coefficients is provided by standard programs for multiple regression, including SAS, SPSS, and SYSTAT. How can it be calculated with R. Thank you very much. pbm Peter B. Mandeville cel: 444 860 3204 tel: 52 444 826 2346-49 ext 532 fax: 52 444 826 2350 P.D. Favor de confirmar la llegada de este correo. Gracias. [[alternative HTML version deleted]]

# Peter B. Mandeville kindly offered to send me code for Hotelling's T^2 # Test. Unfortunately there seems to be no route to his machine. # So i'm trying to reach him via the Mailing List. ------------------------------------------------------------------------ Sir, this morning i recieved your message about...

I respost the message here since I suppose P. Mandeville is on the list, question is rather general and ----- Forwarded message from MAILER-DAEMON at sirio.stat.unipd.it ----- <mandevip at 148.224.17.2>: Sorry, I couldn't find any host named 148.224.17.2. (#5.1.2) --- Below this line is a copy of the message. On Wed, Dec 16, 1998 at 05:01:02A...

...n values to a vector (pvalue). Similar code works in C but not in R. What am I doing wrong? r <- pvalue <- 0 for(i in (1:(k-1))){ for(j in (i+1):k){ r <- r+1 tstat <- (means[i]-means[j])/rms pvalue[r] <- 2*(1-pt(abs(tstat),df)) } } Thank you. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

...tor Pes with 600 elements some of which are NA's. How can I form a vector of the indices of the NA's. for(i in 1:600) if(is.na(Pes[i])) print(i) prints the indices of the NA's but I can't figure out how to put the results in a vector. Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

...the same results. However, I can read the graphic with gimp and if I save the file as an encapsulated postscript file from gimp then the above procedure works. The R document is 4,837 bytes in length and the gimp document is 5,289,399 bytes in length. What am I doing wrong? Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

...t <- (matrix(means, k,k, byrow=TRUE) - matrix(means, k,k,byrow=FALSE))/rms pvalue <- 2*(1-pt(abs(tstat), df=df)) pvalue[upper.tri(pvalue)] S is designed to allow you to work with "whole objects", which means that for loops can (and should) be avoided. Martyn On 05-Oct-99 Peter B. Mandeville wrote: > I am trying to assign values to a vector (pvalue). Similar code works in C > but not in R. What am I doing wrong? > r <- pvalue <- 0 > for(i in (1:(k-1))){ > for(j in (i+1):k){ > r <- r+1 >...

...e an example, Venables and Ripley 1997:381 doesn't either. It isn't mentioned in the scripts nor in the section R Complements nor in Rnotes. I tried aov(c(y1,y2,y3)~x) which seemed reasonable to me but doesn't work. How does it function? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

I am updating my lab manual in Spanish. How should I cite your supplements to Venables and Ripley 1997 and the scripts. Thank you very much Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

Is there any place where the Chambers and Hastie 1993 Statistical Models in S data sets can be found? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

How can I add a least squares line to each panel in coplot? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

.... 7. While running Rw0642 a single click on another minimized application on the task bar highlights the icon for Rw0642 and doesn't change applications. It requieres 2 clicks to change to the other application. Can someone help me with these details? Thank you very much. PBM -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

...d. The command > inflm <- influence.measures(result) reports Warning: NaNs produced in function "qf" and > which(apply(inflm$is.inf,1,any)) reports Error: names attribute must be the same length as the vector What should I do? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

Kjetil Halvorsen mentioned a book call "S-Poetry in R". I checked with amazon.com and Barnes and Noble but they don't have it listed. Can somebody please give me the complete reference and where it can be obtained. Suggestion: Might it be possible to put a list of the most important and new literature about R, S, and S+ which has a direct bearing on R in the Documentation section

...edition. Statistics and Computing. Springer-Verlag New York, Inc., New York, NY, USA. K. A. Weigel, B. A. Craig, T. R. Bidwell and D. M. Bates. 1992. Comparison of Alternative Diphasic Lactation Curve Models under Bovine Somatotropin Administration. J. Dairy Sci. 75:580-589. Peter B. Mandeville cel: 444 860 3204 tel: 52 444 826 2346-49 ext 532 fax: 52 444 826 2352 P.D. Favor de confirmar la llegada de este correo. Gracias. _________________________________________________________________ Discover the new Windows Vista [[alternative HTML version deleted]]

...in qtukey(p, nranges, nmeans, df, lower.tail, log.p) : Non-numeric argument to mathematical function A simple but inelegant workaround is to input mse and df.residual as function arguments in addition to the aov object. Does anyone have a better solution? Thank you, Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

....2330 0.2330 A-C 0.0033 0.0033 0.0168 B-C 0.0850 0.7553 0.0850 Again, why does the contrast contr.sum have distinct results from the other contrasts? Which ones are confiable? Thank you very much Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...

Happy New Year to all, I am using RW0631 and WIN98. 1. The command DATA only works for the base library. 2. The following code on page 312 of Venables and Ripley 1997 sitka.lme <- lme(size~treat*ordered(Time),random=~1,cluster=~tree,data=Sitka, serial.structure="ar1.continuous",serial.covariate=~Time) produced the following error message Error in lme(size ~ treat *

Hi, Does anyone know if the D'agostino test is available with R ? Alex

...that summary(Orthodont.fit2) produces the same fixed effects values (Intercept) 16.340625 Age 0.784375 Sex 1.032102 Age.Sex -0.0304830 as Pinheiro and Bates 1995 (lme and nlme) but everything else is very distinct. Can someone please help? Thank you very much. Peter B. -- Peter B. Mandeville mandevip at deimos.tc.uaslp.mx Jefe del Depto. de Inform?tica y Bioestad?stica rpe1531 at pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Aut?noma de San Luis Potos? Fax: 48 28-23-52 Av. V. Carranza 240...