[R] R cannot be started

2010-09-07 Thread Tao
I used Rtools for installing a package under Windows XP. I used commands like R CMD INSTALL (build, check) in the CMD window. It worked well until I found that my revised package can be updated with R CMD INSTALL. After I restarted my computer, R cannot be started any more. I double click the R

[R] Package build and install under Windows

2010-09-14 Thread Tao
Any other method besides using RTOOLS? I am using R 2.10.1, which version of Rtools should I use? It seems that there are two versions on the website. The commands like R CMD build should be used in the CMD window, right? Thanks. -- View this message in context:

Re: [R] Package build and install under Windows

2010-09-15 Thread Tao
Thanks. It seems that both Rtools210.exe and Rtools211.exe support R 1.10.x, which one is better for R 2.10.1. -- View this message in context: http://r.789695.n4.nabble.com/Package-build-and-install-under-Windows-tp2539509p2540713.html Sent from the R help mailing list archive at Nabble.com.

[R] date calculation

2010-10-29 Thread Shi, Tao
Hi list, Could someone explain to me why the following result is not a integer? difftime(strptime(24NOV2004, format=%d%b%Y), strptime(13MAY2004, format=%d%b%Y), units=days) Time difference of 195.0417 days I'm using R2.12.0 on WinXP. Thanks! ...Tao

Re: [R] date calculation

2010-10-30 Thread Shi, Tao
Hi Phil, Thanks for the reply, but I don't think you have explained where the decimal part is coming from ...Tao - Original Message - From:Phil Spector spec...@stat.berkeley.edu To:Shi, Tao shida...@yahoo.com Cc:r-h...@stat.math.ethz.ch Sent:Friday, October 29, 2010 5:04:23

Re: [R] date calculation

2010-10-30 Thread Shi, Tao
Hi Ben, That must be the case! In fact if I do: difftime(strptime(24NOV2004, format=%d%b%Y), strptime(13MAY2004,format=%d%b%Y), units=days, tz=GMT) Time difference of 195 days which supports your claim. Can someone from the R development team confirm this? Thanks! ...Tao

Re: [R] date calculation

2010-10-30 Thread Shi, Tao
Hi Ben and Josh, No, I don't think that there is a bug. It's just due to my lack of knowledge on Date class, I was simply amazed that it has the daylight saving time built in at first. But on the second thought, it really should. Thank you all for the explanations! ...Tao From: Ben

Re: [R] date calculation

2010-10-30 Thread Shi, Tao
I think this demonstrate on of the differences between Class 'Date' and 'POSIXlt'. Thanks, Marc! ...Tao From: Marc Schwartz marc_schwa...@me.com To: Ben Bolker bbol...@gmail.com Cc: Shi, Tao shida...@yahoo.com; r-h...@stat.math.ethz.ch r-h...@stat.math.ethz.ch Sent: Saturday, October 30

[R] calculating martingale residual on new data using predict.coxph

2010-11-19 Thread Shi, Tao
' and 'predict.Design', but from my reading so far, I'm not sure they can do that. Do you other ways to calculate martingale residuals on a new data? Thank you very much! ...Tao __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] calculating martingale residual on new data using predict.coxph

2010-11-19 Thread Shi, Tao
Hi David, Thank you for the quick reply! resid(fit) only gives the residuals on the training data not on test data. ...Tao - Original Message From: David Winsemius dwinsem...@comcast.net To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org; dieter.me...@menne-biomed.de

Re: [R] calculating martingale residual on new data using predict.coxph

2010-11-21 Thread Shi, Tao
Hi David, Thanks, but I don't quite follow your examples below. The residuals you calculated are still based on the training data from which your cox model was generated. I'm interested in the testing data. Best, ...Tao - Original Message From: David Winsemius dwinsem

[R] question on uniCox

