cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 

Discussions

Solve problems, and share tips and tricks with other JMP users.
Choose Language Hide Translation Bar
acton
Level I

how to calculate relative potency in parallel line model?

Hi,

I fitted the parallel line model LogY=Sample+Log_Dose in fit curve. 

The parallelism test and model parameters are fine. They are exactly same as what I got in R or PLA software.  But the relative potencies are totally different. Any place to find the technical details about its calculation? Thank you! 

 

acton_1-1763246954709.png

 

 

3 REPLIES 3
hogi
Level XIII

Re: how to calculate relative potency in parallel line model?

A small test data table, the different codes used and an overview of the results will facilitate the investigation.

acton
Level I

Re: how to calculate relative potency in parallel line model?

I used the example from PLA. The RP from R is 2.170981, which is the same as what I got from PLA software. But JMP's result is 5.9..

The estimated parameters are the same, so I think both are using the same model, just the RP from JMP is too large.

 

data <- read_excel("PLA Example_parallel line.xlsx", sheet = "Sheet1")
data$Log_Dose = log(data$Dose)
data$LogY = log(data$Y)

fit <- lm(LogY ~ Sample+Log_Dose , data = data) #parallel line model: common slope
summary(fit)
coefficients <- coef(fit)
db <- coefficients["SampleT"] #
beta <- coefficients["Log_Dose"]
LogRP <- db / beta
RP<-exp(logRP)

 

acton_5-1763325132129.png

 

In JMP, I fitted the model in specialized modeling-> fit curve, polynomials->Fit linear, and test parallelism 

acton_0-1763324922076.png

acton_1-1763325005662.png

acton_2-1763325047117.pngacton_3-1763325089853.png

 

acton_4-1763325100944.png

 

 

 

jthi
Super User

Re: how to calculate relative potency in parallel line model?

Recommended Articles