2016-09-26 2 views
0

Comment adapter le polynôme fractionnaire multivariée de la forme suivante:ajustement polynomial multivariée en R

fonction donnée: y = f(x,z), plus précisément

y = (x^2 + x^3)/(z^2 + z^3) 

Numérateur est un polynôme d'un 2e degré d'un prédicteur x , et le dénominateur est aussi un polynôme d'un second degré d'un prédicteur z.

Comment installer une telle fonction?

y = c(-4.10594369806545, -4.23691458506868, -4.24690667936422, -3.53677470254628, -4.30406509320417, -4.19442802077908, -4.66857169733859, -2.82271310942235, -4.19720194766181, 3.52164353473802, -4.3917019001973, -5.41654474791269, 2.87471821731616, -3.85922481986118, -4.25370811223789, -3.57887855889961, -5.33913936106829, -4.11775265312012, -2.89958841300109, -4.18661983833127) 

x = c(8.06526520889773, 9.39897529082673,9.07348918922699,7.5522372875608, 9.17294998275762,5.77455154554441, 9.2005930205213, 8.07309119969315, 7.42177579364465,8.18896686364888, 8.07868822922987, 8.50956416425175,9.71269017726113, 7.98378106897745, 7.69893619981345, 8.49576524400262, 8.02224091680654,8.25400859056484, 7.58171964012531, 8.35655484545343) 

z = c(2.56494935746154, 4.99043258677874, 4.43081679884331,3.66356164612965,4.53259949315326,1.79175946922805,4.23410650459726, 5.38449506278909,3.13549421592915,4.34380542185368, 3.43398720448515,2.77258872223978,6.94985645500077,3.97029191355212, 3.40119738166216,4.39444915467244,2.19722457733622,3.91202300542815,4.06044301054642, 3.87120101090789) 

df = data.frame(cbind(y=y,x=x,z=z)) 
+0

Merci. J'essaie maintenant, je me demande, cependant, comment initialiser les valeurs de départ? – aza07

+0

Je ne peux le faire en 30 '. Cette question je vais supprimer maintenant. – aza07

+0

Désolé pour la confusion. Je vais poster une nouvelle question avec une meilleure explication de ce qui doit être ajusté. – aza07

Répondre

1

Afin d'ajuster une régression des moindres carrés non-linéaire, vous pouvez utiliser nls