similar to: skip the error to continue the logistic regression in a loop

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 plot1 plot2 plot3 plot4 plot5 plot6 plot7 I want them in this order plot1 plot2 plot3 plot4 plot5 plot6 plot7 Could somebody tell me how to do this, please? Thanks so many. Suyan

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 1 12 24 NA 2 20 17 NA 3 NA 19 33 4 NA 20 16 I used setdiff and rbind together and got what I

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

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 this in R? Thanks so many. Suyan

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

Is there any implementation in R for finding the phase shift between two continuous signals. I would like to find the average phase shift for tow signals over two years. thanks Stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted

Hi there I want to pass arguments (i.e. the response variable and the subset argument) in a self-made function to glm. Here is one way I can do this: f.myglm <- function(y,subfact,subval) { glm(d.mydata[,y]~d.mydata[,'x1'],family=binomial,subset=d.mydata[,subfact]==subval) } > str(d.mydata) `data.frame': 15806 obs. of 3 variables: $ y : Factor w/ 2 levels

Suppose we have some glm object such as: myglm <- glm( y ~ x, data=DAT) Is there an elegant way--or the "right way" within the R way of thinking--to obtain the names of the response variable, the predictor variables, and the dataset, as character strings? For instance, suppose the "right way" was to use the (currently fictitious) functions theresponse(), thepredictors(),

Dear all, sorry to bother you all with this but I've been trying to use the loglm in MASS package (v2.8.0) and cannot get any sensible output. I'm wondering am I doing something very foolish or missing something obvious. For example, I tried the documentation help(loglm) example - here's the code # Case 1: frequencies specified as an array. sapply(minn38,

Hello,   I am running a loop to compare some residual deviances obtained from glm, with codes: ........   OUT<-NULL for (i in 1:10){  myglm<-glm(mat ~X1+X2+X3,family = binomial, data =myDATA) OUT<-c(OUT,myglm$deviance) }   .......   In the loop, X1, X2, and X3 chage with i. If X1, X2 and X3 are not highly correlated, OUT is a vector of 10 values.   The problem is:  if there there is

This post is NOT a question, but an answer. For readers please disregard all earlier posts by myself about this question. I'm posting for two reasons. First to say thanks, especially to Dimitris, for suggesting the use of errorest in the ipred library. Second, so that the solution to this problem is in the archives in case it gets asked again. If one wants to run a k-fold cross-validation

Is there a fast way to determine how many documents a term occurs in, besides iterating through every document with TermDocEnum? -- Best regards, Stian Gryt?yr

Hi I have a little R problem. I have created a GLM model in R and now I want to predict some values outside the values I have in the model (extrapolate). I have this code: fitted.model4 <- glm(Yval ~ time, family=gaussian, data=Fuel) The question is - How do I predict a value of Yval ie with a value of time = 340 and also get confidence/prediction intervals for Yvar? I have tried the

Why does the assignment of a 3178x93 object to another 3178x93 object remove the dimension attribute? > GT <- array(dim = c(6,nrow(InData),ncol(InSNPs))) > dim(GT) [1] 6 3178 93 > SNP1 <- InSNPs[InData[,"C1"],] > dim(SNP1) [1] 3178 93 > SNP2 <- InSNPs[InData[,"C2"],] > dim(SNP2) [1] 3178 93 > dim(pmin(SNP1,SNP2)) [1] 3178 93

I have a new version of Rake that I''m (almost) ready to release. But before I make an official release, I would like some Rails folk to give it it try ... Why? Because it contains a workaround for the notorious Ruby 1.8.2 test/unit bug. Using this version of rails, it should be possible to unit test a rails app without patching the Ruby 1.8.2 test/unit library. Where to get it?

http://bugzilla.mindrot.org/show_bug.cgi?id=998 Summary: --without-opensc is not respected Product: Portable OpenSSH Version: -current Platform: All OS/Version: Linux Status: NEW Severity: minor Priority: P2 Component: Build system AssignedTo: openssh-bugs at mindrot.org ReportedBy: stian

http://bugzilla.mindrot.org/show_bug.cgi?id=1157 Summary: ssh-keygen doesn't handle DOS line breaks Product: Portable OpenSSH Version: 3.8.1p1 Platform: All URL: http://openssh.org/txt/draft-ietf-secsh-publickeyfile- 02.txt OS/Version: All Status: NEW Severity: normal

Hi, I have some trouble trying to test the quota plugin. First off all, it says in the wiki: "Dovecot 1.0 supports quota using a quota plugin. It's currently discussed in Dovecot LDA page and it supports only very simple quota calculation which is suitable only for mboxes. Maildir++ quota support will come later." Does that mean that it won't work at all with Maildir++, or

I am pleased to announce the newest release of Locomotive! http://locomotive.sourceforge.net Mostly a bug fix release. If you are using Locomotive on 10.4 you should definitely upgrade. New in 0.2.4 ============ * Rename application entries in status window * Run/Restart/Stop buttons now immediately update * Bundle choosing is fixed for Mac OS X 10.4 * Missing gem dependencies in Rails 0.13.1

Hello, I guess a few of you know about the iptables-p2p project at the http://sourceforge.net/projects/iptables-p2p site. This suite has an excellent filter of which I use today and its running very smoothly. The main reason why I use this is because I desperately need to block torrents. IPP2P does not do this at this time (it seems to me in the doc at least). Is it maybe some way we could merge