aide sur un programme en scilab

[invite4db1c5e9]

Bonjour,
j'ai un problème dans mon programme de scilab, pouvez vous m'aider à trouver un programme qui fonctionne qui ferait la comparaison de deux échantillons gaussiens en disant si ils sont gaussiens, en disant si leurs moyennes et leurs variances sont identiques. Merci pour votre aide, j'y suis presque arrivée je pense mais il manque des ligne qui font que tous mon programme est faux si vous desirer je peux vous envoyer mon programme.
Merci encore pour votre aide c'est urgent, c'est pour mon stage.
-----

[GrisBleu]

Salut

tu peux l'afficher ici et on verra.
- pour la moyenne et la variance ce n'est pas complique avec Scilab
- pour la gaussiannite (mot correct ?), regarde d'abord vers les moments d'ordre 4 (ou lit quelques articles la dessus), puis les tests de normalite

@+

[invite4db1c5e9]

Bonjour, je mets mais programme pour la moyenne et la variance merci d'essayer de trouver mes erreurs.Pour la comparaison de deux échantillons gaussiens , c'est bon je pense que j'ai trouvé.

Voici mon programme:
// calcul de la variance empirique sans biais
function SX2=varempsb(X);
n=length(X);
m=moy(X);
SX2= n/(n-1)*(sum(X.^2)/n- m^2);
endfunction

//comparaison de la moyenne de deux échantillons gaussiens
function pval=compmoy(X1,X3,alpha);
n1=length(X1);
n2=length(X3);
X1bar=sum(X1)/length(X1);
X3bar=sum(X3)/length(X3);
SX21=varempsb(X1);
SX23=varempsb(X3);
k=n1+n2-2; //nbre de degré de liberté
t=(X1bar-X3bar)*sqrt(k)/sqrt(1/n1+1/n2);
T=t/sqrt((n1-1)*SX21+(n2-1)*SX23);
vc=cdft("T",k,1-alpha,alpha);
[p,pval]=cdft("PQ",T,k);
endfunction

//comparaison de la variance de deux échantillons gaussiens
function pval2=compvar(X1,X3,alpha);
nu1=length(X1)-1;
nu3=length(X3)-1;
X1bar=sum(X1)/length(X1);
X3bar=sum(X3)/length(X3);
SX21=varempsb(X1);
SX23=varempsb(X3);
F=SX3/SX21;
vc2=cdff("F",nu1,nu3,1-alpha,alpha);
[pp,qq]=cdff("PQ",F,nu1,nu3);
pval2=qq;
endfunction

Merci pour votre aide.

[GrisBleu]

Salut

quelle est la fonction qui pose probleme ? qu est ce qui ne marche pas ? la compilation ou l execution ?

++

[A voir en vidéo sur Futura]

[invitebb20ecbb]

Moi j'arrive pas ici

Un joueur de fortune initiale a (entier > 0) joue contre un autre de fortune initiale b (entier > 0) à pile ou face avec une pièce éventuellement pipée (p différent de 1/2).

La variable aléatoire Sn représente le gain ou la perte du premier à l'instant n. Le jeu s'arrête dès que l'un des deux est ruiné (Sn == -a ou Sn == b).

On désigne par T ce temps d'arrêt aléatoire, écrire une fonction [T, ST] = ruine(a, b, p) des arguments a, b et p, et retournant les valeurs T et ST. S'en servir pour estimer par des moyennes arithmétiques l'espérance de la variable aléatoire T.

function [T,S]=ruine(a,b,p)
S(1) = 0;
T(1)=[];
for i=2:N
if( rand() < p )
X(i)=a-1;
Y(i)=b+1;
else
Y(i)=b-1
X(i)=a+1
end
S(i)=S(i-1)+X(i)
if(S(i)==-a)
T(i)=string("gagne");
//T(i)=string("nonn")
else
if(S(i)==a+b)
T(i)=string("perdu");
// T(i)=string("ok")
end
end
end
endfunction