Re: [R] pROC

2018-07-24 Thread David Winsemius
> On Jul 23, 2018, at 4:43 AM, Maria Enzian wrote: > > > > Hello, > I'm using the package pROC in RStudio to create my ROC-curves and I have > patients in my data - healthy or sick - in the column "Status" and the value > "SUVmax" to examine it. > I used the following code: > plot.roc(daten

[R] pROC

2018-07-23 Thread Maria Enzian
Hello, I'm using the package pROC in RStudio to create my ROC-curves and I have patients in my data - healthy or sick - in the column "Status" and the value "SUVmax" to examine it. I used the following code: plot.roc(daten$Status,daten$SUVmax,percent=TRUE,ci=TRUE,print.auc=TRUE,main="ROC-Kurve

Re: [R] PROC MIXED RANDOM equivalence in R nlme

2017-08-11 Thread Dennis F. Kahlbaum
Hi David: Thanks for your response and suggestions. --Dennis On 8/11/17 2:03 PM, David Winsemius wrote: On Aug 11, 2017, at 7:56 AM, Dennis F. Kahlbaum wrote: Hi Robert: Thanks for your response, as well. I'm sorry. As you've discovered, I made some capitalization errors when posting my R

Re: [R] PROC MIXED RANDOM equivalence in R nlme

2017-08-11 Thread David Winsemius
> On Aug 11, 2017, at 7:56 AM, Dennis F. Kahlbaum wrote: > > Hi Robert: > > Thanks for your response, as well. > > I'm sorry. As you've discovered, I made some capitalization errors when > posting my R code. The actual R code does use "lme", "data", and "random". > > The dataframe is indeed

Re: [R] PROC MIXED RANDOM equivalence in R nlme

2017-08-11 Thread Dennis F. Kahlbaum
Hi Robert: Thanks for your response, as well. I'm sorry. As you've discovered, I made some capitalization errors when posting my R code. The actual R code does use "lme", "data", and "random". The dataframe is indeed named "emiss" and each item in the formula is a column in the dataframe. I

Re: [R] PROC MIXED RANDOM equivalence in R nlme

2017-08-11 Thread Dennis F. Kahlbaum
Hi Thierry: Thank you for your response. I have been trying to resolve this issue for some time, without success. That's why I've posted the problem hoping that someone with both SAS and R experience would be able to help. As background, this SAS code has been used on many different data sets

Re: [R] PROC MIXED RANDOM equivalence in R nlme

2017-08-11 Thread Robert Baer
On 8/10/2017 8:34 AM, Dennis F. Kahlbaum wrote: -- snip -- I don't have real help, but I'll remind you that R is case sensitive, and it looks like that will be at least one problem in the solution your are working on below: lme not LME data not DATA random = RANDOM -

Re: [R] PROC MIXED RANDOM equivalence in R nlme

2017-08-11 Thread Thierry Onkelinx
Dear Dennis, Your question assumes that people know both SAS PROC MIXED and R nlme. Only a limited number of people do. Add the mathematical formulation of the model. That will increase the number of people that can help you. Adding the number of levels in each categorical variable and the number

[R] PROC MIXED RANDOM equivalence in R nlme

2017-08-10 Thread Dennis F. Kahlbaum
I am trying to reproduce some old SAS PROC MIXED code using R and nlme. The data consists of emission readings from vehicles and fuel properties. All variables are real numbers except "study" and "vehicle", which are character. Unfortunately, since the data are confidential, I'm unable to provi

Re: [R] Proc Mixed variance of random effects in R

2015-06-18 Thread Thierry Onkelinx
Dear Gram, A few things first: Please don't post in HTML, it mangles your text. R-sig-mixed model is a better list for questions on mixed models. Send further replies only to that list and not to r-help. You are probably not fitting the same model in R as the one in SAS. Please provide the equati

[R] Proc Mixed variance of random effects in R

2015-06-17 Thread Grams Robins
Hi, I'm trying to convert the following SAS code in R to get the same result that I get from SAS. Here is the SAS code: DATA plants; INPUT sample $ treatmt $ y ; cards; 1 trt16.426264755 1 trt16.95419631 1 trt16.64385619 1 trt27.34872815

[R] How to deal with multiple class ROC analysis in R (pROC package)?

2013-12-11 Thread Yin Zhao
Here is the question link: http://stackoverflow.com/questions/20507108/how-to-deal-with-multiple-class-roc-analysis-in-r-proc-package Thanks! -- *Yin Zhao <http://www.linkedin.com/in/zhaoyin>* *My Blogger <http://ameenzhao.blogspot.com/>* [[alternative HTML ver

Re: [R] pROC plot.roc - Plotting more than two curves - Is it possible

