# vci: confidence limits for method of moments estimators of... In daewr: Design and Analysis of Experiments with R

## Description

function for getting confidence intervals on variance components estimated by the method of moments

## Usage

 `1` ```vci(confl,c1,ms1,nu1,c2,ms2,nu2) ```

## Arguments

 `confl` input- confidence level `c1` input - linear combination coefficient of ms1 in the estimated variance component `ms1` input - Anova mean square 1 `nu1` input - Anova degrees of freedom for mean square 1 `c2` input - linear combination coefficient of ms2 in the estimated variance component `ms2` input - Anova mean square 2 `nu2` input - Anova degrees of freedom for mean square 2

## Value

returned delta, Lower and Upper limits

John Lawson

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24``` ```vci(.90,.05,.014852,2,.05,.026885,18) ## The function is currently defined as vci<-function(confl,c1,ms1,nu1,c2,ms2,nu2){ delta<-c1*ms1-c2*ms2 alpha<-1-confl Falpha1<-qf(confl,nu1,10000000) Falpha12<-qf(confl,nu1,nu2) Fconf2<-qf(alpha,nu2,10000000) Fconf12<-qf(alpha,nu1,nu2) Falpha2<-qf(confl,nu2,10000000) Fconf1<-qf(alpha,nu1,10000000) Fconf12<-qf(alpha,nu1,nu2) G1<-1-(1/Falpha1) H2<-(1/Fconf2)-1 G12<-((Falpha12-1)**2-G1**2*Falpha12**2-H2**2)/Falpha12 VL<-G1**2*c1**2*ms1**2+H2**2*c2**2*ms2**2+G12*c1*c2*ms1*ms2 H1<-(1/Fconf1)-1 G2<-1-(1/Falpha2) H12<-((1-Fconf12)**2-H1**2*Fconf12**2-G2**2)/Fconf12 VU<-H1**2*c1**2*ms1**2+G2**2*c2**2*ms2**2 L<-delta-sqrt(VL) U<-delta+sqrt(VU) cat("delta=",delta," Lower Limit=",L," Upper Limit=",U,"\n") } ```

daewr documentation built on May 29, 2017, 6:11 p.m.