ANOVA in python

Question

Мой файл csv содержит дискретные и непрерывные переменные, и я хотел бы найти «уравнение модели», которое объясняет мою непрерывную переменную (a) в соответствии с моими дискретными переменными (x, y, z); -> a = f (x, y, z). Проблема в том, что я пытаюсь использовать этот код, но он терпит неудачу. Когда я хочу увидеть результат печати ([f_value, p_value] stats.f_oneway = (x, y, z)), я получаю [nan, nan]. Выписка из моего кода:

from numpy import (genfromtxt,hstack,arange) 
#Pr linear regression 
from scipy import stats 
import scipy 
#Pr ANOVA 
from statsmodels.stats.multicomp import (pairwise_tukeyhsd,MultiComparison) 

from pylab import savefig 
from matplotlib.pyplot import (figure,setp) 

fname="G:/table.csv" 
my_data = genfromtxt(fname,delimiter=',') 

#Transformation of file into table 

x= my_data[:,3] 
y= my_data[:,4] 
z= my_data[:,6] 


#one way anova 
[f_value, p_value] = stats.f_oneway(x, y, z) 

Я хочу оценить модель R = U + F (х, у, г) + е, где R (непрерывная переменная), и (константа), е (погрешность измерения). Я хочу иметь коэффициенты для того, как дискретные переменные влияют на мой непрерывный результат.

Answer 1

Мое лучшее предположение заключается в том, что ваши данные содержат значения NaN или плохие значения. Вы могли бы попытаться обнаружить это с следующим фрагментом кода:

for i in range(1,len(x)): 
    (f_val,p_val) = stats.f_oneway(x[:i],y[:i],z[:i]) 
    if numpy.isnan(f_val) or numpy.isnan(p_val): 
     print i-1,x[i-1],y[i-1],z[i-1],f_val,p_val