2010-11-22 Thread Shi, Tao
a univariate Cox model, the features that give NA beta estimates are actually pretty significant. Could you please let me know what happened and how to avoid this? I’ve attached the outputs of the function calls below. Thank you very much! ...Tao a - uniCox(x=t(dat.ave.train.base), y

Re: [R] calculating martingale residual on new data using predict.coxph

2010-11-22 Thread Shi, Tao
Thank you, Terry! - Original Message From: Terry Therneau thern...@mayo.edu To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org; dieter.me...@menne-biomed.de; r_ting...@hotmail.com Sent: Mon, November 22, 2010 6:11:15 AM Subject: Re: calculating martingale residual on new

Re: [R] calculating martingale residual on new data using

2010-11-22 Thread Shi, Tao
Thank you for the advice, Frank! ...Tao - Original Message From: Frank Harrell f.harr...@vanderbilt.edu To: r-help@r-project.org Sent: Sun, November 21, 2010 5:49:36 AM Subject: Re: [R] calculating martingale residual on new data using The tendency is to use residual-like

Re: [R] calculating martingale residual on new data using predict.coxph

2010-11-22 Thread Shi, Tao
Thank you, David! You've been always helpful! ...Tao - Original Message From: David Winsemius dwinsem...@comcast.net To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org; dieter.me...@menne-biomed.de; r_ting...@hotmail.com Sent: Sun, November 21, 2010 5:50:31 AM Subject: Re

[R] given a sensitivity calculate specificity based on a ROC curve

2010-11-29 Thread Shi, Tao
Hi list, I know this is not hard to implement based on the returned objects from ROC, ROCR or a couple of other roc-related packages. I'm just wondering if there is already such a function exist. Thanks! ...Tao __ R-help@r-project.org mailing list

[R] pairs with same xlim and ylim scale

2010-08-31 Thread Shi, Tao
, they didn't help much. I think this is b/c I have multiple panel functions. Help! Thanks! ...Tao pairs1 - function(x, ...) { f.xy - function(x, y, ...) { points(x, y, ...) abline(0, 1, col = 2) } panel.cor - function(x, y, digits=2, prefix=, cex.cor

Re: [R] pairs with same xlim and ylim scale

2010-09-02 Thread Shi, Tao
Hi Dejian, You're right on this! Do you know how to pass those two argument into lower.panel? Thanks! ...Tao From: Dejian Zhao zha...@ioz.ac.cn To: r-help@r-project.org Sent: Tue, August 31, 2010 6:10:16 PM Subject: Re: [R] pairs with same xlim and ylim

[R] a minor bug in venn from gplots?

2009-10-29 Thread Tao Shi
from the fact there are one grid.newpage() call and two plot function calls when numCircles==4 (see below).  I wonder the grid.newpage() and the second plot call are necessary? Thanks a lot. ...Tao ##=     else

[R] Can't load doSMP from REvolutionR in regular R2.11.0

2010-04-28 Thread Tao Shi
R peripherals were set up to work with the regular R2.11.0.  So, I really want to make this work.  Anyideas? Many thanks in advance! ...Tao library(doSMP) Loading required package: foreach Loading required package: iterators Loading required package: codetools foreach: simple, scalable

Re: [R] Can't load doSMP from REvolutionR in regular R2.11.0

2010-04-29 Thread Tao Shi
Hi David, Thank you for the reply!  Do you know where I can find the source code for these packages?  I can give it a try. ...Tao From: da...@revolution-computing.com Date: Thu, 29 Apr 2010 08:59:08 -0700 Subject: Re: [R] Can't load doSMP from

Re: [R] Can't load doSMP from REvolutionR in regular R2.11.0

2010-04-29 Thread Tao Shi
Thanks, David!  I forgot to check that email... I have built the revoIPC using R2.11.0 and it seemed there were no error messages.  I can load doSMP now, but haven't tested it yet. Tal, I'll send you the file offline, so you can also test it. best! ...Tao

[R] apply vs. foreach vs. foreach with doSMP (multi cores)

2010-04-29 Thread Tao Shi
could you please explain the error message (i.e. Error in ipcTaskSetEnvironment(taskq, envir) : ...) in the second run? The results were recorded on REvolution R 3.2 and I observed the same on regular R2.11.0. Many thanks! ...Tao library(doSMP) Loading required package: foreach Loading

[R] 'matplot' for matrix with NAs: broken lines

2010-05-05 Thread Tao Shi
alternatives so I don't have to use 'plot' + 'lines' and loop through all the columns. Many thanks! ...Tao set.seed(1234) a=b=matrix(rnorm(9), 3,3) b[2,2]=NA matplot(a, type=b) matplot(b, type=b)  ## I want the two 2 connected! matplot(b, type=l)  ## Now my data for the second column are missing from

Re: [R] Read data from .csv file as a matrix

2010-05-05 Thread Tao Shi
) should work. ...Tao I have a csv file that contains weather observation (rows) by days (in columns). I open using: temp = read.csv(Weather.csv, sep=,) and read: X X1.Jan X2.Jan X3.Jan X4.Jan 1 Min 2 3 4 1 2 Max

[R] What is the best way to have R output tables in an MS Word format?