2013-12-04 Thread Murat Tasan
have you looked at the ROCR package? there are at least a few ways to plot multiple ROC curves with the ROCR functions that come to mind. e.g. if you pass lists of 'scores' and 'labels' in to the prediction(...) method, the default plot(...) method will display all of the ROC (or precision-recall,

[R] pROC plot.roc - Plotting more than two curves - Is it possible

2013-12-04 Thread Lopez, Dan
Hi R Experts, The info in this link suggests that plot in the pROC package can only compare two ROC curves at one time. https://stat.ethz.ch/pipermail/r-packages/2011/001220.html But this is from a couple of years ago. Curious if that has changed and if anyone knows how to plot more than one RO

Re: [R] pROC in "R"

2013-01-30 Thread John Kane
ect.org > Subject: [R] pROC in "R" > > > > B Dear, > > I would like to use "pROC" software for my study, butB I could not > uploaded it in "R". Could you please help me to overcome this problem? > This is the message when I write &q

Re: [R] pROC in "R"

2013-01-28 Thread Ista Zahn
And what exactly is the problem? Your code produced no errors (or if it did you have not shown them to us...) Best, Ista On Mon, Jan 28, 2013 at 9:44 AM, Fethi BEZOUBIRI wrote: > > > Dear, > > I would like to use "pROC" software for my study, but I could not > uploaded it in "R". Could you ple

[R] pROC in "R"

2013-01-28 Thread Fethi BEZOUBIRI
 Dear, I would like to use "pROC" software for my study, but  I could not uploaded it in "R". Could you please help me to overcome this problem? This is the message when I write "Library (pROC)" : Le chargement a nécessité le package : plyr Type 'citation("pROC")' for a citation. Attacheme

Re: [R] pROC and ROCR give different values for AUC

2012-12-19 Thread Thomas Lumley
On Thu, Dec 20, 2012 at 6:36 AM, Ivana Cace wrote: > I considered that but on the list more people are likely to see it. So if > they ran into the same thing they may already have figured out what is > going on and have answers. > Or if not, I'll make a reproducible example, ask the maintainers a

Re: [R] pROC and ROCR give different values for AUC

2012-12-19 Thread Steve Lianoglou
On Wed, Dec 19, 2012 at 12:36 PM, Ivana Cace wrote: > I considered that but on the list more people are likely to see it. So if > they ran into the same thing they may already have figured out what is going > on and have answers. > Or if not, I'll make a reproducible example, ask the maintainers

Re: [R] pROC and ROCR give different values for AUC

2012-12-19 Thread Ivana Cace
I considered that but on the list more people are likely to see it. So if they ran into the same thing they may already have figured out what is going on and have answers. Or if not, I'll make a reproducible example, ask the maintainers and post back here where others may find it. (but not til

Re: [R] pROC and ROCR give different values for AUC

2012-12-19 Thread David Winsemius
On Dec 19, 2012, at 4:47 AM, Ivana Cace wrote: > Packages pROC and ROCR both calculate/approximate the Area Under (Receiver > Operator) Curve. However the results are different. > > I am computing a new variable as a predictor for a label. The new variable is > a (non-linear) function of a set

Re: [R] pROC and ROCR give different values for AUC

2012-12-19 Thread Max Kuhn
A reproducible example sent to the package maintainer(s) might yield results. Max On Wed, Dec 19, 2012 at 7:47 AM, Ivana Cace wrote: > Packages pROC and ROCR both calculate/approximate the Area Under (Receiver > Operator) Curve. However the results are different. > > I am computing a new varia

[R] pROC and ROCR give different values for AUC

2012-12-19 Thread Ivana Cace
Packages pROC and ROCR both calculate/approximate the Area Under (Receiver Operator) Curve. However the results are different. I am computing a new variable as a predictor for a label. The new variable is a (non-linear) function of a set of input values, and I'm checking how different parameter

Re: [R] Proc Nnpar1way with D option - equivalent in R

2012-11-13 Thread SASandRlearn
Dan, what you suggested worked out well. This code below also worked out well for me and it matches with SAS output. Ks <- cbind(x,fitted(d1logit)) ks.df <- data.frame(Ks) x <- subset(ks.df,x==0,select=c(V2)) y <- subset(ks.df,x==1,select=c(V2)) ks.test(x[,'V2'], y[,'V2'], alternative = c("two.si

Re: [R] Proc Nnpar1way with D option - equivalent in R

2012-11-13 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of SASandRlearn > Sent: Tuesday, November 13, 2012 4:26 AM > To: r-help@r-project.org > Subject: Re: [R] Proc Nnpar1way with D option - equivalent in R > >

Re: [R] Proc Nnpar1way with D option - equivalent in R

2012-11-13 Thread SASandRlearn
Dan Thank you for your reply. I will try what you recommended. yes.. i have a 1 and 0 as binary. Here is what i have so far d <- read.csv(c:/test.csv", header=T) dlogit <- glm(x ~ a + b + c, data = d, family = "binomial") attach(d) ks.test(x, fitted(values),alternative = c("two.sided"),exact

Re: [R] Proc Nnpar1way with D option - equivalent in R

2012-11-12 Thread Daniel Nordlund
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of SASandRlearn > Sent: Monday, November 12, 2012 5:55 PM > To: r-help@r-project.org > Subject: [R] Proc Nnpar1way with D option - equivalent in R > >

[R] Proc Nnpar1way with D option - equivalent in R

2012-11-12 Thread SASandRlearn
I am trying to match SAS output with R. I am using Proc Npar1way with D option to get KS test statistic. • Here X is a binary dependent variable and Y is the predicted probabilities; proc npar1way data = mydata D; class x; var y; run; When i try this in R ks.test(x, fitted(y),alternative = c

Re: [R] Proc AutoReg (SAS like Output in R)

2012-05-18 Thread R. Michael Weylandt
Well, the easiest way would be to use SAS ;-) In R take a look at the ar() function. You can view documentation by typing ?ar at the console and see some worked examples by running example(ar) Hope this helps, Michael On Fri, May 18, 2012 at 12:11 AM, anil wrote: > Hi, > > I want to find out, h

[R] Proc AutoReg (SAS like Output in R)

2012-05-17 Thread anil
Hi, I want to find out, how can i get a SAS like output for Proc Autoreg with AR(1) and AR(2) terms. Thanks... -- View this message in context: http://r.789695.n4.nabble.com/Proc-AutoReg-SAS-like-Output-in-R-tp4630468.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] [pROC] roc.test returns "NA" p-value...

2012-04-25 Thread O
Dear Uwe: Many thanks for your prompt feedback. You were right about the data. I don't understand why but one map wasn't perfect (i.e., it didn't include count information for each value). If this happens the map seems alright using GIS programs because you can still produce pretty maps based o

Re: [R] [pROC] roc.test returns "NA" p-value...

2012-04-25 Thread Uwe Ligges
Probably you got NAs in some bootstrap results, or you got a 0 variance, or you hat 0 observations in a group you compared to another group? Check the data and the separate results! Uwe Ligges On 25.04.2012 09:28, O wrote: Hello, I am comparing two ROC curves with bootstraping. However, so

[R] [pROC] roc.test returns "NA" p-value...

2012-04-25 Thread O
Hello, I am comparing two ROC curves with bootstraping. However, some runs return "p-value = NA," and I have no clue why. Is this anyhow related to like sample size or no sufficient numbers of bootstraping? I used the default value (i.e. boot.n=2000), and the number of observations are quite bi

Re: [R] Proc Mixed to R

2011-10-28 Thread Abhijit Dasgupta
You need to use either the lme4 or nlme packages for mixed models. (There are some other possibilities as well). See http://glmm.wikidot.com/faq for MUCH more detail On 10/27/2011 7:19 PM, Molly Hanlon wrote: > Hi All, > > I'm working with some SAS code to analyze an experiment set up as follows

[R] Proc Mixed to R

2011-10-27 Thread Molly Hanlon
Hi All, I'm working with some SAS code to analyze an experiment set up as follows: 66 subjects (colonies) treated with a random treatment (1-8) and measured at three time points. The data structure looks like: input colony tmt y1 y2 y3; y=y1; date=*1*; output;

Re: [R] 'PROC CONTENTS' = ?

2011-06-10 Thread Frank Harrell
See the contents function in the Hmisc package. karena wrote: > > Is there any function in R which does the same thing as what 'PROC > CONTENTS' does in SAS? > > thanks, > > karena > - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: htt

Re: [R] PROC NLMIXED what package equivalent in R?

2011-03-09 Thread Bert Gunter
Google on "R mixed effects" ! (or you could use R's search tools -- ?help ) -- Bert On Wed, Mar 9, 2011 at 4:58 PM, Chris Buddenhagen wrote: > To account for likely differences between > families in naturalization rates, we fitted a > generalized linear mixed model, using > PROC NLMIXED in SAS

[R] PROC NLMIXED what package equivalent in R?

2011-03-09 Thread Chris Buddenhagen
To account for likely differences between families in naturalization rates, we fitted a generalized linear mixed model, using PROC NLMIXED in SAS10, with the naturalization rate per genus (that is, the number of naturalized species in a genus as a proportion of the total number of introduced specie

Re: [R] PROC COMPARE

2010-11-23 Thread Joshua Wiley
Hi James, Here is another option: all.equal() is more general, the sapply() option I wrote assumes both data frames are basically the same, but that maybe a few elements are different (e.g., you have two data sets and are not sure which is more recent). y <- x <- iris # x & y copies of iris ## c

Re: [R] PROC COMPARE

2010-11-23 Thread Henrique Dallazuanna
See all.equal, a toy example: all.equal(iris, Orange) On Tue, Nov 23, 2010 at 3:17 PM, Jim Edd Jones wrote: > Is it possible to use R to compare two datasets to look for discrepancies, > as one would with the SAS procedure PROC COMPARE? > > Any help on this would be greatly appreciated. > > Ja

[R] PROC COMPARE

2010-11-23 Thread Jim Edd Jones
Is it possible to use R to compare two datasets to look for discrepancies, as one would with the SAS procedure PROC COMPARE? Any help on this would be greatly appreciated. James E. Jones 135 Salina St. Lafayette, CO 80026 [[alternative HTML version deleted]] __