David:
Its' not going my way today.
I'll try to work with all that you provided.
It's late Friday and I cannot impose on you much more.
Thank you.
If I may, after I rework it, I'll contact you, okay?
Have great weekend. And thanks, again.
Bruce
FYI: Below are where the error messages start:
R> Decile_RespRate <- (No_Resp/No_Inds) Error: object 'No_Resp' not
found R> R> dec_mean_wt_R_nRD <- cbind(dec_mean_wt_R_nR, Cum_RespRate,
Decile_RespRate) Error in cbind(dec_mean_wt_R_nR, Cum_RespRate,
Decile_RespRate) : object 'Decile_RespRate' not found R> R> R> R> avg_RR
<- dec_mean_wt_R_nRD[10,7] Error: object 'dec_mean_wt_R_nRD' not found
R> Cum_Lift <- (Cum_RespRate/avg_RR)*100 Error: object 'avg_RR' not
found R> R> DECILE <- c("top","2","3","4","5","6","7","8","9","bot") R>
R> dec_mean_wt_R_nRDL <- cbind(DECILE, dec_mean_wt_R_nRD, Cum_Lift,
stringsAsFactors=FALSE) Error in cbind(DECILE, dec_mean_wt_R_nRD,
Cum_Lift, stringsAsFactors = FALSE) : object 'dec_mean_wt_R_nRD' not
found R> dec_mean_wt_R_nRDL <- dec_mean_wt_R_nRDL[,c(1,3,4,9,8,10)]
Error: object 'dec_mean_wt_R_nRDL' not found R> R> R>
total_line<-cbind(DECILE="Total", +
as.data.frame(matrix(c(colSums(dec_mean_wt_R_nRDL[ , 2:3]), rep(NA,
3)),nrow=1))) Error in is.data.frame(x) : object 'dec_mean_wt_R_nRDL'
not found R> R> names(total_line)<-names(dec_mean_wt_R_nRDL) Error:
object 'dec_mean_wt_R_nRDL' not found R>
dec_mean_wt_R_nRDLT<-rbind(dec_mean_wt_R_nRDL,total_line) Error in
rbind(dec_mean_wt_R_nRDL, total_line) : object 'dec_mean_wt_R_nRDL' not
found R> decile_table <- dec_mean_wt_R_nRDLT Error: object
'dec_mean_wt_R_nRDLT' not found R> decile_table Error: object
'decile_table' not found R> R> #Install the xtable package:
install.packages("xtable") R> #Load the xtable package: R>
library(xtable) R> R> DECILE_TABLE <-xtable(decile_table) Error in
xtable(decile_table) : object 'decile_table' not found R> DECILE_TABLE
Error: object 'DECILE_TABLE' not found R> R> R>
print.xtable(DECILE_TABLE,
type="html",file="C:/R_Data/DecileTable.html", include.rownames=FALSE)
Error in print.xtable(DECILE_TABLE, type = "html", file =
"C:/R_Data/DecileTable.html", : object 'DECILE_TABLE' not found
R>
Bruce Ratner, Ph.D.
The Significant Statistician™
(516) 791-3544
Statistical Predictive Analtyics -- www.DMSTAT1.com
Machine-Learning Data Mining and Modeling -- www.GenIQ.net
David L Carlson wrote:
I've attached a modification of your script file (called .txt so it doesn't get
stripped). See if this does what you want.
David C
-----Original Message-----
From: Bruce Ratner PhD [mailto:b...@dmstat1.com]
Sent: Friday, April 21, 2017 3:46 PM
To: David L Carlson <dcarl...@tamu.edu>
Cc: r-help@r-project.org
Subject: Re: [R] Looking for a package to replace xtable
David:
Correction: I do need a data frame because from the order Response column I
create a data frame as per all my calculated columns, yielding the five column
decile table.
I used your latest corrections but something buggy is a happenin'.
Bruce
______________
Bruce Ratner PhD
The Significant Statistician™
(516) 791-3544
Statistical Predictive Analytics -- www.DMSTAT1.com
Machine-Learning Data Mining -- www.GenIQ.net
On Apr 21, 2017, at 4:25 PM, Bruce Ratner PhD <b...@dmstat1.com> wrote:
David:
Response=rbinom(50,1,0.2), and yhat=runif(50) are simulating the output of a
say logistic model, where Response is actual 0-1 responses, and yhat is the
predicted
response variable.
I usually resample the original data to get some noise out of the data. I find
it valuable if I can resample from a large sample than the original.
(I know this is viewed by some as unorthodox.)
Your point: I only need Response as a column vector.
That said, what would you alter, please?
Thanks for your time.
Regards,
Bruce
______________
Bruce Ratner PhD
The Significant Statistician™
(516) 791-3544
Statistical Predictive Analytics -- www.DMSTAT1.com
Machine-Learning Data Mining -- www.GenIQ.net
On Apr 21, 2017, at 3:43 PM, David L Carlson <dcarl...@tamu.edu> wrote:
You have an issue at the top with
Resp <- data.frame(Response=rbinom(50,1,0.2), yhat=runif(50))
Resp <- Resp[order(Response$yhat,decreasing=TRUE),]
Since Response$yhat has not been defined at this point. Presumably you want
Resp <- Resp[order(Resp$yhat,decreasing=TRUE),]
The main issue is that you have a variable Response that is located in a data
frame called ResponseX10.
In creating cum_R you need
cum_R <- with(ResponseX10, cumsum(Response))
then dec_mean
dec_mean <- with(ResponseX10, aggregate(Response, by=list(decc), mean))
then dd
dd <- with(ResponseX10, cbind(Response, dd_))
You might consider if Response really needs to be inside a data frame that
consists of a single column (maybe you do if you need to keep track of the row
numbers). If you just worked with the vector Response, you would not have to
use with() or attach().
I'm not sure what the first few lines of your code are intended to do. You
choose random binomial values and uniform random values and then order the
first by the second. But rbinom() is selecting random values so what is the
purpose of randomizing random values? If the real data consist of a vector of
1's and 0's and those need to be randomized, sample(data) will do it for you.
Then those numbers are replicated 10 times. Why not just select 500 values
using rbinom() initially?
David C
-----Original Message-----
From: BR_email [mailto:b...@dmstat1.com]
Sent: Friday, April 21, 2017 1:22 PM
To: David L Carlson <dcarl...@tamu.edu>; r-help@r-project.org
Subject: Re: [R] Looking for a package to replace xtable
David:
I tried somethings and got a little more working.
Now, I am struck at last line provided: "dec_mean <-
aggregate(Response ~ decc, dd, mean)"
Any help is appreciated.
Bruce
*****
Resp <- data.frame(Response=rbinom(50,1,0.2), yhat=runif(50))
Resp <- Resp[order(Response$yhat,decreasing=TRUE),]
ResponseX10 <- do.call(rbind, replicate(10, Resp, simplify=FALSE))
str(ResponseX10)
ResponseX10 <- ResponseX10[order(ResponseX10$yhat,decreasing=TRUE),]
str(ResponseX10)
head(ResponseX10)
ResponseX10[[2]] <- NULL
ResponseX10 <- data.frame(ResponseX10)
str(ResponseX10)
cum_R <- cumsum(Response)
cum_R
sam_size <- n
______________________________________________
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.