2010-05-05 Thread Tao Shi
. Harrell's mentioned on his webpage. http://www.freewaregenius.com/2010/03/06/how-to-convert-pdf-to-word-doc-for-free-a-comparative-test/ Also, Adobe Acrobat 9.0 can do PDF-Word, but I haven't tried it personally. ...Tao On Fri, Apr 30, 2010 at 5:13 PM, Max Gunther max.gunt...@vanderbilt.eduwrote

Re: [R] 'matplot' for matrix with NAs: broken lines

2010-05-05 Thread Tao Shi
sometimes there is only one non-missing value in a particular column and na.approx can't work (well I could selectively impute the NAs ... ) So far, my best solution to this is to use xyplot.  It does this by default, but of course I need some data manipulation first. ...Tao

Re: [R] 'matplot' for matrix with NAs: broken lines

2010-05-06 Thread Tao Shi
to keep track of both trend and missing values.  The original post was really asking whether a simple change of some parameters in matplot can do this.  Now, I guess not. ...Tao From: maech...@stat.math.ethz.ch Date: Thu, 6 May 2010 18:34:22 +0200 To: shi

Re: [R] 'matplot' for matrix with NAs: broken lines

2010-05-06 Thread Tao Shi
values are missing. HTH, Dennis On Thu, May 6, 2010 at 10:44 AM, Tao Shi wrote: I just found out that my does this by default statement (by which I was referring to the ability to automatically connect two points with a NA in the middle in a time series) is wrong! Actually, all plotting

[R] a question about latex in Hmisc

2010-05-12 Thread Shi, Tao
on WinXP. Thanks! ...Tao __ R-help@r-project.org 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

Re: [R] Whiskers on the default boxplot {graphics}

2010-05-12 Thread Shi, Tao
Jason, All these are clearly defined in the help file for 'boxplot' under 'range'. Don't understand how you missed that. ...Tao - Original Message From: Jason Rupert jasonkrup...@yahoo.com To: Dennis Murphy djmu...@gmail.com Cc: R Project Help R-help@r-project.org Sent: Wed

Re: [R] a question about latex in Hmisc

