pairwise (version 0.1.4)

grmSE.poly: Graphical Model Check for polytomous PCM

Description

This function makes the basic calculations for the graphical model check for polytomous itemresponse formats. It is more or less a wraper function, internaly calling the function ipSE.poly. Several splitting options are available (see arguments).

Usage

grmSE.poly(daten, m = max(daten, na.rm = TRUE) + 1,
    teil = "no", splitseed = "no", sortdif = FALSE, ...)

Arguments

daten

daten a data matrix with optionaly named colums (names of items) or a data.frame, potentially with missing values, comprising polytomous responses of respondents (rows) on some items (colums) coded starting with 0 for lowest category to m

m

number of response categories for all items - by default m is defined as m = max(daten,na.rm=TRUE)+1.

teil

Specifies the splitting criterion. Basicly there are three different options available - each with several modes - which are controlled by passing the corresponding character expression to the argument.

1) Using the rawscore for splitting into su

splitseed

numeric, used for set.seed(splitseed) for random splitting - see argument teil

sortdif

logical, In contrast to to the default setting in the function ipSE.poly the argument sortdif here is set to FALSE, so the items will kept in original order - see descripti

...

additional arguments nsample, size, seed, pot for caling ipSE.poly are passed through - see description for

Value

• A (list) object of class grmSEpo containing the item difficulty parameter sigma and their standard errors for two or more subsamples.

A Note on Standard Errors

Estimation of standard errors is done by repeated calculation of item parameters for subsamples of the given data. This procedure is mainly controlled by the arguments nsample and size (see arguments). With regard to calculation time, the argument nsample is the 'time killer'. On the other hand, things (estimation of standard errors) will not necessarily get better when choosing large values for nsample. For example choosing nsample=400 will only result in minimal change for standard error estimation in comparison to (nsample=30) which is the default setting (see examples).

Details

The data is splitted in two or more subsamples and then item thresholds, the parameter (Sigma) and their standard errors (SE) for polytomous items according the PCM are calculatetd for each subsample. Additional arguments (see description of function ipSE.poly) for parameter calculation are passed through.

WARNING: When using data based on booklet designs with systematically missing values (by design) you have to ensure that in each of the booklet the maximum raw value to reach is equal while using the raw value as splitting criterion.

References

description of function ipSE.poly{pairwise}.

Examples

data(bfiN) # loading example data set

# calculating itemparameters and SE for two random allocated subsamples
grmSEpoly<-grmSE.poly(daten=bfiN, teil = "random")

# some examples for plotting options
# plotting item difficulties for two subsamples against each other
# with elipses for a CI = 95\% .
plot(grmSEpoly)

# using triangles as plotting pattern
plot(grmSEpoly,pch=2)

#plotting without CI ellipses
plot(grmSEpoly,ci=0,pch=2)

# plotting with item names
plot(grmSEpoly,itemNames=TRUE)

# Changing the size of the item names
plot(grmSEpoly,itemNames=TRUE, cex.names = 1.3)

# Changing the color of the CI ellipses
plot(grmSEpoly,itemNames=TRUE, cex.names = .8, col.error="green")

###### example from details section 'Some Notes on Standard Errors' ########
# grmSEpoly<-grmSE.poly(daten=bfiN, teil = "random",splitseed=13)
# plot(grmSEpoly)
######
# grmSEpoly_400<-grmSE.poly(daten=bfiN, teil = "random", splitseed=13 ,nsample=400)
# plot(grmSEpoly_400)