Rで順序ロジステック回帰をやったときの平行線検定（複数の式の偏回帰係数が等しいかどうかの検定）のやり方がわからなかったが、Rのbrant パッケージのbrant関数で出来るようだ。
以下サンプルコード

library(brant)
library(MASS)

data(survey)
smoke_data <- survey 
smoke_data$Smoke <- ordered(smoke_data$Smoke,levels=c("Never","Occas","Regul","Heavy"))


model1 <- MASS::polr(Smoke ~ Sex + Height, data=smoke_data, Hess=TRUE)
summary(model1)

# 平行線検定
brant(model1)