2010-05-12 Thread Shi, Tao
. It actually worked for me once (i.e. I got Yap poped out to preview the table I'm generating), but I had no idea what I did differently. ...Tao - Original Message From: Erik Iverson er...@ccbr.umn.edu To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org Sent: Wed, May 12

Re: [R] a question about latex in Hmisc

2010-05-12 Thread Shi, Tao
Hi Richard, Obviously, the list doesn't like my attachment. Here it is: http://i41.tinypic.com/15qz387.jpg I set the path using Windows control panel. And didn't work means I'm still getting the same error window shown above. ...Tao - Original Message From: Richard M

Re: [R] a question about latex in Hmisc

2010-05-12 Thread Shi, Tao
(rather a separate .tex file just for that table). Am I right? Thanks, ...Tao - Original Message From: Erik Iverson er...@ccbr.umn.edu To: Ista Zahn istaz...@gmail.com Cc: r-help@r-project.org Sent: Wed, May 12, 2010 12:34:15 PM Subject: Re: [R] a question about latex in Hmisc

Re: [R] a question about latex in Hmisc

2010-05-12 Thread Shi, Tao
The .dvi file is being created in the working directory for sure, not sure whether it's in the tmp folder or not. (I just tested on my home machine which only has R but no MikTex installed, I can see only .tex file in the tmp folder. No .log file in the tmp folder.) ...Tao

Re: [R] Path to R script

2010-05-12 Thread Shi, Tao
group. ...Tao - Original Message From: Johannes W. Dietrich j.w.dietr...@medizinische-kybernetik.de To: r-help@r-project.org Sent: Wed, May 12, 2010 12:31:19 PM Subject: [R] Path to R script There is certainly a trivial solution for my question, but I can't find the answer

Re: [R] a question about latex in Hmisc

2010-05-12 Thread Shi, Tao
of it. ...Tao = If you would rather have a separate file for your table you can continue to use your original method, e.g., = tmp - latex(x, file=x.tex) @ \include{x.tex} (not sure what your objection to this was in the first place). -Ista

Re: [R] Whiskers on the default boxplot {graphics}

2010-05-12 Thread Shi, Tao
Hi Peter, You're absolutely correct! The description for 'range' in 'boxplot' help file is a little bit confusing by using the words interquartile range. I think it should be changed to the length of the box to be exact and consistent with those in the help file for boxplot.stats. ...Tao

Re: [R] Whiskers on the default boxplot {graphics}

2010-05-13 Thread Shi, Tao
, as I stated earlier. With all these being said, this is up to the R core team to decide. ...Tao - Original Message From: Robert Baer rb...@atsu.edu To: Shi, Tao shida...@yahoo.com; Peter Ehlers ehl...@ucalgary.ca Cc: R Project Help R-help@r-project.org Sent: Thu, May 13, 2010 7:25

Re: [R] a question about latex in Hmisc

2010-05-13 Thread Shi, Tao
directory. So when YAP was called, it was looking for the .dvi file in the R temp folder. Therefore, the error window. http://i41.tinypic.com/15qz387.jpg I wonder how to fix that? Thanks! ...Tao - Original Message From: Duncan Murdoch murdoch.dun...@gmail.com To: Shi, Tao

Re: [R] a question about latex in Hmisc

2010-05-13 Thread Shi, Tao
my 'big' .tex file looks, as it's just a intermediate and I'll do all my editing in the .Rnw file. So, I think I'll use the latex(x, file=) method more. This also will make my working directory cleaner. Thank you very much! ...Tao - Original Message From: Frank E Harrell Jr

Re: [R] Path to R script

2010-05-13 Thread Shi, Tao
In this case, deploy them as a R package is the way to go. - Original Message From: Johannes W. Dietrich j.w.dietr...@medizinische-kybernetik.de To: Shi, Tao shida...@yahoo.com; johannes.dietr...@rub.de Cc: r-help@r-project.org Sent: Thu, May 13, 2010 9:14:24 AM Subject: Re: [R

Re: [R] Multiple plots; single x(y) labels

2010-05-13 Thread Shi, Tao
David said, you're re-inventing the wheels . ...Tao - Original Message From: Xin Ge xingemaill...@gmail.com To: David Winsemius dwinsem...@comcast.net Cc: r-help@r-project.org Sent: Thu, May 13, 2010 10:24:56 AM Subject: Re: [R] Multiple plots; single x(y) labels Hi David

Re: [R] read table for Fisher Exact

2010-05-13 Thread Shi, Tao
Hi Prof. Harrell, Could you please elaborate on why chi-square test is more appropriate in this case? Thank you very much! ...Tao - Original Message From: Frank E Harrell Jr f.harr...@vanderbilt.edu To: r-help@r-project.org Sent: Thu, May 13, 2010 5:35:07 AM Subject: Re: [R

Re: [R] read table for Fisher Exact

2010-05-13 Thread Shi, Tao
Thanks, David! - Original Message From: David Winsemius dwinsem...@comcast.net To: Frank E Harrell Jr f.harr...@vanderbilt.edu Cc: Shi, Tao shida...@yahoo.com; r-help@r-project.org Sent: Thu, May 13, 2010 8:29:31 PM Subject: Re: [R] read table for Fisher Exact On May 13, 2010

Re: [R] SAS for R-users

2010-05-14 Thread Shi, Tao
Thomas, If you're thinking to leverage your R programming skill in learning SAS, you'll be disappointed, as the two have quite different grammar. The book Erik mentioned is a good start. After that, just reading the SAS help file, which is pretty comprehensive, will keep you busy. ...Tao

[R] rpart: how to use each variable only once?

2010-05-14 Thread Shi, Tao
using X1, I guess it's still a reasonable thing to do, right? Many thanks! ...Tao __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Multiple plots; single x(y) labels

2010-05-14 Thread Shi, Tao
Gurmeet, I think Xin is more interested in the common axes, rather than just one single xlab or ylab. Jim's solution is much more fancier than mine :-) ...Tao - Original Message From: Gurmeet r.emailing.l...@gmail.com To: Jim Lemon j...@bitwrit.com.au Cc: r-help@r-project.org

Re: [R] Putting 6 graphs on one page

2010-05-14 Thread Shi, Tao
Without the actual data, it's hard to see what's going on here, but It seems you have to restructure your data object to a long table, then it should be easy to use 'dotplot' to generate your plots. ...Tao - Original Message From: Dimitri Liakhovitski dimitri.liakhovit

Re: [R] Multiple plots; single x(y) labels

2010-05-14 Thread Shi, Tao
ha, I was focusing on the wrong thing! Sorry, Gurmeet. Good job! From: Xin Ge xingemaill...@gmail.com To: Shi, Tao shida...@yahoo.com Cc: Gurmeet r.emailing.l...@gmail.com; Jim Lemon j...@bitwrit.com.au; r-help@r-project.org Sent: Fri, May 14, 2010 10:51:26 AM Subject: Re: [R] Multiple

Re: [R] rpart: how to use each variable only once?

2010-05-14 Thread Shi, Tao
This is just one single decision tree, not forest. If you're asking what package I use to construct one single tree, it's 'rpart'. From: Changbin Du changb...@gmail.com To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org Sent: Fri, May 14, 2010 12:35:20 PM Subject: Re: [R] rpart: how

Re: [R] automate curve drawing on nls() object

2010-05-18 Thread Shi, Tao
I can't directly answer your question regarding 'expression', but can you just replace b, c,d, and e with coef(obj)[1], coef(obj)[2], ... etc. You still can automate the whole process this way, right? - Original Message From: array chip arrayprof...@yahoo.com To:

Re: [R] automate curve drawing on nls() object

2010-05-18 Thread Shi, Tao
In this case, Ben's approach is the way to go. I'm curious how you fit nls without knowing the model formula beforehand? ...Tao - Original Message From: array chip arrayprof...@yahoo.com To: r-help@r-project.org; TaoShi shida...@yahoo.com Sent: Tue, May 18, 2010 5:22:47 PM

Re: [R] save in for loop

2010-05-19 Thread Shi, Tao
Ivan, Try this: eval(parse(text=paste(save(file, i, , file=\file, i, .RData\), sep=))) ...Tao - Original Message From: Ivan Calandra ivan.calan...@uni-hamburg.de To: r-help@r-project.org Sent: Wed, May 19, 2010 7:56:44 AM Subject: [R] save in for loop Dear users, My problem

[R] colored venn diagram

2010-05-19 Thread Shi, Tao
Hi list, This is probably too much to ask, but I'm wondering if there is a ready-to-use function somewhere that allows me to color one area of a venn diagram (e.g. the intersection of two sets)? Thanks! ...Tao __ R-help@r-project.org mailing list

Re: [R] p-values 2.2e-16 not reported

2010-05-20 Thread Shi, Tao
Will, I'm wondering if you have any insights after looking at the cor.test source code. It seems to be fine to me, as the p value is either calculated by your first method or a .C code. ...Tao - Original Message From: Will Eagle will.ea...@gmx.net To: r-help@r-project.org

Re: [R] About the breakpoint when making heatmap with lots of variables

2010-05-20 Thread Shi, Tao
The breakpoint you mentioned is irrelevant here. Clustering 15672 genes (I assume this is a microarray data) requires lots of memory. I suggest you either filter your gene list down to thousands or just plot the column dendrogram without showing the heatmap. plot(hclust(dist(x))) ...Tao

Re: [R] colored venn diagram

2010-05-20 Thread Shi, Tao
Thank you for the suggestions, Peter and David! - Original Message From: Peter Ehlers ehl...@ucalgary.ca To: David Winsemius dwinsem...@comcast.net Cc: Shi, Tao shida...@yahoo.com; r-help@r-project.org Sent: Wed, May 19, 2010 8:45:24 PM Subject: Re: [R] colored venn diagram

Re: [R] Concatenation

2010-05-20 Thread Shi, Tao
Still not clear. Follow the posting guide and some examples always help. Not sure how you concatenate numbers? You can try, for example, apply(df[1:10,], 2, paste, collapse= ) but this will turn everything into strings. Is this what you want? ..Tao - Original Message From

Re: [R] Concatenation

2010-05-21 Thread Shi, Tao
Well, it seems a simple c will do. e.g. df - matrix(rnorm(100), ncol=4) c(t(df[1:10,])) # concatenate rows c(df) # concatenate columns From: santana sarma aimanusa...@gmail.com To: Shi, Tao shida...@yahoo.com; r-help@r-project.org Sent: Thu, May 20, 2010 11:45:25 PM Subject: Re

Re: [R] Fast Matrix Computation

2010-05-22 Thread Shi, Tao
One way to do it: apply(B, 1, function(x) t(apply(A, 1, function(y) abs(y-x) )) ) - Original Message From: David Neu da...@davidneu.com To: r-help@r-project.org Sent: Sat, May 22, 2010 10:35:32 AM Subject: [R] Fast Matrix Computation Hi, I have two (large) matrices A and B

Re: [R] R eat my data

2010-05-25 Thread Shi, Tao
reading. See if any thing strange there.) Also, changing read.table to read.delim often works. ...Tao - Original Message From: Changbin Du changb...@gmail.com To: David Winsemius dwinsem...@comcast.net Cc: r-help@r-project.org Sent: Tue, May 25, 2010 9:12:58 AM Subject: Re

Re: [R] doubt about samr siggenes.table$genes.up

2010-05-25 Thread Shi, Tao
Leo, This question is more suited for the BioC help list. Please post your question over there. Also, offer more info about how you arrived this siggenes.table, what version of R and siggenes package you're using would be helpful. Please read the posting guide! ...Tao - Original

Re: [R] silhouette: clustering labels have to be consecutive intergers starting

2007-10-10 Thread Tao Shi
Thank you very much, Benilton and Prof. Ripley, for the speedy replies! Looking forward to the fix! Tao From: Prof Brian Ripley [EMAIL PROTECTED] To: Benilton Carvalho [EMAIL PROTECTED] CC: Tao Shi [EMAIL PROTECTED], [EMAIL PROTECTED], r-help@r-project.org Subject: Re: [R

[R] RMySQL installation problem

2007-11-15 Thread Tao Shi
was not found. Re-installing the application may fix the problem. I tried the re-installation. It didn't work. The DBI package I have is version 0.2-4, just in case. thanks, ...Tao _ Climb to the top of the charts!  Play Star

[R] SVM cross validation in e1071

2009-07-07 Thread Tao Shi
include lda as a comparison. I'm using WinXP, R-2.9.0, and e1071_1.5-19. Many thanks! ...Tao ##== ## manual ## ## set.seed(1234) dat - data.frame( rbind(matrix(rnorm(1000),ncol = 10), matrix(rnorm(1000, mean=0.6),ncol = 10))) cl - as.factor

Re: [R] SVM cross validation in e1071

2009-07-07 Thread Tao Shi
if the default value for fitted should be set to FALSE to avoid the confusion: i.e. I thought the fitted field in the returning object is the CV results. ...Tao CC: r-help@r-project.org From: mailinglist.honey...@gmail.com To: shi...@hotmail.com Subject: Re: [R] SVM cross validation in e1071 Date

[R] a question regarding package building

2008-06-26 Thread Tao Shi
manually? Thanks, ...Tao _ The i’m Talkathon starts 6/24/08.  For now, give amongst yourselves. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] a question regarding package building

2008-06-26 Thread Tao Shi
Got it! Thank you very much! ...Tao Date: Thu, 26 Jun 2008 20:41:29 +0200 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] CC: r-help@r-project.org Subject: Re: [R] a question regarding package building Tao Shi wrote: Hi List, In Windows

Re: [R] silhouette: clustering labels have to be consecutive integers starting

2009-05-13 Thread Tao Shi
Thanks, Martin! ...Tao From: maech...@stat.math.ethz.ch Date: Wed, 13 May 2009 14:13:48 +0200 To: shi...@hotmail.com CC: rip...@stats.ox.ac.uk; bcarv...@jhsph.edu; r-help@r-project.org; maech...@stat.math.ethz.ch Subject: Re: [R] silhouette: clustering labels have to be consecutive

[R] updating a helper function in a R package

2007-12-06 Thread Tao Shi
kind of stupid but I’m curious to know) I'm using R-2.5.1, on WinXP. Many thanks, …Tao _ 07 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] nlme model specification (revisit)

