R/single.LogSeq.R

In PeFox/DySeq: Functions for Dyadic Sequence Analyses

#'single.LogSeq
#'
#'prints estimates vor single case out of an LogSeq object, see: \code{\link{LogSeq}}
#'
#'
#'@param x a LogSeq object, that should be printed.
#'@param case determines which case should be shown.
#'
#'@examples
#'data(CouplesCope)
#'my.states<-StateExpand(CouplesCope, 2:49, 50:97)
#'my.trans<-StateTrans(my.states, FALSE)
#'my.logseq<-LogSeq(my.trans, single.case=TRUE)
#'my.logseq
#'single.LogSeq(my.logseq, 41)  # prints estimates for case 41
#'
#'@export

single.LogSeq<-function(x, case){

  beta<-c(x[[1]][case,1], x[[1]][case,2],x[[1]][case,3], x[[1]][case,4])

  exp_beta<-exp(beta)

  p.value<-c(c(x[[4]][case,1], x[[4]][case,2],x[[4]][case,3], x[[4]][case,4]))

  output_seq<-cbind(beta,exp_beta, p.value)
  if(attr(x, "firstSeq")){
    rownames(output_seq)<-c("Intercept", "Actor", "Partner", "Interac")
  }else{
    rownames(output_seq)<-c("Intercept","Partner","Actor","interac")
  }

  output_seq<-as.data.frame(output_seq)
  output<-round(output_seq,3)

  Sys.sleep(0.01)
  print(output)
  cat("\nNote: p.value for intercept is not implemented yet;\np.value of 0 means it is lower than .001!\n\n")

  invisible(output)
}