Probit avec R
De Wiki ODR
Révision datée du 26 avril 2007 à 10:13 par Admin (discussion | contributions)
<?
// probit sous R
// Connexion BD test
//http://localhost/test_probit_cedric.php
$t = time() ;
echo "lancement du traitement : probit sous R<br><br>" ;
flush();
echo " Ouverture en lecture du fichier de commande R<BR><BR>";
flush();
if ( !($fichier_R = fopen ('probit2.par', 'r')))
{
print("impossible d'ouvrir le fichier probit2.par");
exit;
}
fclose ($fichier_R);
$mysqldata = "c:\\easyphp\\mysql\\data\\$base\\";
echo "traitement sous r (long) <BR><BR>" ;
flush();
// exécution de R
$ATraiter = "c:\\rw2001\\bin\\Rterm.exe -q --no-restore --no-save < " . "probit2.par" . ">". "fichier_sortie";
exec ($ATraiter);
// affichage
echo "Impression des resultats dans un tableau <BR><BR>" ;
flush();
if ( !($fic_sortie = fopen ('C:\\Documents and Settings\\cedric\\Bureau\\probitessai.txt', 'r')))
{
print("impossible d'ouvrir le fichier fichier_sortie");
exit;
}
while (!feof($fic_sortie))
{
$ligne = fgets ($fic_sortie, 500);
$ligne = trim($ligne) ;
if($ligne == '' ) continue ;
$ligne = explode(";" ,$ligne);
$var[] = $ligne[0] ;
$estim[] = $ligne[1] ;
$stderror[] = $ligne[2];
$zvalue[] = $ligne[3];
$prz[] = $ligne[4];
}
// entête du tableaux
echo "<TABLE border=0 cellspacing='0' cellpading='5' >";
echo "<tr BGCOLOR = #009900><TD colspan='11'></TD></TR>";
echo "<TR bgcolor = #EAFCEA align='center'><TD BGCOLOR = #009900 width = '0px' ></TD><TD></TD><TD BGCOLOR = #009900 ></TD><TD >Estimate</TD> <TD BGCOLOR = #009900 width = '0px' ></TD> <TD>Std. Error</TD><TD BGCOLOR = #009900 width = '0px' ></TD><TD>z value</TD><TD BGCOLOR = #009900 ></TD><TD>Pr(>|z|)</TD><TD BGCOLOR = #009900 width = '0px' ></TD></TR>";
echo "<tr BGCOLOR = #009900 width = '0px' ><TD colspan='11'></TD></TR>";
// lignes du tableau
$i = 0 ;
while ( $i < count($var) )
{
echo "<TR><TD BGCOLOR = #009900 ></TD><TD>".$var[$i]."</TD>";
echo "<TD BGCOLOR = #009900 ></TD><TD Align = 'right' width = '20%' >".sprintf( "%.7f",$estim[$i])."</TD>";
echo "<TD BGCOLOR = #009900 ></TD><TD Align = 'right' width = '20%' >".sprintf( "%.7f",$stderror[$i])."</TD>";
echo "<TD BGCOLOR = #009900 ></TD><TD Align = 'right' width = '20%' >".sprintf( "%.2f",$zvalue[$i])."</TD>";
echo "<TD BGCOLOR = #009900 ></TD><TD Align = 'right' width = '20%' >".sprintf( "%.3f",$prz[$i])."</TD>";
echo "<TD BGCOLOR = #009900></TD></TR>" ;
$i++;
}
// cloture du tableau
echo "<tr BGCOLOR = #009900><TD colspan='11'></TD></TR>";
echo"</TABLE>";
$t1 = time() ;
$temps = $t1 - $t ;
echo "<BR> $temps secondes" ;
?>