2008-01-11 Thread Tao Shi
other attached packages: lattice nlme 0.15-11 3.1-86 Just wondering if you guys have any answers for this since I didn't see any for the original posts. Thanks, ...Tao _ GLM_CPC_VideoChat_distantfamily_012008

[R] follow up on [Rd] NAMESPACE methods guidance, please ( http://tolstoy.newcastle.edu.au/R/e4/devel/08/06/1901.html )

2008-10-07 Thread Tao Shi
, ...Tao _ 50F681DAD532637!5295.entry?ocid=TXT_TAGLM_WL_domore_092008 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] why two diff. se in nlsList?

2009-01-26 Thread Tao Shi
Hi list, In the object returned by summary.nlsList, what's the difference between coefficients and parameters? The have the same Estimate, different se (therefore t value), but same p values. R.2.8.0 on winxp with nlme_3.1-89 Thanks, ...Tao

[R] 2 player Nash Equilbrium game

2008-12-28 Thread tao . wang
Dear all, I am looking to compute a 2-player Nash Equibrium. Has anybody looked into similar things before? I would like to implement it in R instead of doing it in Matlab. Thanks. Tao [[alternative HTML version deleted]] __ R-help@r

[R] latex in Hmisc: cell formating

2009-01-22 Thread Tao Shi
), + cellTexCmds=cell.format, numeric.dollar = FALSE) Thanks! ...Tao __ R-help@r-project.org 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

