[R] cor matrix in multivariate regression

2013-10-21 Thread Suyan Tian
. . . . …. . .. . 0 0 1 0.8 0 0 0.8 1 Does anyone know how to do it? Thanks, Suyan Tian Associate Professor The Jilin University __ R-help@r-project.org mailing list

[R] read wiff extension files into R

2012-03-21 Thread Suyan Tian
Does anyone know how to read or import a wiff (AB SCIEX windows interchange file format) extension into R? Thanks a lot. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-pro

[R] skip the error to continue the logistic regression in a loop

2009-07-02 Thread Suyan Tian
Hi, everyone: I am running logistic regression on a bunch of variables using apply command. But an error occurs, the whole process stops. I am wondering if anyone knows how to skip this error and to continue the regression for the rest of variable. What I did is that first confine a func

[R] save the result into a word file

2009-07-02 Thread Suyan Tian
Hi, everyone: I forget on how to save a result from R as a word document. For example, if I run the linear regression and want to save the result: summary(lm(y~x) in a word file. So I can show it later to my clients. Any idea on how to do this? Thanks a lot, Suyan ___

[R] quantile regression using SemiPar package

2009-01-15 Thread Suyan Tian
Hi everyone: I want to fit the truncated polynomial smoothing to the quantiles instead of means, does someone know how to do it? I am thinking that maybe I can use SemiPar package, but can not find how. Thanks so many, Suyan __ R-help@r-project.

[R] multiple plots in R

2008-11-23 Thread Suyan Tian
Hi, I just try to draw multiple plots in one page using R, I used par command. For example I have 7 plots, but instead of arranging them in the default way plot1plot2 plot3 plot4 plot5 plot6 plot7 I want them in this order plot1 plot2 plot3 plot4 plot5 plot6 plot7

[R] change data into array

2008-10-28 Thread Suyan Tian
Sorry to bother everyone again with one question about array. I have one data set saved in this way patienttime x1 x2 response 1 1 18 1 1 2 27 1 2 1 36 1 2

[R] abind

2008-10-28 Thread Suyan Tian
I am trying to combine two arrays with different dimensions into one. For example The first one is 1 2 3 4 5 6 The second one is 7 8 9 10 The resulted one would be like , , 1 1 2 3 4 5 6 , , 2 7 8 9 10 I used abind to do this, but failed. Could somebody please let me know how to do thi

[R] boxplot

2008-08-22 Thread Suyan Tian
Hi, I just made a boxplot but I want to change the thickness of the line used for the median to look a little thinner. Could anyone please help me figuring out the R-code to do it? Thanks, I really appreciate it. Suyan __ R-help@r-project.org mail

[R] combine two data-frames into one

2008-03-14 Thread Suyan Tian
Hi Everyone: I need some help on combining two datasets into one. Since for those two datasets, there are different variables. Like In dataset 1, I have pt x1 x2 1 12 24 2 20 17 In dataset 2, I have pt x2 x3 3 19 33 4 20 16 The final dataset I want should look as pt x1 x2 x3

[R] kronecker product in mixed model

2008-02-25 Thread Suyan Tian
Hi, does anyone know how to self-define the variance-covariance structure of the within-subject variance for the mixed models , i.g, yi=Xibeta+Zibi+ei For ei I want to use kronecker product so its variance matrix will be R=V*E, here * represents Kronecker product. Thanks so many.

[R] latent class analysis

2008-02-01 Thread Suyan Tian
Hi, does some one know how to run the latent class analysis in R? like which package and command I should use? thanks a lot, Suyan Tian __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] bootstrapping on the growth curve

2007-12-05 Thread Suyan Tian
Hi, I am trying to get 95% CI s around a quantile growth curve, but my result looks strange to me. Here is the function I defined myself boot.qregress<-function(mat1, group, quantile, int, seed.1){ boot.fit<-NULL set.seed(seed.1) for (i in 1:int){ index<-sample((unique(ma