Re: [R] Rcmdr scatter3d

, reproducible code. John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Row-Echelon Form

]))/D[j] } k - 1:(n - 1) D[n] - X[n, n] - sum((L[n, k]^2) * D[k]) if (abs(D[n]) tol) stop(matrix is numerically singular) L %*% diag(sqrt(D)) } John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S

Re: [R] perception of graphical data

://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox

Re: [R] An issue with White's test in Anova

John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED

Re: [R] Rcmdr window border lost

else as well? Thanks in advance, Andy John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox __ R-help@stat.math.ethz.ch

Re: [R] SEM for categorical data

on this. Thank you Upasna -- - Upasna Sharma Research Scholar Shailesh J. Mehta School of Management, Indian Institute of Technology, Bombay Powai, Mumbai - 400076, India John Fox, Professor Department

Re: [R] termplot with uniform y-limits

the comparative contributions of the different components. Notes: the current version of termplot does not allow the user to specify ylim. I checked. the plot tools that come with mgcv do this by default. Thanks Simon. John Fox, Professor Department

Re: [R] SEM model fit

= names(dx)[4:6], par.size = 's', maxiter = 1e3, analytic = F, gradtol = 1e-10) ##set raw to False summary(obj = sem.enf.rq, dig = 3, conf = 0.9) Respectfully, Frank Lawrence John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario

Re: [R] Inverse BoxCox transformation

Dear Des, The following should do the trick: invBoxCox - function(x, lambda) if (lambda == 0) exp(x) else (lambda*x + 1)^(1/lambda) I hope this helps, John John Fox, Professor Department of Sociology McMaster University Hamilton

Re: [R] linear hypothesis test in gls model

this helps, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto

Re: [R] Refactor all factors in a data frame

Dear Hilmar, You could use something like DF - as.data.frame(lapply(DF, function (x) if (is.factor(x)) factor(x) else x)) Where DF is the data frame. I hope this helps, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario

[R] [R-pkgs] Rcmdr 1.3-0 and RcmdrPlugins.TeachingDemos

be appreciated. I'd like to acknowledge Richard Heiberger's suggestions for the design of this plug-in facility. John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] GUI component Margin on tkcanvas, tkframe or tktoplevel

Dear Hao, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Hao Liu Sent: Thursday, May 24, 2007 3:53 PM To: r-help@stat.math.ethz.ch Subject: [R] GUI component Margin on tkcanvas, tkframe or tktoplevel Dear gurus: I have a question on how to

Re: [R] Selecting complementary colours

Dear Chuck, This solution works reasonably well for me. Although it occasionally produces an error, I'm able to trap that. Thank you -- and to everyone else who responded. John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario

Re: [R] Boostrap p-value in regression [indirectly related to R]

, I believe, and differs trivially from the one in Efron and Tibshirani. Regards, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] Selecting complementary colours

Dear Thomas, This seems simpler than the solution that I used, so I'll give it a try. Thanks, John On Tue, 22 May 2007 09:01:01 -0700 (PDT) Thomas Lumley [EMAIL PROTECTED] wrote: On Mon, 21 May 2007, John Fox wrote: In retrospect, I didn't specify the problem clearly: What I want

[R] Selecting complementary colours

Dear r-helpers, I wonder whether, given the #rrggbb representation of a colour, there is a simple way to select the complementary colour, also expressed as a #rrggbb string. Any suggestions would be appreciated. John John Fox, Professor Department of Sociology

Re: [R] Selecting complementary colours

that this is better described as a contrasting than a complementary colour. Your solution, for example breaks down for grays: mycol - #88 do.call(rgb, as.list(1 - col2rgb(mycol) / 255)) [1] #77 Thank you for the suggestion. John John Fox, Professor

Re: [R] Selecting complementary colours

