Napisz odpowiedź

Reklama

Statystyka, prognozowanie, ekonometria, data mining
Forum miłośników statystyki - Portal Statystyczny

FAQ

Szukaj

Użytkownicy

Grupy

Statystyki

Rejestracja

Zaloguj

Chat

Regulamin

Kadra forum

PORTAL STATYSTYCZNY

Ogłoszenie
FORUM STATYSTYCZNE MA JUŻ 10 LAT Znasz statystykę lub ekonometrię, metody prognozowania, data mining i chcesz pomóc w rozwoju forum statystycznego ? Pisz na: administrator(małpa)statystycy.pl Rozpoczął swoją działalność portal statystyczny - masz pomysł na jego rozwój ? Drogi forumowiczu! Zanim napiszesz posta zapoznaj się z regulaminem forum i przedstaw się

Statystyka, prognozowanie, ekonometria, data mining Strona Główna » GRUPA UŻYTKOWNIKÓW R » Biblioteki R, Pakiety R » Problem z nlsLM

Napisz odpowiedź

Użytkownik

Temat

Ikona Tematu/Postu

Treść wiadomości

Ikony Emocji

µ	Ω	Π	φ	Δ	Θ	Λ	Σ	Φ	Ψ	α
β	χ	τ	γ	δ	ε	ζ	η	ψ	θ	λ
ξ	ρ	σ	ω	κ	Γ	♣	♥	ø	§	©
®	•	™	°	«	»	≤	≥	³	²	½
¼	¾	±	÷	×	√	∞	∫	≈	≠	≡
≈	←	→	↑	↓	↔	€	£	¥	¢	ƒ

Kolor:		Cień:		Ogień:
Rozmiar:

[quote="mikos"][code]Dane <- read.csv("dim.csv",sep=";")
head(Dane,3)
  Sales      TV     RD
1 1 222 100 000 30 000
2 1 201  70 000 20 000
3 1 149  35 000 40 000
[/code]
[code]sapply(Dane,class)
sapply(Dane,class)
   Sales       TV       RD 
"factor" "factor" "factor" 
[/code]
Nie można wykonać obliczeń na czynnikach - factorach.
Zmienne muszą być zakodowane jako numeric - zmienne numeryczne/liczby.

* likwidujemy spację w factorach:
[code]Dane02 <- as.data.frame(apply(Dane,2,function(x) gsub('\\s+', '',x)))
[/code]

* konwertujemy wszystkie zmienne typu factor na character (napisy) a następnie na numeric (liczby):
[code]for (i in names(Dane02)) Dane02[,i] <- as.numeric(as.character(Dane02[,i]))
[/code]

[code]head(Dane02,3)
  Sales     TV    RD
1  1222 100000 30000
2  1201  70000 20000
3  1149  35000 40000
[/code]

[code]sapply(Dane02,class)
    Sales        TV        RD 
"numeric" "numeric" "numeric" 
[/code]

[code]power = function(x, B) { 
x^B 
}

nls_model = nlsLM(Sales ~ b0 + a1*power(TV, b1) + a2*power(RD, b2), 
data = Dane02, 
algorithm = "LM", 
start = c(b0 = 500, b1 = 0.5, b2 = 0.5, a1 = 0.5, a2 = 0.5), 
lower = c(b0 = 200, b1 =0, b2 = 0, a1 = 0, a2 = 0), 
upper = c(b0 = 2000, b1 = 1, b2 = 1, a1 = 1, a2 = 1)) 
summary(nls_model)

Formula: Sales ~ b0 + a1 * power(TV, b1) + a2 * power(RD, b2)

Parameters:
    Estimate Std. Error t value Pr(>|t|)    
b0 9.923e+02  8.530e+01  11.633  < 2e-16 ***
b1 4.877e-01  1.653e-01   2.950  0.00355 ** 
b2 1.000e+00  8.850e+00   0.113  0.91014    
a1 1.000e+00  2.172e+00   0.460  0.64575    
a2 1.318e-04  1.313e-02   0.010  0.99200    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 30.22 on 204 degrees of freedom

Number of iterations to convergence: 8 
[/code][/quote]

Freak & Letter styles
Ctrl+Z aby cofnąć

Opcje
HTML: NIE
BBCode: TAK
Uśmieszki: TAK

	Wyłącz BBCode w tym poście
	Wyłącz Uśmieszki w tym poście

Skocz do:

Ideą przyświecającą istnieniu forum statystycznego jest stworzenie możliwości wymiany informacji, poglądów i doświadczeń osób związanych ze statystyką, mierzenie się z różnego rodzaju problemami statystycznymi i aktuarialnymi. Poruszane problemy: Statystyka w badaniach sondażowych rynku, metody reprezentacyjne, Teoria i rachunek prawdopodobieństwa, statystyka opisowa, teoria estymacji, testowanie hipotez statystycznych, ekonometria, prognozowanie, metody data mining.

Strona wygenerowana w 0,51 sekundy. Zapytań do SQL: 8