[R] Call for Abstracts: useR! 2021
Hello all, submissions for useR! 2021 are open until March 15th. You can submit an abstract for a regular talk, an elevator pitch (that would be the equivalent of a poster at an in-person conference), a panel, or an incubator session. Full details here: https://user2021.r-project.org/participation/call-for-abstracts/ Abstract submission: Jan 25–Mar 15 Deadline for video, blogs, scripts: Jun 20 Registration opens: Apr 15 (rates depend on location and employment, from 0 to 100 USD, https://user2021.r-project.org/participation/registration/) Conference dates: July 05–09 (virtual) useR! call for abstracts on Twitter: https://twitter.com/_useRconf/status/1354352011978956806?s=09 Please share the announcement with relevant units. Thank you and hope to see you soon! On behalf of the useR! Program Committee, Dorothea Hug Peter -- *useR! 2021 - user2...@r-project.org * *user2021.r-project.org <http://user2021.r-project.org/> - @_useRConf <https://twitter.com/_useRConf> - LinkedIn Profile <https://www.linkedin.com/company/user-conf>* [[alternative HTML version deleted]] ___ r-annou...@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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 code.
[R] The useR! 2021 Call for Tutorials is open now!
The useR! 2021 Call for Tutorials is open now! *(Due Date Friday, February 5, 2021)* useR! 2021 is going virtual, and we are happy to open the call for the tutorial proposals. The tutorials will be for R users in all sectors, e.g., researchers, government officers and industry representatives who focus on the applicability of R in practical settings and for all levels. Each tutorial's duration can vary between two to four hours, we encourage shorter length, but we expect at least two hours, and we will accept proposals for tutorials to be taught in *English, Spanish, or French.* These tutorials will be organized in *multiple time zones*, so we request that you select a convenient time zone. We will later determine the timezone based on all other tutorials and attendees’ interest to combine them into a session. useR! offers *USD 500* for selected tutorials. Find the full call for tutorials here: https://user2021.r-project.org/participation/call-for-tutorials/ *We would appreciate it if you spread this call*,* and we await your proposals.* useR! 2021 Organizing Committee [[alternative HTML version deleted]] ___ r-annou...@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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 code.
[R] how to merge GRange object?
Hello everyone, I am wondering how to simply merge two GRanges objects by range field and add the value by additional vector. For example, I have two objects below obj1 seqnames ranges strand | Val | [1] chr1_random [272531, 272571] + | 88 [2] chr1_random [272871, 272911] + | 45 obj2 seqnames ranges strand | Val | [1] chr1_random [272531, 272581] + | 800 [2] chr1_random [272850, 272911] + | 450 after merged, it should be an object as the following mergedObject and it would concern the differences in IRANGE data (e.g. 581 and 850 in obj2 above were different from those of obj1, which were 571 and 871 respectively) mergedObject seqnames ranges strand | object2Val object1Val | [1] chr1_random [272531, 272581] + | 800 88 [2] chr1_random [272850, 272911] + | 450 45 On Wednesday, October 16, 2013 8:31 AM, Terry Therneau wrote: On 10/16/2013 05:00 AM, r-help-requ...@r-project.org wrote: > Hello, > > I'm trying to use coxph() function to fit a very simple Cox proportional > hazards regression model (only one covariate) but the parameter space is > restricted to an open set (0, 1). Can I still obtain a valid estimate by > using coxph function in this scenario? If yes, how? Any suggestion would be > greatly appreciated. Thanks!!! Easily: 1. Fit the unrestricted model. If the solution is in 0-1 you are done. 2. If it is outside, fix the coefficient. Say that the solution is 1.73, then the optimal solution under contraint is 1. Redo the fit adding the paramters "init=1, iter=0". This forces the program to give the loglik and etc for the fixed coefficient of 1.0. Terry Therneau __ 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 code. [[alternative HTML version deleted]] __ 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 code.
Re: [R] R web application development
Thanks for reply, but which archives? Thanks again. John From: Stephen Sefick Cc: "r-help@r-project.org" Sent: Monday, March 18, 2013 8:02 PM Subject: Re: [R] R web application development I can't offer any advice, but I feel like you could probably get a good start on this by looking through the archives. On Mon 18 Mar 2013 06:55:33 PM CDT, John linux-user wrote: > > Dear all, > > I am wondering if what would be the simple way to develop a simple web > application that runs R. That is, the web application allows any user upload > a dataframe as a variable to my web server, a linux-based apache, and then > run a R package (my package) on the variable that should ideally be handled > as a variable in memory instead of saving to the disk in my server , for > security concern. After running, the result would be returned to the web. Any > suggestion will be appreciated. > > Best, > > John > [[alternative HTML version deleted]] > > > > __ > 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 code. [[alternative HTML version deleted]] __ 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 code.
[R] R web application development
Dear all, I am wondering if what would be the simple way to develop a simple web application that runs R. That is, the web application allows any user upload a dataframe as a variable to my web server, a linux-based apache, and then run a R package (my package) on the variable that should ideally be handled as a variable in memory instead of saving to the disk in my server , for security concern. After running, the result would be returned to the web. Any suggestion will be appreciated. Best, John [[alternative HTML version deleted]] __ 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 code.
[R] Model Description
Hello dear, I am new to R, Have developed the model for prediction. I dont know exactly about the followed terms residual standard error degrees of freedom, multiple R-squared, adjusted R-squared F-statistics p-values Thanks in advance. [[alternative HTML version deleted]] __ 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 code.
Re: [R] add vectors to multiple objects
Maybe what I previously posted was not clear enough or something else. All vectors L1,L2.. and objects (e.g. a.list, b.list, c.list) already exit or easily to be created in a workspace. Do not worry about those. The key question is how to systematically append/assign these vectors to many objects (hundreds) in a loop without explicitly typing these object names one by one. Please refer to my code posted below to detail my question. I really appreciate your ideas/suggestions. John > listObjects=ls(pattern=".list") #for example, a.list, b.list... > > > for (object in listObjects){ > > #how to sign vectors L1,L2,..Ln to each object (a.list, b.list...) > > > } >lf=ls(pattern=".lst") > > for (x in listfiles) { > dat=read.delim(x,header=F) > > > for (i in lf) { > assign(i$add,as.numeric(dat[,3])) > #or i$add=as.numeric(dat[,3] > names(i)[names(i)=="add"]=substr(x,1,5) > > print (i[1:3,]) > }} From: David Winsemius Cc: "r-help@r-project.org" Sent: Friday, August 3, 2012 2:58 PM Subject: Re: [R] add vectors to multiple objects On Aug 3, 2012, at 8:19 AM, John linux-user wrote: > Hi David, > > Thanks for response, but my key question still remains unsolved. That is, how > to add many vectors (L1,L2,L3) to many list objects (a.list, b.list) > in a workspace? > > listObjects=ls(pattern=".list") #for example, a.list, b.list... > > > for (object in listObjects){ > > #how to sign vectors L1,L2,..Ln to each object (a.list, b.list...) > > > } > > Any suggestions/comments/ideas will be appreciated. Can we have code that creates these entities in the workspace? (That is what is meant as a reproducible example.) I cannot tell whether there are already objects named "L1", "L2",... to which you wish to append the value of the objects found, or that you want to assign values to names that only exist as character vectors. --David. > > John > > > From: David Winsemius > Cc: "r-help@r-project.org" > Sent: Wednesday, August 1, 2012 7:21 PM > Subject: Re: [R] add vectors to multiple objects > > > On Aug 1, 2012, at 8:11 AM, John linux-user wrote: > > > Hi everyone, > > > > I try to add many vectors (L1,L2,L3) to many list objects (a.list, > > b.list) in a workspace. Somethings like below, but it is not working. > > Any suggestions will be appreciated. Best, John > > > > lf=ls(pattern=".lst") > > > > for (x in listfiles) { > > dat=read.delim(x,header=F) > > Presumably that would fail since 'listfiles' has not been defined. did you > mean 'lf'? If you did,then wouldn't the second line overwrite all the early > values of "dat" leaving only the last one? > > Perhaps: > datfils <- list() > for (x in listfiles) { > datfils[x] <- read.delim(x,header=F) > > > > > > for (i in 1: lf) { > > And that would fail because 'lf' is a character vector, and it's not > meaningful to specify such a range. Try instead: > > for (i in names(datfils[x]) ) { > # > # which will then iterate over the names of the files which are now also the > names of the list elements > > > assign(i$add,as.numeric(dat[,3])) > > But since 'i' is a length-1 character vector, the expression `i$add` will be > meaningless. The "$" operator does not do function calling in R unless you do > fancy things with environments, and you cannot "sub-assign" in that manner, > at least not with the assign() function. > > Try instead: > > assign(i, as.numeric(datfils[x][,3])) > names(i)[length(i)] <- "add" > > Or: > i <- transform(i, add=datfils[x][,3] ) > > > > > #or i$add=as.numeric(dat[,3] > > names(i)[names(i)=="add"]=substr(x,1,5) > > I'm not sure these would be doing the same thing. What was your goal here? > > > > > print (i[1:3,]) > > } > > } > > > David Winsemius, MD > Alameda, CA, USA > > > David Winsemius, MD Alameda, CA, USA [[alternative HTML version deleted]] __ 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 code.
Re: [R] add vectors to multiple objects
Hi David, Thanks for response, but my key question still remains unsolved. That is, how to add many vectors (L1,L2,L3) to many list objects (a.list, b.list) in a workspace? listObjects=ls(pattern=".list") #for example, a.list, b.list... for (object in listObjects){ #how to sign vectors L1,L2,..Ln to each object (a.list, b.list...) } Any suggestions/comments/ideas will be appreciated. John From: David Winsemius Cc: "r-help@r-project.org" Sent: Wednesday, August 1, 2012 7:21 PM Subject: Re: [R] add vectors to multiple objects On Aug 1, 2012, at 8:11 AM, John linux-user wrote: > Hi everyone, > > I try to add many vectors (L1,L2,L3) to many list objects (a.list, > b.list) in a workspace. Somethings like below, but it is not working. Any > suggestions will be appreciated. Best, John > > lf=ls(pattern=".lst") > > for (x in listfiles) { > dat=read.delim(x,header=F) Presumably that would fail since 'listfiles' has not been defined. did you mean 'lf'? If you did,then wouldn't the second line overwrite all the early values of "dat" leaving only the last one? Perhaps: datfils <- list() for (x in listfiles) { datfils[x] <- read.delim(x,header=F) > > for (i in 1: lf) { And that would fail because 'lf' is a character vector, and it's not meaningful to specify such a range. Try instead: for (i in names(datfils[x]) ) { # # which will then iterate over the names of the files which are now also the names of the list elements > assign(i$add,as.numeric(dat[,3])) But since 'i' is a length-1 character vector, the expression `i$add` will be meaningless. The "$" operator does not do function calling in R unless you do fancy things with environments, and you cannot "sub-assign" in that manner, at least not with the assign() function. Try instead: assign(i, as.numeric(datfils[x][,3])) names(i)[length(i)] <- "add" Or: i <- transform(i, add=datfils[x][,3] ) > #or i$add=as.numeric(dat[,3] > names(i)[names(i)=="add"]=substr(x,1,5) I'm not sure these would be doing the same thing. What was your goal here? > > print (i[1:3,]) > } > } David Winsemius, MD Alameda, CA, USA [[alternative HTML version deleted]] __ 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 code.
[R] assign vectors to objects
Hi everyone, I try to add many vectors (L1,L2,L3) to multiple list objects (a.list, b.list) in a workspace. Somethings like below, but it is not working. Any suggestions will be appreciated. Best, John lf=ls(pattern=".lst") for (x in listfiles) { dat=read.delim(x,header=F) for (i in lf) { assign(i$add,as.numeric(dat[,3])) #or i$add=as.numeric(dat[,3] names(i)[names(i)=="add"]=substr(x,1,5) print (i[1:3,]) }} From: "r-help-boun...@r-project.org" Sent: Wednesday, August 1, 2012 11:11 AM Subject: Your message to R-help awaits moderator approval Your mail to 'R-help' with the subject add vectors to multiple objects Is being held until the list moderator can review it for approval. The reason it is being held: Post to moderated list Either the message will get posted to the list, or you will receive notification of the moderator's decision. If you would like to cancel this posting, please visit the following URL: https://stat.ethz.ch/mailman/confirm/r-help/4bdc02c25f8cb618a623fce61b61fb51189c52ad [[alternative HTML version deleted]] __ 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 code.
[R] add vectors to multiple objects
Hi everyone, I try to add many vectors (L1,L2,L3) to many list objects (a.list, b.list) in a workspace. Somethings like below, but it is not working. Any suggestions will be appreciated. Best, John lf=ls(pattern=".lst") for (x in listfiles) { dat=read.delim(x,header=F) for (i in 1: lf) { assign(i$add,as.numeric(dat[,3])) #or i$add=as.numeric(dat[,3] names(i)[names(i)=="add"]=substr(x,1,5) print (i[1:3,]) } } [[alternative HTML version deleted]] __ 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 code.
[R] Robust estimation of a geometric random variable
Hi, I have a bunch of data which is assumed to be instances of a geometric random variable with outliers. How can I do a robust estimation of the parameter p so that the effect of outliers is minimized? As a part of the estimation process, I also need to know which are the outliers in the data. I found glmrob which does robust estimation of Poisson and binomial random variables but not geometric random variable. I understand that the maximum likelihood estimate of p of geometric random variable is the mean of the instance values. So if we do robust estimate of mean of the instance values, can we say that we are doing robust estimation of the underlying geometric random variable? If so, which method is most suitable for doing the robust mean estimation. (I am a newbie in statistics and R). Thanks suresh [[alternative HTML version deleted]] __ 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 code.
[R] CONFERENCE: useR! 2011, August 16-18, University of Warwick, early bird deadline
This is a reminder that the deadline for abstract submission and early bird registration for the R User Conference 2011 is April 1st. useR! 2011 http://www.R-project.org/useR-2011 A conference centred on the use of R for data analysis and statistical computing. August 16-18, University of Warwick, Coventry, UK CONTRIBUTED TALKS/POSTERS **DEADLINE APRIL 1** http://www.R-project.org/useR-2011/#Call Abstracts are invited on innovative or interesting applications of R. REGISTRATION **EARLY BIRD DEADLINE APRIL 1** http://www.R-project.org/useR-2011/registration Early bird rates are £100/£200/£300 for student/academic/other respectively. PRE-CONFERENCE TUTORIALS, AUGUST 15 http://www.R-project.org/useR-2011/tutorials A range of half-day tutorials presented by R experts. INVITED SPEAKERS http://www.R-project.org/useR-2011/invited Adrian Bowman, Lee Edlefsen, Ulrike Grömping, Wolfgang Huber, Brian Ripley, Jonathan Rougier, Simon Urbanek, and Brandon Whitcher. **New for useR! 2011** All participants are invited to present a Lightning Talk, for which no abstract is required. http://www.R-project.org/useR-2011/lightning These talks provide a 5-minute platform to speak on any R-related topic and should particularly appeal to R newbies. Participants wishing to give such a talk must provide an informative title on their registration form. Further details including information on accommodation and travel can be found on the conference website. Please send enquiries to: user-2...@r-project.org We hope to see you in Coventry! On behalf of the organizing committee: John Aston, Julia Brettschneider, David Firth, Ashley Ford, Ioannis Kosmidis, Tom Nichols, Jennifer Rogers, Elke Thönnes and Heather Turner [[alternative HTML version deleted]] ___ r-annou...@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce __ 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 code.
[R] R2wd
Hi all, I'm currently trying to get accustomed to R and some of its packages, thus need some help for a package called R2wd. As I am trying to use the wdGET() function I'm getting the following message: Error in if (!(tmp[["ActiveDocument"]][["Name"]] == filename)) tmp$Open(paste("path", : argument is of length zero Does anybody know what that means? Best, -- View this message in context: http://r.789695.n4.nabble.com/R2wd-tp3313527p3313527.html Sent from the R help mailing list archive at Nabble.com. __ 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 code.
[R] useR! 2011: abstract submission & registration open
We are happy to inform you that abstract submission and registration for useR! 2011 is now available online, see http://www.R-project.org/useR-2011 This meeting of the R user community will take place at the University of Warwick, Coventry, UK, August 16-18, 2011. The conference schedule comprises invited lectures and user-contributed sessions. In addition half-day tutorials presented by R experts will run on August 15, 2011, prior to the conference. We invite you to submit abstracts for oral or poster presentations on innovative and exciting applications of R. The call for papers along with the link for abstract submission is available at http://www.R-project.org/useR-2011/#Call In addition to the regular contributed talks, all participants are invited to present a Lightning Talk, for which no abstract is required. These talks provide a 5-minute platform to speak on any R-related topic and should particularly appeal to R newbies. Participants wishing to give such a talk must provide an informative title on their registration form. We hope to see you in Coventry! The organizing committee: John Aston, Julia Brettschneider, David Firth, Ashley Ford, Ioannis Kosmidis, Tom Nichols, Elke Thönnes and Heather Turner ___ r-annou...@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce __ 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 code.
[R] svydesign syntax
This message is for those familiar with the survey package. I need to fit a weighted Cox model to accommodate the sampling weights as I have a case-control study with controls sampled at random from a database in a ratio 2:1 to cases (whom were all sampled). I want to make sure I am using the right svydesign syntax to specify this sampling design. Can anyone please check if the statement below is appropriate for my design? #group represents the case (total of 132) vs control (253 out of the total of 853 controls) groups; prob is 1 for cases and 253/853 for controls and ssize=132 for cases and 853 otherwise; dstr=svydesign(id=~1, strata=~group, prob=~prob, fpc=~ssize, data=noNA) -- View this message in context: http://r.789695.n4.nabble.com/svydesign-syntax-tp2298793p2298793.html Sent from the R help mailing list archive at Nabble.com. __ 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 code.
Re: [R] calibration and validation for svycoxph
Hi Dr. Lumley, Sorry about the succession of the 2 messages. It wasn't intended. This is my first post to the forum and I first signed up to Nabble and submitted the post and later got the message from Nabble that the post has not been accepted by the R forum and that I need separate subscription for the R forum. So I signed up for the R forum and resubmitted the post and both messages ended up being posted. Yes, I was looking for similar functions to those in the Design package to build a calibration plot (Design, calibrate function) and to do bootstrap validation of the model (Design, validate function). I will check the options of the functions you mentioned to see if I can easily get those. Could you please advise on what measure should I use to quantify the optimism of fit (the way c-index was used for Cox models)? If you have existent code that you used for these purposes that would be great. Many thanks. -- View this message in context: http://r.789695.n4.nabble.com/calibration-and-validation-for-svycoxph-tp2249118p2250730.html Sent from the R help mailing list archive at Nabble.com. __ 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 code.
[R] calibration and validation for svycoxph
Hello, This post is for Dr. Thomas Lumley or anybody familiar with the survey package. I am estimating a proportional hazards model with weights using svycoxph. Are there functions already built in the survey package that allow me to do validation and calibration of the model? Thanks -- View this message in context: http://r.789695.n4.nabble.com/calibration-and-validation-for-svycoxph-tp2249118p2249118.html Sent from the R help mailing list archive at Nabble.com. __ 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 code.
[R] Calibration and validation for svycoxph
Hello, This post is for Dr. Thomas Lumley or anybody familiar with the survey package. I am estimating a proportional hazards model with weights using svycoxph. Are there functions already built in the survey package that allow me to do validation and calibration of the model? Thanks -- View this message in context: http://r.789695.n4.nabble.com/Calibration-and-validation-for-svycoxph-tp2249004p2249004.html Sent from the R help mailing list archive at Nabble.com. __ 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 code.
[R] complex--scatter plot
My input "ID" "Label" "*Stype*" "Ntype" "Stype_No""*log*" "S1" "xxx" "A/A" 1 6 2.8 "S1" "xxx" "A/G" 2 2 3 "S1" "xxx" "G/G" 3 1 4 "S2" "yyy" "A/A" 1 1 6.8 "S2" "yyy" "A/G" 2 2 7 "S2" "yyy" "G/G" 3 6 7.4 "S2" "yyy" "NULL""null""null"8 "S3" "zzz" "A/A" 1 3 12 "S3" "zzz" "A/G" 2 3 14 "S3" "zzz" "G/G" 3 3 16 "S3" "zzz" "NULL""null""null"18 "S3" "zzz" "NULL""null""null"20 I would like to draw a single scatter plot for every ID (S1 -plot and S2-Plot)by using *Stype* as X axis and *log* as Y-axis with a header "D_Scatterplot". For example S1-Plot looks like this X-axis consider Stype A/A=1 A/G=2 G/G=3 and their Stype_No (6,2,1)means A/A i.e, 1 repeated 6 times with log values 2.8,3 and 4 and A/G i.e. 2 repeated 2 times and G/G i.e. 3 1 time with the same log values. simply it looks like this 1 2 3 2.8 1 2 3 1 4 1 1 1 Y-axis Label_ 6 | log 5 | 4 | *** 3 | ** ** 2 | ** ** * 1 | 0 |__ A/A A/G G/G ID* Thanx n advance Pearl __ 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 code.
[R] One rather theoretical question about fitting algorithm
Hi, What I'm trying to achieve is very fast algorithm for fitting logistic regression model. I have to estimate regression coeficients using about 10k observations. Once I have coefficients estimated, new 100 rows of data becomes available Now I need to reestimate coeficients using 100 newly arrived observations and removing 100 oldest observations. So, my question is would it be possible to somehow reuse pre iusly calculated coeficients and only adjust them cor newly arrived data? I know it would have to be some aproximation but I suppose it will be good enough. I dont mind doing this in straight C because of of speed perative. Actualy this will have to be cAlculated in a fraction of second. Any ideas would be higly appreciated __ 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 code.
[R] 3d cloud plot with point size reflecting variable value
Hello all, I'm using the cloud function to plot the iris data as per the document: http://www.stat.ucl.ac.be/ISpersonnel/lecoutre/stats/fichiers/_gallery.pdf I'd like to change the point size to reflect a fourth variable, as done here http://www.cis.hut.fi/projects/somtoolbox/download/pics2/shotvs2_origdata.png Does anybody know how to do this? There doesn't seem to be an option in cloud() Thanks, Richie [[alternative HTML version deleted]] __ 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 code.
Re: [R] including tabular information with a plot in R Graphics
Hi all, Thanks very much for the suggestions. My experience for the benefit of others is: - use addtable2plot if you want a fancier legend or small table *within* a plot area - use textplot if you want combinations of tables and text listings in the graphics device. The help for textplot gives a very good example Thanks again, Richie 2009/3/5 Eik Vettorazzi > Hi Richie, > have a look at textplot() from the gplots-package. > hope, it helps! > > R User R User schrieb: > >> Hi all, >> I have a presentation problem that I cannot find a solution to in the >> documetnation. >> I have a nice barplot. Below this I would also like a table with some >> information relating to the plot. >> My idea was to have the plot and table on the same graphics window so I >> can >> output them as an image file for a report. >> >> Does anybody know how to include tabular information from a dataset or >> table >> into the graphics device? >> Is there a better/more correct way to output tables and graphics to a >> single >> file? >> >> Thanks very much, >> Richie >> >>[[alternative HTML version deleted]] >> >> __ >> 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 code. >> >> > > -- > Eik Vettorazzi > Institut für Medizinische Biometrie und Epidemiologie > Universitätsklinikum Hamburg-Eppendorf > > Martinistr. 52 > 20246 Hamburg > > T ++49/40/42803-8243 > F ++49/40/42803-7790 > > [[alternative HTML version deleted]] __ 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 code.
[R] including tabular information with a plot in R Graphics
Hi all, I have a presentation problem that I cannot find a solution to in the documetnation. I have a nice barplot. Below this I would also like a table with some information relating to the plot. My idea was to have the plot and table on the same graphics window so I can output them as an image file for a report. Does anybody know how to include tabular information from a dataset or table into the graphics device? Is there a better/more correct way to output tables and graphics to a single file? Thanks very much, Richie [[alternative HTML version deleted]] __ 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 code.
[R] barplot with specific order of x axis labels
Hi all, I'd appreciate your help with this problem. I need to plot a barplot with the categories in a specific order. My data might be: hours Freq AN 10 MO 14 LU 30 I need the categories to be in the order: MO LU AN Is there some way to pass sort the dataset into this order so that it plots correctly when passed to barplot()? Thanks very much, Richie [[alternative HTML version deleted]] __ 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 code.
[R] vertically aligned X axis labels disappear off R Graphics window
Hi guys, I'm evaluating R for basic data exploration. I produce a bar plot of the data, with the x axis labels aligned vertically. However, the start of labels longer than about 10 characters are cut off by the bottom of the graphics window. I'd appreciate your help in properly spacing the space around a plot reserved for labels so that this problem can be avoided. Thanks, Richie [[alternative HTML version deleted]] __ 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 code.
[R] multiple levels of nesting with trellis plots
Hi guys, I have a tricky problem that I'd appreciate your help with. I have two categorical variables, say varA and varB and an associated frequency Freq for combinations of the levels of varA and varB. This was created with a table() call. I'd now like to make panel plots of the frequency. I can do a panel plot by varA: barchart(data$Freq ~ data$varB | data$varA) How do I achieve two levels of nesting for say a third categorical variable varC? That is, how do I specify 2 'given's in the the panel plot, and preferably arranged the panels vertically within the 'inner' panel/ thanks very much, Richie [[alternative HTML version deleted]] __ 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 code.
[R] how to obtain p values from an ANOVA result
Hi all, I'm tryint to turn my ANOVA result from anova(model) in to p values. The default output is: term ; Df ; Deviance ; resid Df; Resid. Dev Is there a way to add p values to this output so the significance of the terms can be easily determined? Thanks in advance. Richie [[alternative HTML version deleted]] __ 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 code.
[R] barplot() x axes are not updated after removal of categories from the dataframe
Hi all, I'd be grateful for your help. I am a new user struggling with a barplot issue. I am plotting categories (X axis) and their mean count (Y axies) with barplot(). The first call to barplot works fine. I remove records from the dataframe using final=[!final$varname == "some value",] I echo the dataframe and the records are no longer in the dataframe. When I call plot again however, the X axis still contains the removed category values albeit with a count of 0. How can I stop the removed categories appearing in the barplot? Are they stored globally after the first barplot call? thanks in advance for your help. R [[alternative HTML version deleted]] __ 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 code.
[R] Error in storage when attempting to run GLM
Hello, I'd be very grateful for help with the following. I attempt to run a basic GLM as follows: basic_model=glm(factor(A$b) ~ factor(A$c) ) variables b and c are text values and their class is factor. they have no missing values. I get the error: Error: NA/NaN/Inf in foreign function call (arg4) In addition: warning messages: 1: Ops.factor(y,mu) : - not meaningful for factors 2: Ops.factor(eta, offset) : - not meaningful for factors etc What does this mean? R [[alternative HTML version deleted]] __ 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 code.
[R] Help on computing Geary's C statistic to test for Spatial Autocorrelation
Dear Users: I have been trying to use the geary.test() function in *R*, but am having slight difficulty understanding how I am to apply it in my context. I have 2 matrices: 1) *n x p* matrix of *n* observations with *p* measurements each. It may be noted that this matrix has a spatial dimension to it, as the *n*observations are at different geographical locations on a map. 2) *n x n* spatial weights matrix (symmetric) with each matrix element being an inter-locality weight. Using the above 2 matrices, I am trying to compute the Geary's *C* statistic for each of the *p* dimensions and also obtain corresponding *p*-values to determine its statistical significance for spatial autocorrelation. I will appreciate any help you could provide me in helping me compute the same. Thanks. [[alternative HTML version deleted]] __ 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 code.
[R] [R-pkgs] np 0.30-1 (nonparametric kernel smoothing methods for mixed data types) is available on CRAN...
Dear R users, Version 0.30-1 of the np package has been released and uploaded to CRAN. The np package provides nonparametric kernel smoothing methods for mixed data types. We encourage anyone using the package to upgrade to the latest version. Description: This package provides a variety of nonparametric (and semiparametric) kernel methods that seamlessly handle a mix of continuous, unordered, and ordered factor data types. We would like to gratefully acknowledge support from the Natural Sciences and Engineering Research Council of Canada (NSERC:www.nserc.ca), the Social Sciences and Humanities Research Council of Canada (SSHRC:www.sshrc.ca), and the Shared Hierarchical Academic Research Computing Network (SHARCNET:www.sharcnet.ca). License: GPL Note that version 0.30-0 and 0.30-1 provide some much needed functionality and bug fixes to the package. Changes from Version 0.30-0 to 0.30-1 [29-Jan-2009] * predict now supports bandwidth, density, distribution, conbandwidth, condensity, and condistribution objects * Consistently allow predictions for categorical values outside of support of training data Note that predictions based upon unconditional density objects defined over categorical variables that lie outside the support of the training data may no longer be true probabilities (i.e., as defined over the training data and the extended/augmented support -- their sum may exceed one) and may therefore require renormalization by the user * Fixed a numerical issue which could hinder npregbw()'s cross validation with higher-order kernels * Default nmulti in npplregbw() is now set correctly * Fixed a bug with the ridging routine in npscoefbw(), added ridging to npscoef * Fixed trivial i/o issue with "Multistart 1 of" using npscoefbw() Changes from Version 0.20-4 to 0.30-0 [15-Jan-2009] * Added basic user-interrupt checking for all underlying C code so that either (Rterm) or the `STOP' icon (Rgui) will interrupt all running processes. This has a number of desirable side effects in addition to being able to interrupt C-based processes including i) R no longer showing up as `not responding' under the task manager (Windows) or the activity monitor (Mac OS X) and ii) buffered output now being correctly displayed when using Rgui under Windows and Mac OS X Note that repeated interruption of large jobs can reduce available memory under R - if this becomes an issue (i.e., you get a `cannot allocate...' error under R) simply restart R (i.e., exit then run a fresh R session) * Added a function npseed() that allows the user to set/reset the random seed for all underlying C routines * Fixed a bug that caused npplregbw() to ignore any kernel options for the regression of y on z * Refined certain constants used in the normal-reference density bandwidth rule for increased accuracy * Moved from using the maximum likelihood estimate of variance throughout to the degrees of freedom corrected estimate (all variance estimates now change by the factor (n-1)/n) Feedback, comments, bug reports, and suggestions for improvement are always welcome. Best regards, Jeff Racine & Tristen Hayfield -- Professor J. S. Racine Phone: (905) 525 9140 x 23825 Department of EconomicsFAX:(905) 521-8232 McMaster Universitye-mail: raci...@mcmaster.ca 1280 Main St. W.,Hamilton, URL: http://www.economics.mcmaster.ca/racine/ Ontario, Canada. L8S 4M4 `The generation of random numbers is too important to be left to chance.' ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages __ 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 code.
[R] How to see source code for na.omit?
Hi R helpers, I'd like to see source code for some of built-in R functions... for example, I would like too see how "na.omit" was implemented? Thanks? __ 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 code.
Re: [R] Variable shortlisting for the logistic regression
Let's try to bring this discussion back again after Frank made very funny remark! What I'm doing at the moment is: 1. I split dataset in two (development and holdout) 2. I fit single predictor logistic model for every variable and collect following stats: DMaxDeriv=modelD$stats[2] DModelLR=modelD$stats[3] DP=modelD$stats[5] DC=modelD$stats[6] DDxy=modelD$stats[7] DGamma=modelD$stats[8] DTau=modelD$stats[9] DR2=modelD$stats[10] DBier=modelD$stats[11] HMaxDeriv=modelH$stats[2] HModelLR=modelH$stats[3] HP=modelH$stats[5] HC=modelH$stats[6] HDxy=modelH$stats[7] HGamma=modelH$stats[8] HTau=modelH$stats[9] HR2=modelH$stats[10] HBier=modelH$stats[11] where D is prefix for stats on development sample and H is prefix for stats derived from hold out sample 3. Now I screen factor with sommers d grather than 0.3 and relative change on hold out sample is smaller than 5% Any comments are very welcomed On Oct 14, 2:48 pm, John Kane <[EMAIL PROTECTED]> wrote: > --- On Mon, 10/13/08, David Scott <[EMAIL PROTECTED]> wrote: > > > > > From: David Scott <[EMAIL PROTECTED]> > > Subject: Re: [R] Variable shortlisting for the logistic regression > > To: "Frank E Harrell Jr" <[EMAIL PROTECTED]> > > Cc: [EMAIL PROTECTED] > > Received: Monday, October 13, 2008, 6:32 PM > > On Mon, 13 Oct 2008, Frank E Harrell Jr wrote: > > > > useR wrote: > > >> Hi R helpers, > > > >> One rather statistical question? > > > >> What would be the best startegy to shortlist > > thousands of continous > > >> variables automaticaly using R > > >> as the preparation for logistic regression > > modleing! > > > >> Thanks > > > > The easiest approach is to use a random number > > generator. > > > Frank > > > Got a laugh from me Frank! > > > Can I nominate it for a fortune? > > > David > > Seconded. > > __ > [[elided Yahoo spam]] > > __ > [EMAIL PROTECTED] mailing listhttps://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guidehttp://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. __ 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 code.
[R] Variable shortlisting for the logistic regression
Hi R helpers, One rather statistical question? What would be the best startegy to shortlist thousands of continous variables automaticaly using R as the preparation for logistic regression modleing! Thanks __ 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 code.
[R] Stepwise lrm()
Hello, I have the data set of 1 + 49 variables. One of them is binary other are continous. I would like to be able to fit the model with all 49 variables and then run stepwise model selction. I'd appriciate some code snippets... __ 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 code.
[R] TTR getYahooData?
Would it be possible to get data from uk.finance.yahoo.com instead from finance.yahoo.com. Thanks __ 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 code.
[R] Singular information matrix in lrm.fit
Hi R helpers, I'm fitting large number of single factor logistic regression models as a way to immediatly discard factor which are insignificant. Everything works fine expect that for some factors I get error message "Singular information matrix in lrm.fit" which breaks whole execution loop... how to make LRM not to throw this error and simply skip factors with singularity problem... alternative solution also appreciated. Thanks __ 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 code.