Deje que iff forme un triángulo. Entonces y cada . Esto es lo que ha usado para calcular el valor esperado.{ i , j , k } X = ∑ i , j , k Y i j k Y i j k ∼ B e r n o u l l i ( p 3 )Yi j k= 1{ i , j , k }X= ∑i , j , kYi j kYi j k∼ B e r n o u l l i ( p3)
Para la variación, el problema es que no son independientes. De hecho, escriba
Necesitamos calcular , que es la probabilidad de que ambos triángulos estén presentes. Hay varios casos: X 2 = ∑ i , j , k ∑ i ′ , j ′ , k ′ Y i j k Y i ′ j ′ k ′ . E [ Y i j k Y i ′ j ′ k ′ ]Yi j k
X2= ∑i , j , k∑yo′, j′, k′Yi j kYyo′j′k′.
mi[YijkYi′j′k′]
- Si (los mismos 3 vértices) entonces . Habrá tales términos en la suma doble.E [ Y i j k Y i ′ j ′ k ′ ] = p 3 ( n{ i , j , k } = { i′, j′, k′}mi[ Yi j kYyo′j′k′] = p3( n3)
- Si los conjuntos y tienen exactamente 2 elementos en común, entonces necesitamos 5 aristas presentes para obtener los dos triángulos, de modo que . habrá tales términos en la suma.{ i ′ , j ′ , k ′ } E [ Y i j k Y i ′ j ′ k ′ ] = p 5 12 ( n{ i , j , k }{ i′, j′, k′}mi[ Yi j kYyo′j′k′] = p5 512(n4)
- Si los conjuntos y tienen 1 elemento en común, entonces necesitamos 6 aristas presentes, de modo que . Habrá tales términos en la suma.{ i ′ , j ′ , k ′ } E [ Y i j k Y i ′ j ′ k ′ ] = p 6 30 ( n{i,j,k}{i′,j′,k′}E[YijkYi′j′k′]=p630(n5)
- Si los conjuntos y tienen 0 elementos en común, entonces necesitamos 6 aristas presentes, de modo que . Habrá tales términos en la suma.{ i ′ , j ′ , k ′ } E [ Y i j k Y i ′ j ′ k ′ ] = p 6 20 ( n{i,j,k}{i′,j′,k′}E[YijkYi′j′k′]=p620(n6)
Para verificar que hemos cubierto todos los casos, tenga en cuenta que la suma se suma a .(n3)2
(n3)+12(n4)+30(n5)+20(n6)=(n3)2
Recordando restar el cuadrado de la media esperada, poner todo esto junto da:
E[X2]−E[X]2=(n3)p3+12(n4)p5+30(n5)p6+20(n6)p6−(n3)2p6
Usando los mismos valores numéricos que su ejemplo, el siguiente código R calcula la desviación estándar, que es razonablemente cercana al valor de 262 de su simulación.
n=50
p=0.6
sqrt(choose(n, 3)*p^3+choose(n, 2)*(n-2)*(n-3)*p^5+(choose(n, 3)*choose(n-3, 3)+n*choose(n-1, 2)*choose(n-3, 2))*p^6-4233.6^2)
298.7945
El siguiente código de Mathematica también calcula la desviación estándar, que da el mismo resultado.
mySTD[n_,p_]:=Sqrt[Binomial[n,3]p^3+12Binomial[n,4]p^5+30 Binomial[n,5]p^6+20Binomial[n,6]p^6-(Binomial[n,3]p^3)^2]
mySTD[50,0.6] // gives 298.795