Lui estimators

Raaed

How can Liu's estimations be calculated using JMP pro?!

 

Attached
Real data obtained from the Central Child Hospital in Baghdad - Iraq, representing the number of cases of congenital defects of children in the heart and circulatory system for the period (2006-2011).

 

regards

Poi Reg .xlsx

Mark_Bailey

You send the individual vectors but then you refer to the data table in the liureg function call. You didn't send the data table as a data frame. You didn't make a data frame in R from the vectors. There are several approaches but it seems to me that the easiest way is to send the JMP data table as a data frame instead of individual vectors. I should you how to do it both ways but you mixed them together and left out the important parts in the process.

 

Try this script:

 

Names Default to Here( 1 );

// make sure JMP can find the desired installed version of R
Set Environment Variable( "R_HOME", "C:\Program Files\R\R-4.1.3" );

// Open Data Table: PoiReg.jmp
dta = Open( "/D:/Rdata/PoiReg.jmp" );

// connect to R session
R Init( );

// send JMP data as vectors
R Send( dta );

// use centered function in liureg package with entire data set
R Submit("
library(Matrix)
library(glmnet)
library(liureg)
library (readxl)
## liu regression

lu.reg <- liu(y~., data=dta, d=seq(0, 0.1, 0.1), scaling= \!"centered\!")
print(lu.reg)
");

// close R connection
R Term();