import sys
import pandas as pd

prod = pd.read_csv('datos_produccion.csv', index_col = 0, encoding = 'utf-8')
prod.apply(lambda x: pd.lib.infer_dtype(x.values))

lista_productos= prod[['planta', 'producto', 'velocidad de producción']]

#antes de generar un nuevo archivo deseo agrupar los datos y rellenar los faltantes
prod.to_csv('listra_produccion.csv', encoding = 'utf-8')

import pandas as pd

prod = pd.read_csv('datos.csv', index_col=0, encoding='utf-8')
groups_means = prod.groupby(["planta", "producto"])['velocidad de producción'].transform('mean')
prod.loc[prod['velocidad de producción'].isnull(), 'velocidad de producción'] = groups_means

>>> prod

         planta    producto  velocidad de producción
index
0      Planta 1  producto 1                     40.0
1      Planta 3  producto 2                     20.0
2      Planta 1  producto 1                     45.0
3      Planta 1  producto 1                     50.0
4      Planta 3  producto 2                     25.0