[R] latex in Hmisc: cell formating

2009-01-23 Thread Tao Shi
, the problem persisted. ...Tao Tao Shi hotmail.com writes: ## I'm using R 2.8.0 on WinXP, Hmisc_3.4-3 table1 - matrix(10, 180,7) cell.format - matrix(, ncol=7, nrow=180) cell.format[c(seq(3,180,6),seq(4,180,6)),] - color{red} cell.format[c(seq(5,180,6),seq(6,180,6)),] - color{green} latex

[R] wrap long lines in table using latex in Hmisc

2010-02-26 Thread Tao Shi
Hi list, Is there a way to control long-line wrapping in a table using latex function in Hmisc or any other functions?  It seems I can't find any examples. Thank you very much! ...Tao

Re: [R] text editors

2010-02-26 Thread Tao Shi
If you do everything in Windows, Tinn-R is one of the best and also the one I use.  I also tried WinEdt. It's very good, but it is not free.  If you want a cross-platform editor, Emacs+ESS is the one.  Like others said, the learning curve is steep, but worth it. ...Tao

Re: [R] wrap long lines in table using latex in Hmisc

2010-02-26 Thread Tao Shi
Hi Ista, Thanks!  I missed that. ...Tao From: istaz...@gmail.com Date: Fri, 26 Feb 2010 17:48:32 -0500 Subject: Re: [R] wrap long lines in table using latex in Hmisc To: shi...@hotmail.com CC: r-help@r-project.org Hi Tao, Just set

