R/CTT.R
In leandromarino/CesgTools: Educational Assessment Toolkit (Cesgranrio Foundation)
#' Calcula as Estatisticas Classicas
#'
#' @description
#' Esta funcao calcula as Estatisticas Classicas, score, score por bloco,
#' percentual de acertos para um conjunto de dados e indica itens que possam
#' ter algum tipo de problema como, por exemplo: gabarito incorreto, item
#' mal elaborado, problema de impressao, entre outros.
#'
#' @param data conjunto de dados a ser utilizado nas analises
#' @param ids vetor com os nomes das colunas de \code{data} que serao utilizadas
#' como identificacao na saida dos scores.
#' @param vars vetor com os nomes das colunas de \code{data} na ordem
#' \code{c(caderno,respostas,peso (se houver))}.
#' @param nitems Numero de itens total.
#' @param nitemform Numero de itens em cada caderno.
#' @param nforms Numero de cadernos diferentes.
#' @param nblform Numero de blocos em que a prova estua organizada. Se eh uma prova
#' fixa, \code{nblform=1}.
#' @param nbl numero de blocos totais.7
#' @param tbl Tamanho do bloco, ou seja, numero de itens existentes no bloco.
#' Se eh uma prova fixa, \code{tbl=nitemform}.
#' @param gab Vetor com a string de gabaritos, cada item eh um elemento deste
#' vetor.
#' @param resp_possible Respostas possiveis (por exemplo:
#' \code{resp_possible = c(LETTERS[1:5],' ','*')}).
#' @param items Objeto do tipo matriz com \code{rownames} e
#' \code{colnames} definido.
#' @param ndec Numero de decimais para a saida das estatisticas classicas.
#' Default: \code{ndec=2}.
#' @param acer Codigo para a representacao do acerto no arquivo de saida.
#' Default: \code{acer='1'}
#' @param erro Codigo para a representacao do erro no arquivo de saida.
#' Default: \code{erro='0'}
#' @param napres Codigo para a representacao do item, quando nao apresentado
#' ao estudante, no arquivo de saida. Default: \code{napres='9'}
#' @param peso vetor de pesos: Default: \code{peso = FALSE}
#' @param mostra_napres indicador sobre como os itens nao apresentados aos alunos
#' devem ser exibidos. Default: \code{mostra_napres = 0}, ou seja o não apresentado
#' é desconsiderado das médias dos alunos. \code{mostra_napres = 1} o não apresentado
#' entra nas contas como se fosse uma resposta possível.#'
#'
#' @return O arquivo de saida desta funcao eh composto por 4 objetos
#' dispostos em uma lista.
#' @return CTT[[1]] \code{data.frame} com as estatisticas classicas para
#' todos os itens. Cada item eh uma linha deste objeto.
#' @return CTT[[2]] \code{data.frame} contendo: padrao de respostas:
#' acerto = \code{acer}, erro = \code{erro} e
#' item nao apresentado = \code{napres}; numero de acertos, numero de itens
#' realizados pelo aluno, percentual de acertos e numero de acertos por bloco.
#' @return CTT[[3]] \code{data.frame} com as estatiticas classicas para os
#' itens que podem ter apresentado algum problema. No caso, itens com 2 ou
#' mais coeficientes bisseriais positivos e itens com coeficiente biserial
#' abaixo de 0.15
#'
#' @examples
#' data("quest_dados")
#' data("gabpar05P")
#'
#'
#' item_prova <- matrix(c(1:36,13:36,1:12,25:36,1:24),byrow = T,ncol=36)
#' colnames(item_prova) <- paste('it',sprintf("%02d",1:36),sep='')
#' rownames(item_prova) <- 1:nrow(item_prova)
#'
#' dadosP <- quest_dados[nchar(quest_dados$rsp_por)==36,]
#' dadosP <- quest_dados[nchar(quest_dados$rsp_por)==36 & !is.na(quest_dados$cad_por),]
#'
#' CTT(
#' data = dadosP,
#' ids = c('codesc','turma','id'),
#' vars = c('cad_por','rsp_por'),
#' peso = FALSE,
#' nitems = 36,
#' nitemform = 36,
#' nforms = 3,
#' gab = gabpar05P$gab,
#' resp_possible = c(LETTERS[1:4]," ","*"),
#' items = item_prova,
#' ndec = 2,
#' nblform = 3,
#' tbl = 12,
#' nbl = 3,
#' acer = '1',
#' erro = '0',
#' napres = '9',
#' calc_normit = F)
#'
#' @seealso \code{\link{RespItem}},\code{\link{Escore}},
#' \code{\link{PontoBisserial}}, \code{\link{ItemPos}},
#' \code{\link{write.ctt}}
#' @export
# data = conjunto de dados a ser utilizados
# ids = nome das colunas que identificam os dados
# vars = variaveis que serao utilizadas na ordem: caderno e respostas
CTT <- function (data, ids, vars, peso = FALSE, nitems, nforms, nitemform,
nblform, nbl, tbl, gab, resp_possible, items, ndec = 2,
acer = "1", erro = "0", napres = "9", mostra_napres = 0,
calc_normit = FALSE)
{
if (sum(is.element(colnames(data), ids)) != length(ids))
stop("Alguma das colunas de ids não existem nos dados")
if (sum(is.element(colnames(data), vars)) != length(vars))
stop("Alguma das colunas de vars não existem nos dados")
if (sum(nchar(data[, vars[2]]) == 0) > 0)
stop("sum(nchar(responses)==0)>0 \n Um ou mais elementos de responses tem tamanho 0")
if (length(unique(data[, vars[1]])) > nforms)
stop("length(unique(vecform)) > nforms \n Existem mais tipos de cadernos para os alunos (vecform) do que cadernos informados (form)")
if (peso)
peso = data[, vars[3]]
if (length(peso) == 1) {
if (!peso)
peso = rep(1, nrow(data))
}
nalu <- nrow(data)
if (mostra_napres > 0) {
resp_possible <- c(resp_possible, napres)
}
data <- data[, c(ids, vars)]
colnames(data)[which(names(data) == vars[1])] <- "caderno"
colnames(data)[which(names(data) == vars[2])] <- "respostas"
data[1, ]
nalt <- length(resp_possible)
CTTest <- list()
nCTTest <- 3 * nalt + 5 + 2 + 1
CTTest[[1]] <- data.frame(matrix(NA, ncol = nCTTest, nrow = nitems))
colnames(CTTest[[1]]) <- c("IT", "GAB", "DIFI", "DISCR",
"ABAI", "ACIM", "BISE", "PBISE", paste(rep(c("Perc",
"Bise", "PBise"), c(nalt, nalt, nalt)), rep(resp_possible,
3), sep = ""))
CTTest[[1]]$IT <- 1:nitems
CTTest[[1]]$GAB <- gab
Gabarito <- data.frame(caderno = 1:nrow(items), gabarito = apply(matrix(gab[items],
nrow = nrow(items)), 1, paste0, collapse = ""), stringsAsFactors = F)
CTTest[[2]] <- cbind(data[, c(ids, "caderno", "respostas")],
Escore(respostas = data[, c("caderno", "respostas")],
gabarito = Gabarito, NumCad = nforms, NumItens = nitemform,
CodAcer = acer, CodErro = erro, CodNaoAp = napres,
nblform = nblform, tbl = tbl))
itemposNum <- ItemPos(items, nitems, tipo = "integer")
if (calc_normit) {
(CTTest[[2]]$normit <- normit(scores = CTTest[[2]]$nacer,
caderno = CTTest[[2]]$caderno, nitemform = nitemform,
peso = peso))
}
writeLines("\n Calculando as estatisticas classicas...")
for (i in 1:nitems) {
setTxtProgressBar(txtProgressBar(min = 1, max = nitems,
style = 3, width = 70, initial = 1), i)
auxcols <- paste(rep(c("Perc", "PBise"), c(nalt, nalt)),
rep(resp_possible, 2), sep = "")
CTTest[[1]][i, auxcols] <- unlist(PontoBisserial(respostas = data[,
c("caderno", "respostas")], scores = CTTest[[2]][,
ifelse(calc_normit, "normit", "pacer")], itempos = itemposNum[[i]],
resposta_possivel = resp_possible, CodNaoAp = napres,
Peso = peso, mostra_napres = mostra_napres))
(tp <- CTTest[[1]][i, paste(rep("Perc", nalt), resp_possible,
sep = "")])
(tb <- CTTest[[1]][i, paste(rep("PBise", nalt), resp_possible,
sep = "")])
(CTTest[[1]][i, paste(rep("Bise", nalt), resp_possible,
sep = "")] <- tb * sqrt(tp * (1 - tp))/dnorm(qnorm(as.numeric(tp))))
CTTest[[1]][i, c("DIFI", "BISE", "PBISE")] <- if (gab[i] !=
"X")
CTTest[[1]][i, paste(c("Perc", "Bise", "PBise"),
gab[i], sep = "")]
else rep(NA, 3)
dummyitem <- RespItem(respostas = data[, c("caderno",
"respostas")], itempos = itemposNum[[i]])
if (mostra_napres == 0) {
scoreitem <- CTTest[[2]][, ifelse(calc_normit, "normit",
"pacer")][dummyitem != napres & dummyitem !=
""]
pesomodi <- peso[dummyitem != napres & dummyitem !=
""]
dummyitem <- dummyitem[dummyitem != napres & dummyitem !=
""]
aux <- dummyitem == gab[i]
}
if (mostra_napres > 0) {
scoreitem <- CTTest[[2]][, ifelse(calc_normit, "normit",
"pacer")][dummyitem != ""]
pesomodi <- peso[dummyitem != ""]
dummyitem <- dummyitem[dummyitem != ""]
aux <- dummyitem == gab[i]
}
dummyitem[aux] <- 1
dummyitem[!aux] <- 0
dummyitem <- as.integer(dummyitem)
ginf <- scoreitem <= Hmisc::wtd.quantile(x = scoreitem,
weights = pesomodi, probs = 0.27)
gsup <- scoreitem >= Hmisc::wtd.quantile(x = scoreitem,
weights = pesomodi, probs = 1 - 0.27)
CTTest[[1]][i, "ABAI"] <- sum(dummyitem[ginf] * pesomodi[ginf])/sum(pesomodi[ginf])
CTTest[[1]][i, "ACIM"] <- sum(dummyitem[gsup] * pesomodi[gsup])/sum(pesomodi[gsup])
CTTest[[1]][i, "DISCR"] <- CTTest[[1]][i, "ACIM"] -
CTTest[[1]][i, "ABAI"]
}
cols <- c("DIFI", "DISCR", "ABAI", "ACIM", "BISE", "PBISE",
paste(rep(c("Perc", "Bise", "PBise"), c(nalt, nalt,
nalt)), rep(resp_possible, 3), sep = ""))
CTTest[[1]][, cols] <- round(CTTest[[1]][, cols], ndec)
CTTest[[1]]$BL <- rep(c(1:nbl), rep(tbl, nbl))
CTTest[[1]]$OB <- rep(1:tbl, nbl)
CTTest[[1]]$GAB <- gab
ordemcols <- c("IT", "BL", "OB", "GAB", cols)
CTTest[[1]] <- CTTest[[1]][, ordemcols]
#browser()
aux_bisepos <- CTTest[[1]][rowSums(CTTest[[1]][, paste("Bise",resp_possible[is.letter(resp_possible)], sep = "")] > 0, na.rm = T) > 1, ]
aux_bisepeq <- CTTest[[1]][CTTest[[1]][, "BISE"] <= 0.15 & !is.na(CTTest[[1]][, "BISE"]), ]
aux_probbise <- unique(rbind(aux_bisepos, aux_bisepeq))
if (nrow(aux_probbise) > 0) {
aux_probbise <- aux_probbise[order(aux_probbise$IT),
]
rownames(aux_probbise) <- aux_probbise$IT
aux_probbise$prob <- ""
aux_probbise[as.character(intersect(aux_bisepos$IT,
aux_bisepos$IT)), "prob"] <- "2 Bis > 0 & Bis < .15"
aux_probbise[as.character(setdiff(aux_bisepos$IT, aux_bisepeq$IT)),
"prob"] <- "2 Bis > 0"
aux_probbise[as.character(setdiff(aux_bisepeq$IT, aux_bisepos$IT)),
"prob"] <- "Bis < .15"
}
CTTest[[3]] <- aux_probbise
names(CTTest) <- c("Classical Stats (using score)", "Scores, Block result and vector of corret pattern",
"Classical Stats of Items with problems")
CTTest
}
leandromarino/CesgTools documentation built on May 25, 2022, 5:03 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
#' Calcula as Estatisticas Classicas
#'
#' @description
#' Esta funcao calcula as Estatisticas Classicas, score, score por bloco,
#' percentual de acertos para um conjunto de dados e indica itens que possam
#' ter algum tipo de problema como, por exemplo: gabarito incorreto, item
#' mal elaborado, problema de impressao, entre outros.
#'
#' @param data conjunto de dados a ser utilizado nas analises
#' @param ids vetor com os nomes das colunas de \code{data} que serao utilizadas
#' como identificacao na saida dos scores.
#' @param vars vetor com os nomes das colunas de \code{data} na ordem
#' \code{c(caderno,respostas,peso (se houver))}.
#' @param nitems Numero de itens total.
#' @param nitemform Numero de itens em cada caderno.
#' @param nforms Numero de cadernos diferentes.
#' @param nblform Numero de blocos em que a prova estua organizada. Se eh uma prova
#' fixa, \code{nblform=1}.
#' @param nbl numero de blocos totais.7
#' @param tbl Tamanho do bloco, ou seja, numero de itens existentes no bloco.
#' Se eh uma prova fixa, \code{tbl=nitemform}.
#' @param gab Vetor com a string de gabaritos, cada item eh um elemento deste
#' vetor.
#' @param resp_possible Respostas possiveis (por exemplo:
#' \code{resp_possible = c(LETTERS[1:5],' ','*')}).
#' @param items Objeto do tipo matriz com \code{rownames} e
#' \code{colnames} definido.
#' @param ndec Numero de decimais para a saida das estatisticas classicas.
#' Default: \code{ndec=2}.
#' @param acer Codigo para a representacao do acerto no arquivo de saida.
#' Default: \code{acer='1'}
#' @param erro Codigo para a representacao do erro no arquivo de saida.
#' Default: \code{erro='0'}
#' @param napres Codigo para a representacao do item, quando nao apresentado
#' ao estudante, no arquivo de saida. Default: \code{napres='9'}
#' @param peso vetor de pesos: Default: \code{peso = FALSE}
#' @param mostra_napres indicador sobre como os itens nao apresentados aos alunos
#' devem ser exibidos. Default: \code{mostra_napres = 0}, ou seja o não apresentado
#' é desconsiderado das médias dos alunos. \code{mostra_napres = 1} o não apresentado
#' entra nas contas como se fosse uma resposta possível.#'
#'
#' @return O arquivo de saida desta funcao eh composto por 4 objetos
#' dispostos em uma lista.
#' @return CTT[[1]] \code{data.frame} com as estatisticas classicas para
#' todos os itens. Cada item eh uma linha deste objeto.
#' @return CTT[[2]] \code{data.frame} contendo: padrao de respostas:
#' acerto = \code{acer}, erro = \code{erro} e
#' item nao apresentado = \code{napres}; numero de acertos, numero de itens
#' realizados pelo aluno, percentual de acertos e numero de acertos por bloco.
#' @return CTT[[3]] \code{data.frame} com as estatiticas classicas para os
#' itens que podem ter apresentado algum problema. No caso, itens com 2 ou
#' mais coeficientes bisseriais positivos e itens com coeficiente biserial
#' abaixo de 0.15
#'
#' @examples
#' data("quest_dados")
#' data("gabpar05P")
#'
#'
#' item_prova <- matrix(c(1:36,13:36,1:12,25:36,1:24),byrow = T,ncol=36)
#' colnames(item_prova) <- paste('it',sprintf("%02d",1:36),sep='')
#' rownames(item_prova) <- 1:nrow(item_prova)
#'
#' dadosP <- quest_dados[nchar(quest_dados$rsp_por)==36,]
#' dadosP <- quest_dados[nchar(quest_dados$rsp_por)==36 & !is.na(quest_dados$cad_por),]
#'
#' CTT(
#' data = dadosP,
#' ids = c('codesc','turma','id'),
#' vars = c('cad_por','rsp_por'),
#' peso = FALSE,
#' nitems = 36,
#' nitemform = 36,
#' nforms = 3,
#' gab = gabpar05P$gab,
#' resp_possible = c(LETTERS[1:4]," ","*"),
#' items = item_prova,
#' ndec = 2,
#' nblform = 3,
#' tbl = 12,
#' nbl = 3,
#' acer = '1',
#' erro = '0',
#' napres = '9',
#' calc_normit = F)
#'
#' @seealso \code{\link{RespItem}},\code{\link{Escore}},
#' \code{\link{PontoBisserial}}, \code{\link{ItemPos}},
#' \code{\link{write.ctt}}
#' @export
# data = conjunto de dados a ser utilizados
# ids = nome das colunas que identificam os dados
# vars = variaveis que serao utilizadas na ordem: caderno e respostas
CTT <- function (data, ids, vars, peso = FALSE, nitems, nforms, nitemform,
nblform, nbl, tbl, gab, resp_possible, items, ndec = 2,
acer = "1", erro = "0", napres = "9", mostra_napres = 0,
calc_normit = FALSE)
{
if (sum(is.element(colnames(data), ids)) != length(ids))
stop("Alguma das colunas de ids não existem nos dados")
if (sum(is.element(colnames(data), vars)) != length(vars))
stop("Alguma das colunas de vars não existem nos dados")
if (sum(nchar(data[, vars[2]]) == 0) > 0)
stop("sum(nchar(responses)==0)>0 \n Um ou mais elementos de responses tem tamanho 0")
if (length(unique(data[, vars[1]])) > nforms)
stop("length(unique(vecform)) > nforms \n Existem mais tipos de cadernos para os alunos (vecform) do que cadernos informados (form)")
if (peso)
peso = data[, vars[3]]
if (length(peso) == 1) {
if (!peso)
peso = rep(1, nrow(data))
}
nalu <- nrow(data)
if (mostra_napres > 0) {
resp_possible <- c(resp_possible, napres)
}
data <- data[, c(ids, vars)]
colnames(data)[which(names(data) == vars[1])] <- "caderno"
colnames(data)[which(names(data) == vars[2])] <- "respostas"
data[1, ]
nalt <- length(resp_possible)
CTTest <- list()
nCTTest <- 3 * nalt + 5 + 2 + 1
CTTest[[1]] <- data.frame(matrix(NA, ncol = nCTTest, nrow = nitems))
colnames(CTTest[[1]]) <- c("IT", "GAB", "DIFI", "DISCR",
"ABAI", "ACIM", "BISE", "PBISE", paste(rep(c("Perc",
"Bise", "PBise"), c(nalt, nalt, nalt)), rep(resp_possible,
3), sep = ""))
CTTest[[1]]$IT <- 1:nitems
CTTest[[1]]$GAB <- gab
Gabarito <- data.frame(caderno = 1:nrow(items), gabarito = apply(matrix(gab[items],
nrow = nrow(items)), 1, paste0, collapse = ""), stringsAsFactors = F)
CTTest[[2]] <- cbind(data[, c(ids, "caderno", "respostas")],
Escore(respostas = data[, c("caderno", "respostas")],
gabarito = Gabarito, NumCad = nforms, NumItens = nitemform,
CodAcer = acer, CodErro = erro, CodNaoAp = napres,
nblform = nblform, tbl = tbl))
itemposNum <- ItemPos(items, nitems, tipo = "integer")
if (calc_normit) {
(CTTest[[2]]$normit <- normit(scores = CTTest[[2]]$nacer,
caderno = CTTest[[2]]$caderno, nitemform = nitemform,
peso = peso))
}
writeLines("\n Calculando as estatisticas classicas...")
for (i in 1:nitems) {
setTxtProgressBar(txtProgressBar(min = 1, max = nitems,
style = 3, width = 70, initial = 1), i)
auxcols <- paste(rep(c("Perc", "PBise"), c(nalt, nalt)),
rep(resp_possible, 2), sep = "")
CTTest[[1]][i, auxcols] <- unlist(PontoBisserial(respostas = data[,
c("caderno", "respostas")], scores = CTTest[[2]][,
ifelse(calc_normit, "normit", "pacer")], itempos = itemposNum[[i]],
resposta_possivel = resp_possible, CodNaoAp = napres,
Peso = peso, mostra_napres = mostra_napres))
(tp <- CTTest[[1]][i, paste(rep("Perc", nalt), resp_possible,
sep = "")])
(tb <- CTTest[[1]][i, paste(rep("PBise", nalt), resp_possible,
sep = "")])
(CTTest[[1]][i, paste(rep("Bise", nalt), resp_possible,
sep = "")] <- tb * sqrt(tp * (1 - tp))/dnorm(qnorm(as.numeric(tp))))
CTTest[[1]][i, c("DIFI", "BISE", "PBISE")] <- if (gab[i] !=
"X")
CTTest[[1]][i, paste(c("Perc", "Bise", "PBise"),
gab[i], sep = "")]
else rep(NA, 3)
dummyitem <- RespItem(respostas = data[, c("caderno",
"respostas")], itempos = itemposNum[[i]])
if (mostra_napres == 0) {
scoreitem <- CTTest[[2]][, ifelse(calc_normit, "normit",
"pacer")][dummyitem != napres & dummyitem !=
""]
pesomodi <- peso[dummyitem != napres & dummyitem !=
""]
dummyitem <- dummyitem[dummyitem != napres & dummyitem !=
""]
aux <- dummyitem == gab[i]
}
if (mostra_napres > 0) {
scoreitem <- CTTest[[2]][, ifelse(calc_normit, "normit",
"pacer")][dummyitem != ""]
pesomodi <- peso[dummyitem != ""]
dummyitem <- dummyitem[dummyitem != ""]
aux <- dummyitem == gab[i]
}
dummyitem[aux] <- 1
dummyitem[!aux] <- 0
dummyitem <- as.integer(dummyitem)
ginf <- scoreitem <= Hmisc::wtd.quantile(x = scoreitem,
weights = pesomodi, probs = 0.27)
gsup <- scoreitem >= Hmisc::wtd.quantile(x = scoreitem,
weights = pesomodi, probs = 1 - 0.27)
CTTest[[1]][i, "ABAI"] <- sum(dummyitem[ginf] * pesomodi[ginf])/sum(pesomodi[ginf])
CTTest[[1]][i, "ACIM"] <- sum(dummyitem[gsup] * pesomodi[gsup])/sum(pesomodi[gsup])
CTTest[[1]][i, "DISCR"] <- CTTest[[1]][i, "ACIM"] -
CTTest[[1]][i, "ABAI"]
}
cols <- c("DIFI", "DISCR", "ABAI", "ACIM", "BISE", "PBISE",
paste(rep(c("Perc", "Bise", "PBise"), c(nalt, nalt,
nalt)), rep(resp_possible, 3), sep = ""))
CTTest[[1]][, cols] <- round(CTTest[[1]][, cols], ndec)
CTTest[[1]]$BL <- rep(c(1:nbl), rep(tbl, nbl))
CTTest[[1]]$OB <- rep(1:tbl, nbl)
CTTest[[1]]$GAB <- gab
ordemcols <- c("IT", "BL", "OB", "GAB", cols)
CTTest[[1]] <- CTTest[[1]][, ordemcols]
#browser()
aux_bisepos <- CTTest[[1]][rowSums(CTTest[[1]][, paste("Bise",resp_possible[is.letter(resp_possible)], sep = "")] > 0, na.rm = T) > 1, ]
aux_bisepeq <- CTTest[[1]][CTTest[[1]][, "BISE"] <= 0.15 & !is.na(CTTest[[1]][, "BISE"]), ]
aux_probbise <- unique(rbind(aux_bisepos, aux_bisepeq))
if (nrow(aux_probbise) > 0) {
aux_probbise <- aux_probbise[order(aux_probbise$IT),
]
rownames(aux_probbise) <- aux_probbise$IT
aux_probbise$prob <- ""
aux_probbise[as.character(intersect(aux_bisepos$IT,
aux_bisepos$IT)), "prob"] <- "2 Bis > 0 & Bis < .15"
aux_probbise[as.character(setdiff(aux_bisepos$IT, aux_bisepeq$IT)),
"prob"] <- "2 Bis > 0"
aux_probbise[as.character(setdiff(aux_bisepeq$IT, aux_bisepos$IT)),
"prob"] <- "Bis < .15"
}
CTTest[[3]] <- aux_probbise
names(CTTest) <- c("Classical Stats (using score)", "Scores, Block result and vector of corret pattern",
"Classical Stats of Items with problems")
CTTest
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.