Deutsches R-Forum

Hallo,
ich kenne zwar stepAIC, aber gibt es eine Möglichkeit, automatisch das beste Modell zu finden, die auch Wechselwirkungen/Quadrierte Parameter und ähnliches beachtet?

Mit freundlichen Grüßen
Gwinzy

Den p-Wert habe ich jetzt doch schon alleine geschafft, aber ich hänge bei der nächsten Aufgabe.

Eine Funktion, die einen F-Test für die Nullhypothese H0: L*beta=r mit selbstgewählten L und r (bei vernünftig gewählter Voreinstellung) erzeugt.
Ich scheitere schon daran, dass ich keine Ahnung habe ...

Hallo,

ich soll den p-wert eines Regressionsmodell für alle Koeffizienten berechnene, so wie es summary macht.
Leider finde ich da keinen passenden Ansatz. Hat da jemand eine Idee?

Mit freundlichen Grüßen
Gwinzy

Danke für die Tipps, so schaut es echt gleich viel schöner aus

attach(swiss)
meinlogit <- function (x,y)
{ idx <- (1:length(x))
nenner <- 0
zaehler <- 0

beta <- 0
for(i in idx){
zaehler <- zaehler + (x[i]-mean(x))*(y[i]-mean(y))
nenner <- nenner + (x[i]-mean(x))^2

}
beta[2] <- zaehler/nenner
beta[1] <- mean(y)-beta[2]*mean(x)

return (list(beta ...

Hallo Ruedi und Jörg,
danke, aber ich habe jetzt eine andere Formel benutzt, mit der ich auf das richtige Ergebnis gekommen bin.

Hallo,

ich soll eine kleinste-Quadrate Schätzung programmieren.

Meine Funktion ist die folgende:
attach(swiss)
meinlogit <- function (x,y)
{
(solve(t(x)%*%x)[1,1]*t(x))%*%y
}
meinlogit(Education, Fertility)

Allerdings ist meine Ergebnis nicht das gleiche wie das von
summary(lm(Fertility ...