similar to: Merging RData files

?load Read this carefully. Pay attention to its instructions re: overwriting existing objects. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Jan 16, 2018 at 12:43 AM, Steven Yen <styen at ntu.edu.tw> wrote: >

Hi all, This is great. Again, here is what I need. I run two separate jobs (a.R and b.R) with results (say regression outputs) going to a.RData and b.RData. I like to put all results in one place (where I can retrieve them in one place, ab.RData). The following codes do it (I am not sure if line 2 is needed but I am happy). Thank you all. load("a.RData")

Huh. I may by completely wrong but you cannot do such "merging". .RData files are AFAIK places where all objects from given session are stored. However you could load each .RData file and save/export result (one object). BTW, what do you mean exactly by "combine/consolidate"? And finally, post your questions in plain text not html, otherwise they can be mangled. Cheers

On 16/01/2018 6:33 AM, Steven Yen wrote: > Hi all, > This is great. Again, here is what I need. I run two separate jobs (a.R > and b.R) with results (say regression outputs) going to a.RData and > b.RData. I like to put all results in one place (where I can retrieve > them in one place, ab.RData). The following codes do it (I am not sure > if line 2 is needed but I am happy).

Understood. In my case, a.RData and b.RData contain identical variables/data, plus simulation outputs from separate runs. The codes deliver what I need. Good to know the three lines work. Thank you. On 1/16/2018 8:06 PM, Duncan Murdoch wrote: > On 16/01/2018 6:33 AM, Steven Yen wrote: >> Hi all, >> This is great. Again, here is what I need. I run two separate jobs (a.R >>

Hey, I need some help. I want to obtain a cross validation for a regression model (binary response) but I got an error with CVbinary. Well I did this: fit <- lm(resp ~ PC1 + PC2 + PC3 + PC4 + PC5 + PC6 + PC7 + PC8 + PC9+PC10+PC11+PC12+PC13+PC14+PC15+PC16+PC17+PC18+PC19+PC20+PC21+PC22+PC23+PC24+PC25+PC26+PC27+PC28, data = dexp.cp, family=binomial()) CVbinary(fit) Error in sample(nfolds, m,

Dear Martin, Have you try to create a new variable for age squared, say agesq? If you fit the model using this new variable you should get the coefficients. So your new model is something like height~age+agesq I hope this helps, Saghir > -----Original Message----- > From: MARTIN Ludovic [SMTP:martinl@mathinfo.ens.univ-reims.fr] > Sent: Wednesday, 07 May, 2003 2:02 PM > To:

hello everyone, I working on the dataset for my project in class and got stuck on trying to run logistic regression. here is my code: data <- read.csv(file="C:/Users/fieder.data.2000.csv") # creating subset of men fieder.male<-subset(data,data[,8]==1) unmarried.male<-subset(data,data[,8]==1&data[,6]==1) # glm fit agesq.male<-(unmarried.male[,5])^2

Dear all, Does any one knows why should I get the following error message, when trying to do a simple data.frame?? DataF<-data.frame(Subject,BiomR,Spp,Capas,Litter,Herbs,LitterD,MaxCanH,DDifS p,DSSp,Slope, CanDens,NearestSp) Erro em data.frame(Subject, BiomR, Spp, Capas, Litter, Herbs, LitterD, : arguments imply differing number of rows: 202, 0 The data I am using

Dear Forum, I have a question about interaction estimate in the Cox model: why the hazard ratio of the interaction is not produced in the summary of the model? (Instead, the estimate of the coefficient is given in the print of the model.) # Example: modINT <-cph( Surv(T_BASE, T_FIN,STATUS)~ NYHA + ASINI + RFP + FE_REC + XX_PR*XX_DISF) print(modINT) coef se(coef) z

Hello -- I would like to know of a more efficient way of writing the following piece of code. Thanks. options(stringsAsFactors=FALSE) orig <- c(rep('11111111',100000),rep('22222222',200000),rep('33333333',300000),rep('44444444',400000)) orig.unique <- unique(orig) system.time(df <- as.data.frame(sapply(orig.unique, function(x) ifelse(orig==x, 1, 0))))

Hi, I use glht to make multcomp, using Tukey, from a glm model. It is possible to get a more simplified output of result? Somethink like ordering by letters. Thanks Ronaldo -- Human kind cannot bear very much reality. -- T. S. Eliot, "Four Quartets: Burnt Norton" -- > Prof. Ronaldo Reis J?nior | .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia | : :' :

Hi, The fact that apparently the last tigervnc update from upstream was missed triggered me to check for missing updates and packages in 6.3. Here are my results. Sorry for any false positives that might have crept in, but note that some of the 6_x updates actually are updates and not a parsing error. And perhaps an occasional false positive due to having to compare upstream SRPMS vs downstream

Just as a followup even g++ is broken: [piot at paris ~]$ g++ bash: g++: command not found... Install package 'gcc-c++' to provide command 'g++'? [N/y] y * Waiting in queue... The following packages have to be downgraded: cpp-8.2.1-3.5.el8.x86_64 The C Preprocessor glibc-2.28-42.el8.1.i686 The GNU libc libraries glibc-2.28-42.el8.1.x86_64 The GNU libc libraries

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1 t2 t3 ................ 0.00208 0.00417 0.00625 ................. a1 a2 a3 ...................

Does R graphics have a way to easily label the horizontal axis by the reciprocal of the scaled value? -- View this message in context: http://r.789695.n4.nabble.com/Graphics-Reciprocal-labeling-tp3949054p3949054.html Sent from the R help mailing list archive at Nabble.com.

I performed the update yesterday that included installing kernel 2.6.9-42.0.2.EL and replacing ethereal with wireshark. First, the trivial... Wireshark would not run from the KDE menu without editing the menu to change the command from "wireshark" to "kdesu wireshark" because dumpcap is at usr/sbin/dumpcap instead of /usr/bin/dumpcap . Making the change causes wireshark to

CentOS Errata and Enhancement Advisory 2012:0564 Upstream details at : https://rhn.redhat.com/errata/RHEA-2012-0564.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 952425f3d4da6243fc5dbafe35279ddad98e91244b7a13fe0a5fc69c38929e76 kmod-tg3-3.122-1.el6_2.i686.rpm x86_64:

CentOS Errata and Enhancement Advisory 2012:0576 Upstream details at : https://rhn.redhat.com/errata/RHEA-2012-0576.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) i386: 6b7cb779728387373da59d1cc5d6e7d298dbe5f06b0c5fc2efa8d09e90cc3430 kmod-tg3-3.122-2.el6_2.i686.rpm x86_64:

Hi Elisa, I hope this is what you wanted. dat1<-read.csv("peaks.csv",sep=",") #Subset dat2<-dat1[1:5,] res1<-do.call(cbind,lapply(seq_len(nrow(dat2)),function(i) do.call(rbind,lapply(split(rbind(dat2[i,],dat2[-i,]),1:nrow(rbind(dat2[i,],dat2[-i,]))), function(x) {x1<-rbind(dat2[i,],x);