Nothing
#'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)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.