John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: Achim Zeileis [mailto:[EMAIL PROTECTED] Sent: Monday, May 21, 2007 7:07 PM

Re: [R] Is it possible to pass a Tcl/Tk component as argument to afunction

Dear Hao, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: Thursday, May 17, 2007 8:29 AM To: John Fox Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Is it possible to pass a Tcl/Tk component as argument to afunction John: Thanks for your reply, I

Re: [R] Is it possible to pass a Tcl/Tk component as argument to afunction

Dear Hao, You might take a look at how the Rcmdr package is implemented with many reusable elements. There is, for example, an initializeDialog function. I hope this helps, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario

Re: [R] Some questions on repeated measures (M)ANOVA mixed modelswith lme4

Dear Marco, You might also take a look at ?Anova (or ?Manova) in the car package; the last examples are for a repeated-measures ANOVA using both MANOVA and univariate approaches, the latter with GG and HF corrections. I hope this helps, John John Fox, Professor

Re: [R] Weighted least squares

Dear Adai, -Original Message- From: Adaikalavan Ramasamy [mailto:[EMAIL PROTECTED] Sent: Tuesday, May 08, 2007 8:38 PM To: S Ellison Cc: [EMAIL PROTECTED]; [EMAIL PROTECTED]; R-help@stat.math.ethz.ch Subject: Re: [R] Weighted least squares

Re: [R] Weighted least squares

Dear Hadley, -Original Message- From: hadley wickham [mailto:[EMAIL PROTECTED] Sent: Wednesday, May 09, 2007 2:21 AM To: John Fox Cc: R-help@stat.math.ethz.ch Subject: Re: [R] Weighted least squares Thanks John, That's just the explanation I was looking for. I had hoped

Re: [R] Weighted least squares

John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf

Re: [R] Library Package for Tobit regression

John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf

Re: [R] R News, volume 7, issue 1 is now available --error in AMMI article

- Biplots. Why is this an error? Regards, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From

Re: [R] R News, volume 7, issue 1 is now available --error in AMMI article

Dear Yuandan, The function definition given in the article doesn't produce an error for me. Regards, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] intersect more than two sets

intersection(a, b, c) [1] c d intersection(a, b, c, d) [1] d intersection(a, b, c, d, e) character(0) do.call(intersection, list(a, b, c)) [1] c d I hope this helps, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S

Re: [R] intersect more than two sets

, a, b) [1] b c d Fold(intersect, a, b, c) [1] c d Fold(intersect, a, b, c, d) [1] d Fold(intersect, a, b, c, d, e) character(0) do.call(Fold, list(intersect, a, b, c, d)) [1] d Regards, John John Fox, Professor Department of Sociology McMaster University

Re: [R] Fitting multinomial response in structural equation

well to find someone with whom you can about your research problem. Regards, John John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] how to convert the lower triangle of a matrix to a symmetricmatrix

John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf

Re: [R] Suggestions for statistical computing course

Dear Duncan and Giovanni, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Duncan Murdoch Sent: Friday, April 20, 2007 10:13 AM To: Giovanni Petris Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Suggestions for statistical computing course On

Re: [R] general question about plotting multiple regression results

read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox

Re: [R] Fit sem model with intercept

that it includes no covariance components; moreover, if this is the complete structural part of the model, then I think it is underidentified, and the two parts of the model (those involving eta1 and eta2) appear entirely separate. I hope this helps, John John Fox Department

Re: [R] Fit sem: problem of 'Error in solve.default(C[ind, ind]) : Lapack routine dgesv: system is exactly singular.'

research. John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED

Re: [R] Windows Vista issues

Dear Brian, I don't yet have a Vista machine, but it occurs to me that some of these problems might be avoided if R installed by default into c:\R rather than c:\Program Files\R. Is that the case? Thank you, by the way, for pursuing these issues. John John Fox

Re: [R] Windows Vista issues

Dear Brian, -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Saturday, April 14, 2007 9:44 AM To: John Fox Cc: [EMAIL PROTECTED] Subject: RE: [R] Windows Vista issues On Sat, 14 Apr 2007, John Fox wrote: Dear Brian, I don't yet have a Vista

Re: [R] creating a path diagram in sem

mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Fox Department of Sociology McMaster University Hamilton

Re: [R] sem is not taking the model

/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada http

Re: [R] sem vs. LISREL: sem fails

Dear Dimitri, -Original Message- From: John Smith [mailto:[EMAIL PROTECTED] Sent: Tuesday, April 10, 2007 8:11 AM To: John Fox Cc: r-help@stat.math.ethz.ch Subject: Re: [R] sem vs. LISREL: sem fails Dear Scott, (John, actually.) thanks a lot - all your comments are extremely

Re: [R] sem vs. LISREL: sem fails

model fit to the covariance matrix by setting the argument par.size to startvalues (as suggested in ?sem): sem.anxiety - sem(model, COVAR, N=150, par.size=startvalues) Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada

Re: [R] Dealing with large nominal predictor in sem package

will produce the error that you've encountered. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message

Re: [R] sem vs. LISREL: sem fails

, since you didn't appear to use anything in it; (2) comments in R are prefixed by #, not !. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] Dealing with large nominal predictor in sem package

Dear adschai, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: Monday, April 09, 2007 7:33 PM To: John Fox Cc: r-help@stat.math.ethz.ch Subject: Re: RE: [R] Dealing with large nominal predictor in sem package Hi John, Thank you. I think (2) from

Re: [R] Dealing with large nominal predictor in sem package

Dear adschai, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: Monday, April 09, 2007 8:30 PM To: John Fox; r-help@stat.math.ethz.ch Subject: Re: RE: [R] Dealing with large nominal predictor in sem package Hi John, Additional two questions

Re: [R] p value for coefficients in multinomial model

to contr.treatment; it calculates sums-of-squares in a manner that's independent of the parametrization, while Anova() does not. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http

Re: [R] about systemfit

the Wald tests)? I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED

Re: [R] about systemfit

. This suggests to me that the first place to look for a problem is in ftest.systemfit(). Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] MANOVA with repeated measurements

