Description Usage Arguments Details Value Author(s) References See Also Examples
Hypothesis test for the proportionality parameter of the semiparametric proportional odds rate model.
1 2 | test.theta(x, y, alternative = c("two.sided", "less", "greater"),
theta = 1, B = 1000, conf.level = 0.95)
|
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. |
Using the Monte Carlo simulation method to approximate the p-value of the test statistic which is distribution-free.
theta |
Maximum rank-based likelihood estimate of θ |
p-value |
The p-vakue of the test statistic |
Zhong Guan <zguan@iusb.edu>
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.
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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.