test.theta: Hypothesis test for proportionality parameter

In sporm: Semiparametric Proportional Odds Rate Model

Description

Hypothesis test for the proportionality parameter of the semiparametric proportional odds rate model.

Usage

test.theta(x, y, alternative = c("two.sided", "less", "greater"), 
    theta = 1, B = 1000, conf.level = 0.95)

Arguments

x,y	Vectors containing the data values of the two samples x1,...,xm and y1,...,yn.
alternative	A character string specifying the alternative hypothesis, must be one of "two.sided" (default), "greater" or "less". You can specify just the initial letter.
theta	The hypothesized values of θ.
B	The number of number Monte Carlo trials for simulation approach.
conf.level	Confidence level of the test.

Details

Using the Monte Carlo simulation method to approximate the p-value of the test statistic which is distribution-free.

Value

theta	Maximum rank-based likelihood estimate of θ
p-value	The p-vakue of the test statistic

Author(s)

Zhong Guan <[e-mail]>

References

Zhong Guan and Cheng Peng (2011), "A rank-based empirical likelihood approach to two-sample proportional odds model and its goodness-of-fit", Journal of Nonparametric Statistics, to appear.

See Also

mrle.sporm.

Examples

# Radar tube life data
z<-RadarTube$Days
v<-RadarTube$Type
x<-z[v==1]; y<-z[v==2]
test.theta(x,y,B=100)

sporm documentation built on Jan. 27, 2021, 5:08 p.m.