Dear Mark, Take a look at the Anova() function in the car package. There's an example of a repeated-measures MANOVA in ?Anova (the last example). Also see ?anova.mlm. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario

Re: [R] rcmdr help

hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] partial R

rownames(R) - colnames(R) - colnames(X) R } cr.plots() in the car package produces partial-residual plots (component+residual plots) but returns neither partial correlations nor partial-regression coefficients. Regards, John John Fox Department of Sociology

Re: [R] partial R

correlation between a particular X and Y controlling for all the other X's. This is exactly what the partial.cor() function computes. If I follow correctly what you want, you just have to square the partial correlations. Regards, John John Fox Department of Sociology

Re: [R] Bonferroni p-value greater than 1

that there are no unusually large studentized residuals (and indeed that the largest studentized residual is smaller than one would expect under the standard linear-model assumptions). I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario

Re: [R] fitted probabilities in multinomial logistic regression are identical for each level

. John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada http://socserv.mcmaster.ca/jfox/ __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Effect display of proportional odds model

). Second, the appendix to the paper and the originally posted examples didn't include the code for Figure 8 (which is Figure 10 in the published version of the paper). I'll add that. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario

Re: [R] Effect display of proportional odds model

provided do not construct graphs automatically, as those in the effects package do for linear and generalized linear models. Eventually, I'd like to incorporate this material in the effects package. Regards, John John Fox Department of Sociology McMaster University

Re: [R] how to get lsmeans?

to improve the appearance of the default graphs. But Brian's suggestion is very straightforward, and I guess that I shouldn't wait to implement it; I'll do so very soon. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905

Re: [R] how to get lsmeans?

John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Muenchen

Re: [R] how to get lsmeans?

Dear Frank, -Original Message- From: Frank E Harrell Jr [mailto:[EMAIL PROTECTED] Sent: Thursday, March 22, 2007 10:51 AM To: John Fox Cc: 'Prof Brian Ripley'; 'r-help'; 'Chuck Cleland'; 'Liaw, Andy' Subject: Re: [R] how to get lsmeans? John Fox wrote: Dear Frank

Re: [R] how to get lsmeans?

as the number of levels of the factors held constant increase. On the scale of the linear predictor, these sets would differ only by constants. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http

Re: [R] arrowhead styles

Dear Hendrik, Arrows() is the first entry turned up by RSiteSearch(arrows, restrict=functions), so its effort to hide isn't very successful. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604

Re: [R] Rcmdir does not work in SciVews-R

, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf

Re: [R] No fit statistics for some models using sem

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada http

Re: [R] sem: standardized covariance estimates

Dear Mike, I don't believe that I've provided a way to get these correlations automatically, but you can, of course, just divide the covariance by the product of the standard deviations. I hope this helps, John John Fox Department of Sociology McMaster

Re: [R] SEM - standardized path coefficients?

Dear Tim, See ?standardized.coefficients (after loading the sem package). Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] enhanced question / standardized coefficients

standardize the quantitative variables (Y and X's) in the regression equation. I hope that this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] Reference to dataframe and contents

Dear Rene, There are several ways to do what you want, including get(s1)[,S2]. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] efficient code. how to reduce running time?

thing that seems particularly striking in your results is the large difference between elapsed time and user CPU time, making me wonder what else was going on when you ran these examples. Regards, John John Fox Department of Sociology McMaster University Hamilton

Re: [R] efficient code. how to reduce running time?

Dear Brian, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Prof Brian Ripley Sent: Monday, January 22, 2007 11:06 AM To: Charilaos Skiadas Cc: John Fox; r-help@stat.math.ethz.ch Subject: Re: [R] efficient code. how to reduce running time

Re: [R] efficient code. how to reduce running time?

(). (BTW, is maxF really a matrix?) You also needlessly recompute max(F), never appear to use MSSid, and end lines with unnecessary semicolons. I hope that this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525

Re: [R] efficient code. how to reduce running time?

40.61NANA system.time(mods - lapply(as.list(X), function(x) lm(y ~ x))) [1] 53.29 0.37 53.94NANA In cases such as this, I don't even find the code using *apply() easier to read. Regards, John John Fox Department of Sociology McMaster University

