در این تحقیق، از 701212 رکورد روزآزمون تولید شیر از دوره شیردهی اول، 199903 رأس گاو شیری هلشتاین که از سال 1385 تا 1389 توسط مرکز اصلاح نژاد کرج جمع­آوری شده بود، استفاده شد. مدل­ رگرسیون تصادفی با توابع چند جمله­ای ­ های لژاندر با درجات مختلف 2 تا 5 برای برازش اثرات ژنتیک افزایشی و محیط دائمی، ویلمینک و علی- شیفر تحت واریانس باقی­مانده ثابت در کل دوره شیردهی مقایسه شدند. پارامترهای ژنتیکی با استفاده از روش حداکثر درستنمائی محدود شده (REML) ­ ­، برآورد شدند. مقایسه مدل­ها با استفاده از -2Logl ، معیار اطلاعات آکایک ( AIC )، معیار اطلاعات بیزی ( BIC )، واریانس باقیمانده ( RV ) و آزمون نسبت درستنمائی ( LRT ) انجام گرفت. با توجه به نتایج بدست آمده، مدل رگرسیون تصادفی با چندجمله­ای لژاندر (2،5) به عنوان بهترین مدل انتخاب شد . واریانس باقی­مانده همراه با افزایش درج ه برازش اثرات محیط دائمی در چند جمله ­ای­های لژاندر، کاهش یافت. واریانس محیط دائمی در اوایل دوره شیردهی بالاتر از دیگر مراحل دوره شیردهی و واریانس ژنتیک افزایشی در اوایل دوره نسبت به انتهای دوره شیردهی پایین­تر برآورد شد. واریانس فنوتیپی تولید شیر در طول دوره شیردهی ثابت نبود و در ابتدا و انتهای دوره شیردهی بالاتر بود. دامنه وراثت­پذیری در طول دوره شیردهی در بین توابع مختلف بین 08/0 تا 23/0 متغیر بود .