Dev Faq
7
votes

Calculer à l'intérieur et entre les variances et les intervalles de confiance en r

Je dois calculer l'intérieur et entre les variances d'exécution de certaines données dans le cadre de la mise en place d'une nouvelle méthode de chimie analytique. J'ai également besoin d'intervalles de confiance à partir de ces données à l'aide de la langue r

Je suppose que je dois utiliser l'anova ou quelque chose? P>

Mes données sont comme p> 

> variance
   Run Rep Value
1    1   1  9.85
2    1   2  9.95
3    1   3 10.00
4    2   1  9.90
5    2   2  8.80
6    2   3  9.50
7    3   1 11.20
8    3   2 11.10
9    3   3  9.80
10   4   1  9.70
11   4   2 10.10
12   4   3 10.00

0 commentaires

4 Réponses :

4
votes

Si vous souhaitez appliquer une fonction (telle que var code>) sur un facteur tel que code> exécuter code> ou rep code>, vous pouvez utiliser TAPPLY CODE>: 

> with(variance, tapply(Value, Run, var))
          1           2           3           4 
0.005833333 0.310000000 0.610000000 0.043333333 
> with(variance, tapply(Value, Rep, var))
          1          2          3 
0.48562500 0.88729167 0.05583333

1 commentaires

Agréable! C'est un code élégant à mon avis.

10
votes

Vous avez quatre groupes de trois observations: 

> anova(fit)["group", "Mean Sq"]
[1] 0.8586111

1 commentaires

Êtes-vous sûr de votre formule pour la variance entre groupe? Il me semble que ce n'est pas la variance, juste la somme moyenne de la place?

3
votes

Je vais faire une fente à cela quand j'ai plus de temps, mais en attendant, voici le dput () pour la structure de données de Kiar: xxx

... au cas où vous voudriez faire un coup rapide dessus.

0 commentaires
1
votes

J'ai regardé un problème similaire. J'ai trouvé référence à des intervalles de confiance calclatants de Burdick et de Graybill (Burdick, R. et Graybill, F. 1992, Intervalles de confiance sur les composants de la variance, Press CRC)

Utilisation du code que j'ai essayé de recevoir ces valeurs < / p> xxx

ici, l'intervalle de confiance plus faible pour le CV exécuté est inférieur à zéro, donc rapporté comme nan.

J'adorerais avoir une meilleure façon de Faites cela. Si je reçois le temps d'essayer de créer une fonction pour le faire.

Paul.

-

Edit: J'ai finalement écrit une fonction, ici c'est (cavalier vide) xxx

1 commentaires

Pourriez-vous élaborer sur la formule que vous utilisez pour la variance? Cherchait cela mais trouvé des informations rares

Articles qui pourrait vous intéresser :

Comment créer une barre de couleur dégradée discrète avec geom_contour_filled?
Variable de bande dans R
Mélanger le vecteur dans R, mais les éléments identiques doivent avoir une distance minimale
Pour imprimer la valeur de l'année bissextile suivante en utilisant la variable dans la programmation R