search for: yuelin

...a vector: letters[1:10]. I want to get a matrix of 120 rows that contains all possible subsets of 3 letters out of the 10. And I would like to be able to generalize to any r out of n, like function(x=letters[1:10], r=3) that returns: a, b, c a, b, d a, b, e ... etc. Many thanks in advance, - Yuelin Li. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at...

...tin_va_alist' undeclared (first use in this function) errors.c: In function `Rf_WarningMessage': errors.c:562: `__builtin_va_alist' undeclared (first use in this function) *** Error code 1 make: Fatal error: Command failed for target `errors.o' Current working directory /export/home/yuelin/R-1.1.0/src/main *** Error code 1 make: Fatal error: Command failed for target `R' Current working directory /export/home/yuelin/R-1.1.0/src/main *** Error code 1 make: Fatal error: Command failed for target `R' Current working directory /export/home/yuelin/R-1.1.0/src *** Error code 1 make...

...ntially meaningless once the experiment has been done". But they seem to be used, and sometimes encouraged (e.g., http://www.informaworld.com/smpp/content~content=a787469931). I do not intend to start a debate. I simply would like to learn more. Can anyone suggest a few articles? Thanks, Yuelin. Tukey JW. Tightening the clinical trial. Control Clin Trials. 1993;14:266-285. ===================================================================== Please note that this e-mail and any files transmitted with it may be privileged, confidential, and protected from disclosu...

.../local/lib/libreadline.a,libreadline.so@, and libreadline.so.4. The attached errors are from config.log. ./configure finds readline header files in /usr/local/include/readline. I can make and install R, but command line editing and savehistory() do not work. Suggestions are appreciated. Yuelin Li. -------- configure:11366: checking for rl_callback_read_char in -lreadline configure:11397: gcc -o conftest -g -O2 -I/usr/local/include -L/usr/local/lib conftest.c -lreadline -ldl -lncurses -lm >&5 ld: fatal: symbol `_init' is multiply-defined: (file /usr/local/lib/gcc-l...

In the hypothetical example below, how do I add two segments() into the two panels, respectively? Say segments(x0=5, y0=10, x1=5, y1=20) on the left and segments(x0=15, y0=-10, x1=15, y1=-2) on the right? Many thanks in advance, Yuelin Li. ps. part of the code came from a solution given by Deepayan Sarkar. ------------------- library(lattice) set.seed(12345) x <- 0:20 y.male.obs <- - 1.2 * x + 22 + rnorm(length(x), sd = 3) y.male.prd <- - 1.2 * x + 22 y.fema.obs <- - 2.2 * x + 30 + rnorm(length(x), sd = 2) y.fema.pr...

Why is as.numeric("3/2") NA plus a warning? Thanks, Mike. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch

...vel is profound disability. My questions are: 1. How do I interpret the coefficients in the third model (the intercept is removed)? 2. Does glm3 make sense? Does anyone ever want to construct a model like glm3? If so, when? 3. What is the algebraic specification of glm3? Many thanks, -- Yuelin Li. ------------------- categ <- factor(rep(c(1, 2, 3), times=c(20, 20, 20)), labels=c("P", "S", "M")) Y <- c(1,1,1,1,1,1,2,1,2,1,2,1,1,1,1,1,1,2,2, 1,2,2,1,2,2,1,2,1,2,1,1,2,2,2,2,2,2,2,2,1, 2,2,2,2,2,1,2,1,1,1,2,2,2,2,1,2,1,1,1,2...

...to approximate statistical power for a univariate test in repeated-measures ANOVAs when the sphericity assumption is not met. If not, perhaps the MANOVA approach by O'Brien & Kaiser (1985, Psych Bull: 97, 316-33). library(pwr) does not seem to support either. Many thanks in advance, Yuelin Li. ===================================================================== Please note that this e-mail and any files transmitted with it may be privileged, confidential, and protected from disclosure under applicable law. If the reader of this message is not the inten...

I have been trying to analyse an unreplicated repeated measures 2-level factorial design with 11 factors using aov() with Error(), similar to that described in "Notes on the use of R for psychology experiments and questionnaires" by Jonathan Baron and Yuelin Li (see the example of Hay's) I have found that there seems to be a limit, in the number of factors in Error() . For example, the following worked : aov(CV ~(Twb+Plb+Psb+tinsb+tseedb+Tmb+Mcb+Mom0b+Mom1b+Mom2b+Mom3b)^3 + Error(Batch / (Twb+Plb+Psb+tinsb+tseedb+Tmb+Mcb)),simulation.6batches...

...r MI, then how do I get the MI-adjusted coefficients and standard errors? I can't seem to find answers in the documentations, including rereading section 8.10 of the Harrell (2001) book Googling located a thread in R-help back in 2003, which seemed dated. Many thanks in advance for the help, Yuelin. http://idecide.mskcc.org ------------------------------- > library(Hmisc) Loading required package: survival Loading required package: splines > data(kyphosis, package = "rpart") > kp <- lapply(kyphosis, function(x) + { is.na(x) <- sample(1:length(x), size = 10); x...

...t;>>> 2-level >>>>> factorial design with 11 factors using aov() with Error(), similar to >>>>> that described in "Notes on the use of R for psychology >>>>> experiments and >>>>> questionnaires" by Jonathan Baron and Yuelin Li (see the example of >>>> >>> Hay's) >> >>>>> >>>>> I have found that there seems to be a limit, in the number of >>>>> factors in >>>>> Error() . For example, the following worked : >>>>&g...

......, dependence=TRUE). The script complains about not finding "/usr/local/lib/R/library/rggobi/libs/*", but that directory is there, with one file "rggobi.so" (possibly from an earlier successful compilation). Any suggestions on what I am doing wrong? Many thanks in advance. Yuelin. --------- R output --------------- > install.packages("rggobi", repos = "http://lib.stat.cmu.edu/R/CRAN", dependencies = TRUE, clean = TRUE) trying URL 'http://lib.stat.cmu.edu/R/CRAN/src/contrib/rggobi_2.1.6.tar.gz' Content type 'application/x-gzip' length...

I just successfully compiled R-1.6.0 on a Red Hat 7.2 machine. When I tried to install "grid" from source, I got this error. Any suggestions? I don't think it is a dependency problem. I also attach the output from library(). Yuelin Li. ------------- > install.packages("grid") trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 84585 bytes opened URL .......... .......... .......... .......... .......... .......... .......... .......... .. downl...

For those new to this list, I point out that Yuelin Li and I have drafted an introduction to R for psychology, trying to emphasize the things the psychologists usually do. It is at http://www.psych.upenn.edu/~baron/rpsych.pdf and http://www.psych.upenn.edu/~baron/rpsych.htm Because it is a draft, with frequent changes, we have not yet asked for it...

We now have a draft of: Notes on the use of R for psychology experiments and questionnaires by Jonathan Baron and Yuelin Li at http://www.psych.upenn.edu/~baron/rpsych.htm and http://www.psych.upenn.edu/~baron/rpsych.pdf and http://www.psych.upenn.edu/~baron/rpsych.tex It is intended for students and others who are doing research in psychology. What makes it "for psychology" is that it tries to cover the...

...--> 0 (better) 6 12 ---------------------------------------------------- A fisher.test(matrix(c(2,6,0,12), ncol=2)) shows a p-value of 0.15. I'd appreciate suggestions to alternative methods, perhaps a test of conditional independence in loglin()? But I am not sure how to do that. Yuelin Li. --------- Table ------------- Control Group (baseline by 6 months) Frequency| Percent | Row Pct | 0| 1| Total ---------+--------+--------+ 0 | 20 | 2 | 22 | 71.43 | 7.14 | 78.57 | 90.91 | 9.09 | ---------+--------+--------+...

Yuelin Li and I have made some revisions in our introductory document, "Notes on the use of R for psychology experiments and questionnaires," which is still available through my R page (below) in html or pdf. We fixed typos. We updated and organized the list of commands. ("Re"-organi...

...a coefficients associated with items i5 to i17, but summary() does not automatically summarize the log odds associated with the strata "id". I suppose the log odds associated with strata(id) gives me \alpha_k. But I can't seem to get them with summary(). Help is greatly appreciated. Yuelin Li. ---- my R syntax is ----------- > exam2.clog<- clogit(resp ~ i5+i6+i7+i8+i9+i10+i11+i12+i13+i14+ i15+i16+i17+ strata(id), data=exam1.1) > summary(exam2.clog) Call: clogit(resp ~ i5 + i6 + i7 + i8 + i9 + i10 + i11 + i12 + i13 + i14 + i15 + i16 + i17 +...

I'd like to build R using my optimized blas and lapack libraries. It seems know matter what I do, the configure script uses the blas supplied with the source. My blas and lapack libraries are in /usr/local/atlas/lib. How can I get configure to use these? [[alternative HTML version deleted]]

I am totally confused as to how to use anova. I have three vectors and would like to use anova on them but I don't understand how lm or glm comes into play. In matlab, you just give the three vectors. Why isn't it the same in R? Any help would be greatly appreciated. Anna