[R] Looking for a package to replace xtable

[BR_email]

David:
This is my last query on this issue.
I thank you, again, for your patience,
without perhaps your not understanding what drives (bugs) me.

Is there a more efficient way of writing the code, below,
which is replicating the original dataset Response 5 times?

Response5x <- rbind(Response, Response, Response, Response, Response)

Regards,
Bruce

  

BR_email wrote:
> 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:br at dmstat1.com]
>> Sent: Friday, April 21, 2017 3:46 PM
>> To: David L Carlson <dcarlson at tamu.edu>
>> Cc: r-help at 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 <br at 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 <dcarlson at 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:br at dmstat1.com]
>>>> Sent: Friday, April 21, 2017 1:22 PM
>>>> To: David L Carlson <dcarlson at tamu.edu>; r-help at 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
>