Search Archives

Dear R fans, I wrote a trivial wordpress plugin for users of wordpress and Rweb (or other web interface ) Plugin URL: http://lixiaoxu.lxxm.com/rwebfriend/ Xiaoxu

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20081223/9be0081f/attachment.pl>

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20081224/13618218/attachment.pl>

Dear R fans, I am preparing a lecture discussing the paper of Wilkinson and APA Task Force on Statistical Inference (1999, American Psychologist, 54, 594-604.) I'd like to demo their Figure 3B (p. 602) with R. That is...

Dear R fans, MASS::mvrnorm with empirical=T is very useful to reproduce practice data of published papers. Now I want to reproduce a binomial variable with multiple normal variables. The target paper gave a big corr table with sd...

You can choose to save your function in the file ..\etc\Rprofile.site Xiaoxu On Fri, Oct 31, 2008 at 1:45 AM, Douglas M. Hultstrand <dmhultst at metstat.com> wrote: > Hello all, > > I am fairly new to R and...

##Q1. confint.glm(...) fails for an example of HSAUR data("womensrole", package = "HSAUR"); ## summary(womensrole); womensrole_glm_2 <- glm(fm2, data = womensrole,family = binomial()) ## summary(womensrole_glm_2); confint(womensrole_glm_2); ## -------Fail--------- # Waiting for profiling to be done... # Error...

I think you should try one X for A/X=sum(exp(c_i/d -log(X))) The optional X could be as the following ... ## to test ### C = c(c_i/d; i=1,2,...,n) n<-1000; C<-runif...

I should have read the following page on R_Extension_for_MediaWiki http://mars.wiwi.hu-berlin.de/mediawiki/sk/index.php/R_Extension_for_MediaWiki_v0.06#New_tags_and_attributes Has anybody seen an <Rform>...</Rform> online example...

## Use model.matrix ## Data is the same ## continue m <- model.matrix(lm(rep(0,length(y)) ~ disease + drug +disease:drug)); ## Model.matrix(lm(y~...)) will drop is.na(y) rows. That result will be Type II rather than III ## for...

Ajay and R fans, I'll never think my several lines of simplest PHP codes can be something great. They only shortcut Rweb and wordpress. But, SaaS is really wonderful. If only some R web interface can be smartly deployed...

## I got it. IV(s) of interaction should be orthogonal to main effect IV(s). ## Type III ANOVA / Interaction alone x_interaction<-cbind( (drug==2)&(disease==2) ,(drug==3)&(disease==2) ,(drug==4)&(disease==2) ,(drug==2)&(disease==3) ,(drug...

## Question1: How to define IV with interaction alone, without main effects? ## Question2: Should Type III ANOVA in package car be independent of the factor level order? ## data from http://www.otago.ac.nz/sas/stat/chap30/sect52.htm drug <- c...