Re: [R] wrap long lines in table using latex in Hmisc

2010-03-01 Thread Tao Shi
:  http://nepsweb.co.uk/docs/tableTricks.pdf Thanks again. ...Tao Subject: Re: [R] wrap long lines in table using latex in Hmisc From: marc_schwa...@me.com Date: Fri, 26 Feb 2010 18:03:33 -0600 CC: r-help@r-project.org; shi...@hotmail.com To: ch

[R] lineplot.CI in sciplot: option ci.fun can't be changed?

2010-04-02 Thread Tao Shi
hi List and Manuel, I have encounter the following problem with the function lineplot.CI.  I'm running R 2.10.1, sciplot 1.0-7 on Win XP.  It seems like it's a scoping issue, but I couldn't figure it out. Thanks! ...Tao lineplot.CI(x.factor = dose, response = len, data = ToothGrowth

Re: [R] lineplot.CI in sciplot: option ci.fun can't be changed?

2010-04-02 Thread Tao Shi
: ex.fn - function(x, fun = mean, fun2 = function(x) fun(x)+sd(x)) { list(fun=fun(x), fun2=fun2(x)) } data - rnorm(10) ex.fn(data) #works ex.fn(data, fun=median) #works ex.fn(data, fun2=function(x) fun(x)+3) #error with fun(x) not found On Fri, 2010-04-02 at 17:36 +, Tao Shi wrote

[R] a warning message from heatmap.2

2010-04-05 Thread Tao Shi
Hi List, I want to show the heatmap of a correlation matrix using heatmap.2, however always get this warning message (see below) and the column dendrogram is not showing.  It's not really a big deal, but curious how to suppress it and still let R show what I want to show (i.e. a symmetrical

[R] read password-protected files

