spsurvey (version 4.1.4)

revisit_bibd: Create a Balanced Incomplete Block Panel Revisit Design

Description

Create a revisit design for panels in a survey that specifies the time periods for the units of each panel to be sampled based on searching for a D-optimal block design that is a member of the class of generalized Youden designs. The resulting design need not be a balanced incomplete block design. Based on algorithmic idea by Cook and Nachtsheim (1989) and implemented by Robert Wheeler.

Usage

revisit_bibd(n.period, n.pnl, n.visit, nsamp, panel_name = "BIB",
  begin = 1, skip = 1, iter = 30)

Arguments

n.period

Number of time periods for the survey design. Typically, number of periods if sampling occurs once per period or number of months if sampling occurs once per month. (v, number of varieties/treatments in BIBD terms)

n.pnl

Number of panels (b, number of blocks in BIBD terms)

n.visit

Number of time periods to be visited in a panel (k, block size in BIBD terms)

nsamp

Number of samples in each panel.

panel_name

Prefix for name of each panel

begin

Numeric name of first sampling occasion, e.g. a specific period.

skip

Number of sampling occasions to skip between planned sampling periods, e.g., sampling will occur only every 5 periods if skip = 5.

iter

Maximum number of iterations in search for D-optimal Generallized Youden Design.

Value

A two-dimensional array of sample sizes to be sampled for each panel and each sampling occasion.

Details

The function uses find.BIB function from crossdes package to search for a D-optimal block design. crossdes uses package AlgDesign to search balanced incomplete block designs.

References

Cook R. D. and C. Nachtsheim. (1989). Computer-aided blocking of factorial and response-surface designs. Technometrics 31(3), 339-346.

See Also

revisit_dsgn

create a panel revisit design

revisit_rand

create a revisit design with random assignment to panels and time periods

panel_summary

summarize characteristics of a revisit panel design

power.dsgn

power calculation for multiple panel designs

cov.panel.dsgn

covariance matrix for a panel design

plot_powerpaneldesign

plot power curves for panel designs

Examples

# NOT RUN {
# Balanced incomplete block design with 20 sample occasions, 20 panels,
# 3 visits to each unit, and 20 units in each panel.
revisit_bibd(n.period = 20, n.pnl = 20, n.visit = 3, nsamp = 20)

# }