Running ‘Becker-Clogg1989.R’ [0s/1s]
  Running ‘Clogg-Shihadeh1994.R’ [2s/3s]
  Running ‘Goodman1979.R’ [2s/3s]
  Running ‘Goodman1991.R’ [2s/3s]
  Running ‘Wong2010-2.7.R’ [0s/1s]
  Running ‘Wong2010-4.7.R’ [0s/1s]
  Running ‘Yamaguchi1990.R’ [2s/3s]
  Running ‘hmskewL.R’ [4s/5s]
  Running ‘iac.R’ [7s/8s]
  Running ‘maor.R’ [5s/6s]
  Running ‘ras.R’ [2s/3s]
  Running ‘supplementary.R’ [2s/3s]
  Running ‘vanderHeijden-Mooijaart1995.R’ [2s/3s]
Running the tests in ‘tests/supplementary.R’ failed.
Complete output:
  > ## Variations based on Becker & Clogg (1989), cf. Becker-Clogg1989.R
  > 
  > library(logmult)
  Loading required package: gnm
  
  Attaching package: 'logmult'
  
  The following object is masked from 'package:gnm':
  
      se
  
  > data(color)
  > 
  > caithness <- color[,,"Caithness"]
  > 
  > 
  > ## Adding the table again as supplementary rows
  > rsup <- caithness
  > csup <- caithness
  > rownames(rsup) <- paste(rownames(caithness), "2")
  > colnames(csup) <- paste(colnames(caithness), "2")
  > 
  > names(dimnames(rsup)) <- names(dimnames(csup)) <- names(dimnames(caithness))
  > 
  > mA1r <- rc(caithness, 1, rowsup=rsup, weighting="u", start=NA)
  Running base model to find starting values...
  Running real model...
  Initialising
  Running main iterations.....................
  Done
  > mA1c <- rc(caithness, 1, colsup=csup, weighting="u", start=NA)
  Running base model to find starting values...
  Running real model...
  Initialising
  Running main iterations.....................
  Done
  > mA1 <- rc(caithness, 1, rowsup=rsup, colsup=csup, weighting="u", start=NA)
  Running base model to find starting values...
  Running real model...
  Initialising
  Running main iterations.....................
  Done
  > 
  > rscoresA1 <- mA1$assoc$row[1:4, 1, 1] * sqrt(mA1$assoc$phi[1, 1])
  > cscoresA1 <- mA1$assoc$col[1:5, 1, 1] * sqrt(mA1$assoc$phi[1, 1])
  > 
  > rmA1 <- gnm(Freq ~ Eye + Hair + Eye:cscoresA1[col(rsup)],
  +                data=as.table(rsup), family=poisson)
  > cmA1 <- gnm(Freq ~ Eye + Hair + Hair:rscoresA1[row(csup)],
  +                data=as.table(csup), family=poisson)
  > 
  > rscA1 <- getContrasts(rmA1, pickCoef(rmA1, "cscores"), ref="mean")$qvframe[,1]
  > cscA1 <- getContrasts(cmA1, pickCoef(cmA1, "rscores"), ref="mean")$qvframe[,1]
  > 
  > stopifnot(all.equal(rscA1, rscoresA1, mA1$assoc$row[5:8,,] * sqrt(mA1$assoc$phi[1, 1]),
  +                     tolerance=1e-6, check.attributes=FALSE))
  Error in all.equal(rscA1, rscoresA1, mA1$assoc$row[5:8, , ] * sqrt(mA1$assoc$phi[1,  : 
    all(scale > 0) is not TRUE
  Calls: stopifnot ... tryCatch -> tryCatchList -> tryCatchOne -> <Anonymous>
  Execution halted