getPred {hbmem}R Documentation

Function getPred

Description

getPred is a helper function that takes a block and its indexes, and returns the predicted value from that block for each trial. It is the same thing as taking $(XtX)^(-1) Xty$ where X is the design matrix for 1 (mu),sub, item, and lag; y is the vector of effect estimates, block.

Usage

getPred(block,cond,sub,item,lag,N,I,J,R)

Arguments

block Block of parameter values starting with the grand mean, followed by participant effects, and finially item effects.
cond Vector of condition index, starting at zero.
sub Vector of subject index, starting at zero.
item Vector of item index, starting at zero.
lag Vector of lag index, zero-centered.
N Number of conditions.
I Number of subjects.
J Number of items.
R Total number of trials.

Value

Returns predicted means for each trial.

Author(s)

Michael S. Pratte

See Also

hbmem

Examples

library(hbmem)
sim=uvsdSim()
N=length(levels(as.factor(sim@cond[sim@Scond==0])))
I=length(levels(as.factor(sim@subj)))
J=length(levels(as.factor(sim@item)))
RN=sum(sim@Scond==0)

block=c(sim@muN,sim@alphaN,sim@betaN,0,0,-.001)

means=getPred(block,sim@cond[sim@Scond==0],sim@subj[sim@Scond==0],
sim@item[sim@Scond==0],sim@lag[sim@Scond==0],N,I,J,RN)

multi=sim@resp[sim@Scond==0] 

plot(multi~means,ylab="Response to New Item",xlab="Mean of New-Item Distribution")


[Package hbmem version 0.2 Index]