2009-11-23 9 views
6

Jestem programistą w Portland, OR. Zastanawiam się, czy ktoś może pomóc:Loess Fit and Resulting Equation

pracuję nad modeli nadających lessowe używając R, raz mam pasowanie osiągnięte, szukam do back-out równanie dopasowanej krzywej nieliniowej , zastanawiasz się, czy istnieje sposób, aby ustalić to równanie w R? Szukałem, ale nie mogę znaleźć żadnej literatury. Dla mnie wykres funkcji jest świetny, ale bez równania wykresu, jestem trochę martwy w wodzie.

Odpowiedz

1

Chcesz odzyskać formułę z obiektu less? Możesz być w stanie zrobić coś takiego:

> cars.lo <- loess(dist ~ speed, cars) 
> formula(unclass(cars.lo)$terms) 
dist ~ speed 

Edit: Niestety ... myślę, że błędnie interpretowane, co chciał. Nie ma prostego sposobu na wyrażenie modelu lessowego w postaci równania.

+0

Dzięki, sprawdzam to i sprawdzam, czy to działa! –

7

Loess nie daje równania [1]. Jeśli chcesz po prostu uzyskać wartości zwracanych przez funkcję lessowej użyć predict(loess.object, new.data)

[1] Z Wikipedii:

Inną wadą lessowe pada fakt, że nie wywołuje funkcję regresji że jest łatwo reprezentowany przez formułę matematyczną. Może to utrudnić przeniesienie wyników analizy do innych osób . Aby przekazać funkcję regresji innej osobie, , potrzebowaliby oni zestawu danych i oprogramowania do obliczeń LOESS.

+0

poprawne, nie daje wzoru matematycznego, ale daje wykres. Chcę znaleźć formułę dla podanego wykresu. –

+0

Nie możesz. Loess stosuje iteracyjne podejście lokalne i jest oparte na danych, zobacz dokumentację referencyjną, aby uzyskać szczegółowe informacje. –

7

Nie ma formuły. Less to metoda nieparametryczna. Nie można go wyrazić jako prostego równania.

1

Nie można wyodrębnić formuły z samego LOESS. Możesz jednak po prostu uruchomić inną metodę na punktach znalezionych przez LOESS. Jeśli jest to prosty wykres 2D, nie powinno być tak trudno znaleźć dobrą formułę. Jedną z metod jest regresja symboliczna (patrz: wiki).

Należy pamiętać, że prawdopodobnie nie jest to optymalne i lepiej zastosować inną metodę niż LOESS.