Problemi con Regressione Lineare
« il: Marzo 27, 2018, 22:07 »
Salve a tutti, sono nuovo qui e vi faccio vedere questo problema che non riesco a risolvere:

https://imgur.com/5bik8Lu   <------- immagine grafico

Le linee di regressione come vedete ( ne ho messe 2 ma dovrebbe essercene una terza ) , sono corrette ma hanno una posizione rispetto all'asse y sballata. Non riesco a capire il perchè. Il codice ( ho tolto le parti inutili )

Citazione
import numpy as np
import matplotlib.pyplot as plt

dat15 = np.genfromtxt('misure15.txt')   #plot misure 15 gradi
t15 = dat15[:,0]
v15 = dat15[:,1]
ev15 = dat15[:,2]
plt.errorbar(t15, v15, yerr=ev15, fmt='or', label = 'err(t) 15', ecolor='blue')

dat30 = np.genfromtxt('misure30.txt')   #plot misure 30 gradi   
t30 = dat30[:,0]
v30 = dat30[:,1]
ev30 = dat30[:,2]
plt.errorbar(t30, v30, yerr=ev30, fmt='or', label = 'err(t) 30', ecolor='red')

dat45 = np.genfromtxt('misure45.txt')   #plot misure 45 gradi
t45 = dat45[:,0]
v45 = dat45[:,1]
ev45 = dat45[:,2]
plt.errorbar(t45, v45, yerr=ev45, fmt='or', label = 'err(t) 45', ecolor='green')

N = len(t15)    #( = len(t30) = len(t45) )

sumt15 = np.sum(t15)    #regressione lineare 15 gradi
sumv15 = np.sum(v15)
sum2t15 = np.sum(t15*t15)
sumtv15 = np.sum(t15*v15)
d15 = N*sum2t15 - sumt15*sumt15
A15 = (sum2t15*sumv15 - sumt15*sumv15)/d15
B15 = (N*sumtv15 - sumt15*sumv15)/d15
Y15 = A15 + B15*t15
plt.plot(t15,Y15)

sumt30 = np.sum(t30)    #regressione lineare 30 gradi
sumv30 = np.sum(v30)
sum2t30 = np.sum(t30*t30)
sumtv30 = np.sum(t30*v30)
d30 = N*sum2t30 - sumt30*sumt30
A30 = (sum2t30*sumv30 - sumt30*sumv30)/d30
B30 = (N*sumtv30 - sumt30*sumv30)/d30
Y30 = A30 + B30*t30
plt.plot(t30,Y30)

Grazie a chi potrà aiutarmi :)