obsnoMrss: observation numbers based on classical and modified ranked...

In RSSampling: Ranked Set Sampling

Description

The obsno.Mrss function gives the observation numbers to sample from a target population by using modified ranked set sampling methods. Ranking is done using the concomitant variable Y.

Usage

  obsno.Mrss(Y,m,r=1,type="r",p)

Arguments

Y	A vector of concomitant variable from target population
m	Size of units in each set
r	Number of cycles
type	type of the modified RSS method. "r" for traditional RSS, "p" for Percentile RSS, "m" for Median RSS, "bg" for Balanced Groups RSS, "e" for Extreme RSS. Default value is "r"
p	Value of percentile for Percentile RSS method

Details

Concomitant variable Y must be a vector.

References

McIntyre, G. A. (1952). A method for unbiased selective sampling, using ranked sets. Australian Journal of Agricultural Research, 3(4), 385-390.

Dell, T. R., & Clutter, J. L. (1972). Ranked set sampling theory with order statistics background. Biometrics, 28, 545-553.

Samawi, H. M., Ahmed, M. S., & Abu-Dayyeh, W. (1996). Estimating the population mean using extreme ranked set sampling. Biometrical Journal, 38(5), 577-586.

Muttlak, H. A. (1997). Median ranked set sampling. Journal of Applied Statistical Sciences, 6(4), 245-255.

Muttlak, H. A. (2003). Modified ranked set sampling methods. Pakistan Journal Of Statistics, 19(3), 315-324.

Jemain, A. A., Al-Omari, A., & Ibrahim, K. (2008). Some variations of ranked set sampling. Electronic Journal of Applied Statistical Analysis, 1(1), 1-15.

See Also

con.Mrss, Mrss, rss

Examples

  y=rexp(10000)
  ## Determining the observation numbers of the units which are chosen to sample

  y=rexp(10000)
  obsno.Mrss(y,m=3,r=5)
  obsno.Mrss(y,m=5,r=6,type="m")
  obsno.Mrss(y,m=7,r=3,type="e")
  obsno.Mrss(y,m=4,r=5,type="p",p=0.3)
  obsno.Mrss(y,m=6,r=2,type="bg")

Example output

      m = 1        m = 2        m = 3       
r = 1 "Obs.  9359" "Obs.  505"  "Obs.  7580"
r = 2 "Obs.  9512" "Obs.  8612" "Obs.  9678"
r = 3 "Obs.  3654" "Obs.  6857" "Obs.  3801"
r = 4 "Obs.  9230" "Obs.  7039" "Obs.  2058"
r = 5 "Obs.  796"  "Obs.  349"  "Obs.  3270"
      m = 1        m = 2        m = 3        m = 4        m = 5       
r = 1 "Obs.  385"  "Obs.  2529" "Obs.  2880" "Obs.  6767" "Obs.  9038"
r = 2 "Obs.  8337" "Obs.  6987" "Obs.  6589" "Obs.  1188" "Obs.  2543"
r = 3 "Obs.  1181" "Obs.  231"  "Obs.  7355" "Obs.  1493" "Obs.  3041"
r = 4 "Obs.  8288" "Obs.  8009" "Obs.  6616" "Obs.  455"  "Obs.  4846"
r = 5 "Obs.  1141" "Obs.  245"  "Obs.  8679" "Obs.  5896" "Obs.  3974"
r = 6 "Obs.  7636" "Obs.  8853" "Obs.  6039" "Obs.  4655" "Obs.  8275"
      m = 1        m = 2        m = 3        m = 4        m = 5       
r = 1 "Obs.  961"  "Obs.  1944" "Obs.  6795" "Obs.  513"  "Obs.  3523"
r = 2 "Obs.  4997" "Obs.  4285" "Obs.  6495" "Obs.  9008" "Obs.  3313"
r = 3 "Obs.  6671" "Obs.  5928" "Obs.  8344" "Obs.  3635" "Obs.  3221"
      m = 6        m = 7       
r = 1 "Obs.  2243" "Obs.  9315"
r = 2 "Obs.  9668" "Obs.  6164"
r = 3 "Obs.  6219" "Obs.  5650"
      m = 1        m = 2        m = 3        m = 4       
r = 1 "Obs.  1827" "Obs.  6551" "Obs.  8065" "Obs.  4760"
r = 2 "Obs.  6527" "Obs.  7624" "Obs.  1056" "Obs.  8952"
r = 3 "Obs.  31"   "Obs.  4650" "Obs.  4008" "Obs.  8132"
r = 4 "Obs.  2231" "Obs.  2087" "Obs.  1496" "Obs.  4796"
r = 5 "Obs.  3744" "Obs.  1032" "Obs.  9315" "Obs.  1122"
      m = 1        m = 2        m = 3        m = 4          m = 5       
r = 1 "Obs.  7183" "Obs.  3188" "Obs.  5852" "Obs.  8126"   "Obs.  7162"
r = 2 "Obs.  8793" "Obs.  685"  "Obs.  6523" "Obs.  3833.5" "Obs.  268" 
      m = 6       
r = 1 "Obs.  2995"
r = 2 "Obs.  9492"

RSSampling documentation built on May 2, 2019, 4:28 a.m.