2011-03-31 Thread Shi, Tao
Hi list, I have a bunch of .csv files that are password-protected. I wonder if there is a way to read them in in R without manually removing the password protection for each file? Thank you very much! ...Tao [[alternative HTML version deleted

Re: [R] How to update R?

2011-03-31 Thread Shi, Tao
! ...Tao From: Li, Yunfei yunfei...@wsu.edu To: r-help@r-project.org Sent: Thu, March 31, 2011 10:16:23 AM Subject: [R] How to update R? Hi, My R version is old, and I would like to update it. How can I update R and keep all the libraries? Thanks Best, Yunfei

Re: [R] How to update R?

2011-03-31 Thread Shi, Tao
Forgot to mention that my way of updating is for Windows only. ...Tao From: Joshua Wiley jwiley.ps...@gmail.com To: Li, Yunfei yunfei...@wsu.edu Cc: r-help@r-project.org Sent: Thu, March 31, 2011 11:05:29 AM Subject: Re: [R] How to update R? Hi Yunfei

Re: [R] read password-protected files

2011-03-31 Thread Shi, Tao
this? Thanks! ...Tao From: Prof Brian Ripley rip...@stats.ox.ac.uk Cc: r-help@r-project.org Sent: Thu, March 31, 2011 11:15:35 AM Subject: Re: [R] read password-protected files On Thu, 31 Mar 2011, Shi, Tao wrote: Hi list, I have a bunch of .csv files

Re: [R] read password-protected files

2011-03-31 Thread Shi, Tao
Thanks, Josh. I use xlsReadWrite routinely. It would be nice if it has a password option. ...Tao - Original Message From: Joshua Wiley jwiley.ps...@gmail.com To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org Sent: Thu, March 31, 2011 11:50:27 AM Subject: Re: [R] read

Re: [R] read password-protected files

2011-03-31 Thread Shi, Tao
Thanks, Henrique! I'll try that. ...Tao - Original Message From: Henrique Dallazuanna www...@gmail.com To: Shi, Tao shida...@yahoo.com Cc: Prof Brian Ripley rip...@stats.ox.ac.uk; r-help@r-project.org Sent: Thu, March 31, 2011 12:00:45 PM Subject: Re: [R] read password-protected

Re: [R] a bug in write.csv?

2011-04-08 Thread Shi, Tao
As pointed out by Ista, please read the help file. write.csv and write.csv2 provide convenience wrappers for writing CSV files. They set sep, dec and qmethod, and col.names to NA if row.names = TRUE and TRUE otherwise. From: John Kane jrkrid...@yahoo.ca

[R] barplot with both color and shading

2012-02-21 Thread Shi, Tao
Hi list, I want to draw a bar plot with color indicating one grouping and different shading on top of the color indicating another grouping.  How should I proceed? Thanks! ...Tao __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] barplot with both color and shading

2012-02-21 Thread Shi, Tao
: barplot(1:10, angle=20, density=c(0,20), col=rep(0:1, each=5)) ...Tao - Original Message - From: R. Michael Weylandt michael.weyla...@gmail.com To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org r-help@r-project.org Sent: Tuesday, February 21, 2012 6:49 PM Subject: Re: [R] barplot

Re: [R] barplot with both color and shading

2012-02-22 Thread Shi, Tao
Thank you both, Michael and Jim, for the answers! - Original Message - From: Michael Bibo michael_b...@health.qld.gov.au To: r-h...@stat.math.ethz.ch Cc: Sent: Wednesday, February 22, 2012 12:46 AM Subject: Re: [R] barplot with both color and shading Shi, Tao shidaxia

Re: [R] UniCox in R

2010-12-17 Thread Shi, Tao
[[1]]-aa[[3]], aa[[2]], aa[[1]]+aa[[3]]) Can't figure out how he made the horizontal bars at the ends of each error bar vary in length. How do you like uniCox? Have seen my previous post regarding the error message? https://stat.ethz.ch/pipermail/r-help/2010-November/260857.html ...Tao

[R] best subset selection on random effects model

2012-02-12 Thread Tao Zhang
Hi, I know leaps() computes the best subset selection for linear model, and the bestglm() computes the best subset selection for generalized linear model. Is there any package for best subset selection on random effects model, or mixed effects model? Thank you so much. Tao

Re: [R] R versions for Red Hat Linux

2011-05-11 Thread Shi, Tao
I second Marc on this. I just installed R2.12.2 on my linux box running RHEL. Now I also see R2.13.0 available by typing yum info R. ...Tao - Original Message From: Marc Schwartz marc_schwa...@me.com To: Marta Avalos marta.ava...@isped.u-bordeaux2.fr Cc: r-help@r-project.org

[R] changes in coxph in survival from older version?

2011-05-11 Thread Shi, Tao
which one I should trust? Thanks, ...Tao # R2.13.0, survival 2.36-9 = dat=read.csv(tmp1.csv, header=T) fit2 - coxph(Surv(tt, cens) ~., data=dat) Warning message: In fitter(X, Y, strats, offset, init, control, weights = weights

Re: [R] changes in coxph in survival from older version?

2011-05-16 Thread Shi, Tao
didn't even converge? Thanks! ...Tao - Original Message From: Terry Therneau thern...@mayo.edu To: Shi, Tao shida...@yahoo.com Cc: r-help@r-project.org Sent: Thu, May 12, 2011 6:42:09 AM Subject: Re: changes in coxph in survival from older version? On Wed, 2011-05-11 at 16:11

Re: [R] changes in coxph in survival from older version?

2011-05-16 Thread Shi, Tao
, I believe I can always narrow down the variables based on univarate analysis and build a multivariate model from that. Many thanks in advance. ...Tao - Original Message From: Frank Harrell f.harr...@vanderbilt.edu To: r-help@r-project.org Sent: Mon, May 16, 2011 11:25:20 AM

  1   2   >