Re: [R] extend summary.lm for hccm?

Dear Achim, -Original Message- From: Achim Zeileis [mailto:[EMAIL PROTECTED] Sent: Monday, December 25, 2006 8:38 AM To: John Fox Cc: 'Dirk Eddelbuettel'; 'ivo welch'; r-help@stat.math.ethz.ch Subject: Re: [R] extend summary.lm for hccm? On Sun, 24 Dec 2006, John Fox wrote

Re: [R] extend summary.lm for hccm?

] - table[,1]/table[,2] table[,4] - 2*pt(abs(table[,3]), df.residual(model), lower.tail=FALSE) sumry$coefficients - table print(sumry) cat(Note: Heteroscedasticity-consistant standard errors using adjustment, type, \n) } John Fox

Re: [R] extend summary.lm for hccm?

. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED

Re: [R] effect plot

Dear Ambrogi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ambrogi Federico Sent: Tuesday, December 19, 2006 12:03 PM To: r-help@stat.math.ethz.ch Subject: [R] effect plot Dear R users, Is there a simple way to use the effect function

Re: [R] unique columns of a matrix

Dear Roman, You can use unique(X, MARGIN=2). See ?unique for details. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] How to load Rcmd without Commander() ?

of the function in the global environment via reliability - Rcmdr:::reliability. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] Small Rcmdr issue

=no) } If anyone knows how I can fix the problem, I'd be happy to so do. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] Box Tidwell / Error Message / Error in parse(file, n, text, prompt) : syntax error in

. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto

Re: [R] Multilevel Modeling in R

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. John Fox Department of Sociology McMaster University Hamilton, Ontario, Canada

Re: [R] lmer and a response that is a proportion

Dear Cameron, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Cameron Gillies Sent: Sunday, December 03, 2006 1:58 PM To: r-help@stat.math.ethz.ch Subject: [R] lmer and a response that is a proportion Greetings all, I am using lmer (lme4

Re: [R] lmer and a response that is a proportion

, if you have many values of 0, 1, or both, then a transformation isn't promising (and, indeed, the logit wouldn't be defined for these values). In any event, I'd check diagnostics after a preliminary fit. I hope this helps, John John Fox Department of Sociology McMaster

Re: [R] Box Tidwell / Error Message

hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] error in hetcor function (polycor package)?

John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of BRENDAN KLICK

Re: [R] String question

Dear Carmen, length(grep([^01], bits)) == 0 should do the trick, returning TRUE if the string contains only 0's and 1's. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http

Re: [R] Random effects ANOVA

and covariance components, comparing each null model with the full model. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Re: [R] how to create this design matrix?

suggestion won't give you the answer, but I guess that at this point I should leave to someone else who understands more clearly what you want to sort it out. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905

Re: [R] how to create this design matrix?

Dear Michael, This looks like a multivariate simple regression -- that is, 12 response variables, one predictor. If the data are in the matrix X, then lm(X[,1:12] ~ X[,13]) should do the trick. I hope this helps, John John Fox Department of Sociology McMaster

Re: [R] Making a case for using R in Academia

, hosted by the Department of Politics and International Relations at Oxford, has mostly employed R and Stata. Of course, my own preference is for R. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604

Re: [R] poly() question

Dear Jonathan, An alternative is to use the raw argument to poly(), possibly after centering the predictor. This probably will work OK as long as the degree of the polynomial isn't high. I hope this helps, John John Fox Department of Sociology McMaster

[R] correlation argument for lmer?

Dear r-help members, Can lmer() in the lme4 package fit models that have a specified within-group correlation structure, as provided, for example, by the correlation argument to lme() in the nlme package? Thanks, John John Fox Department of Sociology McMaster

Re: [R] coefficient constraints

hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] Error using White-corrected Anova (white.adjust=hc3)

several days. Sorry, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto

Re: [R] Qualitative Data??(String command)

Dear Davia, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Davia Cox Sent: Friday, October 27, 2006 3:37 PM To: r-help@stat.math.ethz.ch Subject: [R] Qualitative Data??(String command) I am using the read.table function to load an Excel data

Re: [R] Recoding categorical variables

Dear Murray, How about as.numeric(factor(y)) ? I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original

Re: [R] Help with recode and factor functions

(see ?car) that has several examples of the use of the function. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original

Re: [R] How to get multiple Correlation Coefficients

upload the new version of the package to CRAN shortly. Regards, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original

Re: [R] Differences in estimates calculated in R and SPSS

about how R handles factors in linear models). I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original

Re: [R] regression analyses using a vector of means anda variance-covariance matrix

. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto

  1   2   3   4   5   6   7   >