|
|
@@ -159,4 +159,90 @@ desc <- function(html="desc", titre=NULL, table=NULL, variables=NULL, variables_ |
|
|
|
dec() |
|
|
|
HTML("</table>") |
|
|
|
HTMLEnd() |
|
|
|
} |
|
|
|
|
|
|
|
#' Diagnostics de la table courante |
|
|
|
#' |
|
|
|
#' Produit un tableau descriptif de la table courante et des graphiques pour la verification des conditions d'utilisation des tests satistiques |
|
|
|
#' |
|
|
|
#' Permet de produire un tableau descriptif des variables contenues dans la table courante. |
|
|
|
#' Si les labels et formats sont definis et charges ils seront utilises pour peupler le tableau. |
|
|
|
#' La fonction attribue un paramètre "param" à chaque variable à partir de tests de normalité, ce paramètre est modifiable ensuite, il conditionne l'utilisation des tests statistiques dans la fonction comparer(). |
|
|
|
#' Tous les parametres ont une valeur par defaut, il est donc possible de simplement lancer diag() pour obtenir le diagnostic par defaut. |
|
|
|
#' @encoding UTF-8 |
|
|
|
#' @param table Table a utiliser, par defaut la derniere chargee |
|
|
|
#' @export |
|
|
|
diag <- function(table=NULL) |
|
|
|
{ |
|
|
|
if (is.null(table)) table<-eval(parse(text=.attachd)) |
|
|
|
|
|
|
|
HTMLInit() |
|
|
|
assign(".param", character(), pos=1) |
|
|
|
|
|
|
|
for (var in names(table)) |
|
|
|
{ |
|
|
|
# Creation du label de la variable |
|
|
|
etiq <- label(table[[var]]) |
|
|
|
if (etiq == "table[[var]]") etiq <- var |
|
|
|
|
|
|
|
HTML("<div class='varblock'>",etiq) |
|
|
|
inc() |
|
|
|
HTML("<div class='varstat'>") |
|
|
|
inc() |
|
|
|
HTML("<table>") |
|
|
|
inc() |
|
|
|
HTML("<tr>", paste("<th>", names(summary(table[[var]])), "</th>", collapse=""), "</tr>") |
|
|
|
HTML("<tr>", sep="") |
|
|
|
for (valeur in summary(table[[var]])) |
|
|
|
{ |
|
|
|
HTML("<td>", valeur, "</td>", sep="") |
|
|
|
} |
|
|
|
HTML("</tr>") |
|
|
|
dec() |
|
|
|
HTML("</table>") |
|
|
|
if (is.numeric(table[[var]])) |
|
|
|
{ |
|
|
|
HTML("IC95 : [ ", format(t.test(table[[var]])$conf.int[1],digits=2,nsmall=2), " - ", format(t.test(table[[var]])$conf.int[2],digits=2,nsmall=2), " ]<br/>") |
|
|
|
HTML("Test de normalité de Shapiro-Wilk : p = ", format(shapiro.test(table[[var]])$p.value,digits=2,nsmall=2), "<br/>") |
|
|
|
if ((shapiro.test(table[[var]])$p.value > 5) & (length(na.omit(table[[var]])) >= 30)) |
|
|
|
{ |
|
|
|
assign(".param", c(.param,var)) |
|
|
|
HTML("<b>Paramétrique</b>") |
|
|
|
} |
|
|
|
} |
|
|
|
dec() |
|
|
|
HTML("</div>") |
|
|
|
|
|
|
|
HTML("<div class='varplot'>") |
|
|
|
inc() |
|
|
|
if (is.numeric(table[[var]])) |
|
|
|
{ |
|
|
|
pngfile=tempfile(pattern="figure",fileext=".png") |
|
|
|
png(filename=pngfile) |
|
|
|
hist(table[[var]],freq=F,col="lightblue",border="darkblue",xlab=etiq,ylab="Densité",main=paste0("Distribution de ",etiq)) |
|
|
|
lines(density(table[[var]],na.rm=T),col="red") |
|
|
|
dev.off() |
|
|
|
HTML("<img src='", pngfile, "'/>") |
|
|
|
|
|
|
|
pngfile=tempfile(pattern="figure",fileext=".png") |
|
|
|
png(filename=pngfile) |
|
|
|
qqnorm(table[[var]], xlab="Quantiles théoriques",ylab="Quantiles de l'échantillon") |
|
|
|
qqline(table[[var]]) |
|
|
|
dev.off() |
|
|
|
HTML("<img src='", pngfile, "'/>") |
|
|
|
} |
|
|
|
else if (is.factor(table[[var]])) |
|
|
|
{ |
|
|
|
pngfile=tempfile(pattern="figure",fileext=".png") |
|
|
|
png(filename=pngfile) |
|
|
|
barplot(table(table[[var]]), col="lightblue", border="darkblue") |
|
|
|
dev.off() |
|
|
|
HTML("<img src='", pngfile, "'/>") |
|
|
|
} |
|
|
|
dec() |
|
|
|
HTML("</div>") |
|
|
|
dec() |
|
|
|
HTML("</div>") |
|
|
|
} |
|
|
|
HTMLEnd() |
|
|
|
} |