2010-04-19 9 views
2

Comment puis-je extraire les paramètres que la fonction loess a ajustés pour la fonction polynomiale qu'elle utilise, pour un x particulier de ma Les données?Extraction des termes ajustés dans la fonction polynomiale locale d'un loess (dans R) (NOT predict())

Par exemple, dans:

cars.lo <- loess(dist ~ speed, cars) 
cars.lo 

Qu'est-ce qu'il convient quand cars.lo $ x == 5?

Mise à jour: Je veux les paramètres de la fonction polynomiale, pas la prédiction (prédire) du loess.

Je lui demande d'obtenir une estimation de la pente à ce point.

Merci, Tal

Répondre

1

Je dirais

predict(cars.lo, data.frame(speed=5)) 
[1] 7.797353 
+0

oops - Je vais reformuler - je veux les paramètres polynomiaux ** montés **, pas la valeur ajustée. –

+0

Avez-vous essayé 'str (summary (cars.lo))'? Pouvez-vous obtenir de l'aide là-bas? –

+0

Je l'ai essayé, c'est pourquoi je demande ici :) (la sortie est complexe, et je préfère avoir l'aide de quelqu'un). Merci de toute façon :) –

2

Je calculerait les valeurs prévues à x, x + EPS, x-EPS, puis adapter à un second degré aux résultats. Ce n'est pas très efficace en temps d'ordinateur, mais si vous n'avez pas à le faire très souvent, c'est très efficace en temps de programmation.

+0

Merci Rob, je pensais que de le faire de cette façon, mais un statisticien principal en milieu de travail a dit qu'il causerait des résultats instables (par rapport à l'ajustement polynomiale qui repose sur de nombreux points de données autour X0) ... –

0

Je pense que vous allez être seul pour celui-ci. Je commence par lire les références dans lowess, puis lire la source de la fonction C lowess. Je recommande de commencer là au lieu de loess parce que c'est un algorithme plus ancien et légèrement plus simple que vous pourriez trouver plus facile à adapter à vos besoins.

+0

Merci Hadley. Je craignais que c'était le cas. Je vais retourner chez le client et voir à quel point c'est important ... Best, Tal –

+0

Salut Tal Galili, avez-vous trouvé une solution? Si oui, pourriez-vous l'afficher? Merci. – Annamarie