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)

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




