2016-12-24 46 views
0

我有日期温度数据集,按日期排序,我需要在scikit-learn中使用[SVR] [1]预测未来温度。scikit学习svr的时间序列预测

我坚持选择XY的培训和X测试 设置。例如,如果我想在时间t预测Y然后我需要 训练集包含了X & Yt-1, t-2, ..., t-N其中N是用于在t预测Y以前的天数。

我该怎么做?

这是它。

df=daily_temp1 
# define function for create N lags 
def create_lags(df, N): 
    for i in range(N): 
     df['datetime' + str(i+1)] = df.datetime.shift(i+1) 
     df['dewpoint' + str(i+1)] = df.dewpoint.shift(i+1) 
     df['humidity' + str(i+1)] = df.humidity.shift(i+1) 
     df['pressure' + str(i+1)] = df.pressure.shift(i+1) 
     df['temperature' + str(i+1)] = df.temperature.shift(i+1) 
    df['vism' + str(i+1)] = df.vism.shift(i+1) 
    df['wind_direcd' + str(i+1)] = df.wind_direcd.shift(i+1) 
    df['wind_speed' + str(i+1)] = df.wind_speed.shift(i+1) 
    df['wind_direct' + str(i+1)] = df.wind_direct.shift(i+1) 

    return df 

# create 10 lags 
df = create_lags(df,10) 


# the first 10 days will have missing values. can't use them. 
df = df.dropna() 

# create X and y 
y = df['temperature'] 
X = df.iloc[:, 9:] 

# Train on 70% of the data 
train_idx = int(len(df) * .7) 

# create train and test data 
X_train, y_train, X_test, y_test = X[:train_idx], y[:train_idx], X[train_idx:], y[train_idx:] 


# fit and predict 
clf = SVR() 
clf.fit(X_train, y_train) 

clf.predict(X_test) 

回答

1

下面是构建特征矩阵X作为LAG1简单的解决方案 - lagN其中LAG1是前几天的温度和lagN是温度N天前。

# create fake temperature 
df = pd.DataFrame({'temp':np.random.rand(500)}) 

# define function for create N lags 
def create_lags(df, N): 
    for i in range(N): 
     df['Lag' + str(i+1)] = df.temp.shift(i+1) 
    return df 

# create 10 lags 
df = create_lags(df,10) 

# the first 10 days will have missing values. can't use them. 
df = df.dropna() 

# create X and y 
y = df.temp.values 
X = df.iloc[:, 1:].values 

# Train on 70% of the data 
train_idx = int(len(df) * .7) 

# create train and test data 
X_train, y_train, X_test, y_test = X[:train_idx], y[:train_idx], X[train_idx:], y[:train_idx] 

# fit and predict 
clf = SVR() 
clf.fit(X_train, y_train) 

clf.predict(X_test) 
+0

谢谢!我试图让它符合我的目的,因为我还需要除了以前的温度之外的其他功能,但我一直有错误TypeError:float()参数必须是字符串或数字 –

+0

scikit学习不能处理字符串。您必须将所有字符串转换为数字。如果你有一个熊猫数据框,使用'pd.get_dummies'。如果您严格使用sklearn,请在预处理模块中使用“LabelBinarizer”。 –

+0

我已经改变了弦的,现在我有浮动和日期时间数据类型 –