estou exercitando um pouco com as listas e tentando criar operações estatísticas com elas e estou com dificuldades para calcular a covariância
Basicamente estou tentando criar elementos separados e realizar o cálculo como um produto entre dois somatórios de mesmo índice i.
Basicamente, uma lista para as amostras X, ump para as amostras Y Uma para calcular o total somado do comprimento de cada lista (O somatório da lista de X com a lista Y dá o total N). Iniciei os contadores i em 1 para equivaler a operação (n-1) E basicamente tentei criar o desvio padrão para X e o desvio padrão para Y, divido-lo pela soma entre o comprimento das listas X e Y e depois exibir a covariância como um resultado numérico
Esse é meu código
#Covariância
#Índice referente ao número da amostra.
n=int(input("Insira o número total de amostras"))
#Contador para percorrer os índices
i=1 #Amostras X e Amostras Y
Amostrax=[]
Amostray=[]
#Tamanho total da lista de cada variável
for i in range(n):
#Adiciona o valor da amostra X dentro da lista de amostras
Amostrasx.append(int(input("Insira o valor das amostras associadas a X")))
#Escreve a lista de amostras X na tela para o usuário
print(Amostrasx)
#Soma das variáveis para gerar as médias de X e Y
Xi = 0
Nx = len(Amostrax)
j=1
for j in range(n):
#Adiciona o valor da amostra X dentro da lista de amostras
Amostrasy.append(int(input("Insira o valor das amostras associadas a Y")))
#Escreve a lista de amostras X na tela para o usuário
print(Amostrasy)
Yi = 0
Ny = len(Amostray)
#Preenchendo as listas de amostras
for amostraix in Amostrax:
Xi +=amostraix
for amostraiy in Amostray:
Yi +=amostraiy
Médiax=Xi/Nx
Médiay=Yi/Ny
print ("Esse é o valor da soma das amostras", Xi)
print ("Esse é o valor da soma das amostras", Yi)
print ("Este é a quantidade de amostras X que você possui", Nx)
print ("Este é a quantidade de amostras Y que você possui", Ny)
print ("Esta é a média da variável X", Médiax)
print ("Esta é a média da variável Y", Médiay)
DesviosMédiosx=[]
DesviosMédiosy=[]
for i in range (n):
DesviosMédiosx.append(Amostrax[i]-Médiax)
XiMenosXmédio=0
for desviosmediosi in DesviosMédiosx:
XiMenosXmédio+=DesviosMédios[i]
for j in range (n):
DesviosMédiosy.append(Amostray[j]-Médiay)
YiMenosXmédio=0
for desviosmediosi in DesviosMédiosy:
YiMenosYmédio+=DesviosMédios[j]
Covariância=(XiMenosXmédio*YiMenosXmédio)/(Nx+Ny)
print ("Este é cada desvio médio individual:", DesviosMédiosx)
print ("Este é cada desvio médio individual:", DesviosMédiosy)
print ("Este é o desvio médio:",Covariância)
O problema foi solucionado pelo Elton Nunes no post abaixo, só queria editar para anexar uma solução alternativa utilizando uma outra fórmula que é matematicamente equivalente e que consegui desenvolver (Referências - Estimativas e Erros experimentais em Física - Editora UERJ, 3 edição, página 26)
#Covariância
#Índice referente ao espaço amostral associadas as variáveis X e Y.
n=int(input("Insira o número total de amostras"))
#Contador para percorrer os índices
i=1 #Amostras X e Amostras Y
Amostrax=[]
Amostray=[]
Amostraxy=[]
#Tamanho total da lista de cada variável
for i in range(n):
#Adiciona o valor da amostra X dentro da lista de amostras
Amostrax.append(int(input("Insira o valor das amostras associadas a X")))
#Escreve a lista de amostras X na tela para o usuário
print("Estas são as amostras da variável X:", Amostrax)
#Soma das variáveis para gerar as médias de X e Y
Xi = 0
Nx = len(Amostrax)
j=1
for j in range(n):
#Adiciona o valor da amostra X dentro da lista de amostras
Amostray.append(int(input("Insira o valor das amostras associadas a Y")))
#Escreve a lista de amostras X na tela para o usuário
print(Amostray)
Yi = 0
Ny = len(Amostray)
#Preenchendo as listas de amostras
for amostraix in Amostrax:
Xi +=amostraix
for amostraiy in Amostray:
Yi +=amostraiy
Médiax=Xi/Nx
Médiay=Yi/Ny
for k in range(n):
Amostraxy.append(Amostrax[k]*Amostray[k])
Nxy=len(Amostraxy)
XiYi=0
for amostraixy in Amostraxy:
XiYi +=amostraixy
Médiaxy=XiYi/Nxy
Covariância=Médiaxy-(Médiax * Médiay)
print ("Esse é o valor da soma das amostras X", Xi)
print ("Esse é o valor da soma das amostras Y", Yi)
print ("Este é o produto das amostras X e Y:", Amostraxy)
print ("Este é a quantidade de amostras X que você possui", Nx)
print ("Este é a quantidade de amostras Y que você possui", Ny)
print ("Esta é a média da variável X", Médiax)
print ("Esta é a média da variável Y", Médiay)
print ("Esta é a média do produto X com Y:", Médiaxy)
print ("Esta é a covariância:", Covariância)