test.theta: Hypothesis test for proportionality parameter

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/test.theta.R

Description

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

Usage

1
2
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 <zguan@iusb.edu>

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

1
2
3
4
5
# 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 May 2, 2019, 8:56 a.m.