Usa questo tag per qualsiasi domanda * sull'argomento * che (a) coinvolga `R` come parte critica della domanda o risposta prevista, e (b) non è * solo * su come usare` R`.
Ho adattato un modello con diverse variabili indipendenti, una delle quali è il ritardo della variabile dipendente, usando il pacchetto dynlm. Supponendo che io abbia previsioni di un passo avanti per le mie variabili indipendenti, come posso ottenere previsioni di un passo avanti per le mie variabili dipendenti? Ecco un …
Ecco cosa voglio fare, ma non sembra esserci alcun predictmetodo per il mlogit. Qualche idea? library(mlogit) data("Fishing", package = "mlogit") Fish <- mlogit.data(Fishing, varying = c(2:9), shape = "wide", choice = "mode") Fish_fit<-Fish[-1,] Fish_test<-Fish[1,] m <- mlogit(mode ~price+ catch | income, data = Fish_fit) predict(m,newdata=Fish_test)
Diciamo che ho dei dati: x1 <- rnorm(100,2,10) x2 <- rnorm(100,2,10) y <- x1+x2+x1*x2+rnorm(100,1,2) dat <- data.frame(y=y,x1=x1,x2=x2) res <- lm(y~x1*x2,data=dat) summary(res) Voglio tracciare il continuo per interazione continua in modo tale che x1 si trovi sull'asse X e x2 sia rappresentato da 3 linee, una che rappresenta x2 con un …
Chiunque segua il baseball probabilmente ha sentito parlare delle prestazioni fuori dal nulla del tipo MVP del Jose Bautista di Toronto. Nei quattro anni precedenti, ha segnato circa 15 corse in casa per stagione. L'anno scorso ha raggiunto 54, un numero superato da soli 12 giocatori nella storia del baseball. …
Supponendo quindi che sia utile testare l'assunto di normalità per anova (vedere 1 e 2 ) Come può essere testato in R? Mi aspetterei di fare qualcosa del tipo: ## From Venables and Ripley (2002) p.165. utils::data(npk, package="MASS") npk.aovE <- aov(yield ~ N*P*K + Error(block), npk) residuals(npk.aovE) qqnorm(residuals(npk.aov)) Il che …
Sto cercando di spiegare (visivamente) la semplice correlazione lineare agli studenti del primo anno. Il modo classico di visualizzare sarebbe quello di dare un diagramma a dispersione Y ~ X con una linea di regressione diritta. Di recente, mi è venuta l'idea di estendere questo tipo di grafica aggiungendo alla …
Voglio implementare (in R) seguente dinamica molto semplice modello lineare per cui Ho 2 variabili nel tempo parametri incogniti (la varianza dell'errore di osservazione e la varianza dell'errore stato ε 2 t ).ϵ1tϵt1\epsilon^1_tϵ2tϵt2\epsilon^2_t Ytθt+1==θt+ϵ1tθt+ϵ2tYt=θt+εt1θt+1=θt+εt2 \begin{matrix} Y_t & = & \theta_t + \epsilon^1_t\\ \theta_{t+1} & = & \theta_{t}+\epsilon^2_t \end{matrix} Voglio stimare …
Ho due anni di dati che assomigliano sostanzialmente a questo: Data _ __ Violenza S / N? _ Numero di pazienti 1/1/2008 _ ___ 0 __ _ __ _ ____ 11 2/1/2008 _ __ _ 0 _ __ _ __ _ __ 11 3/1/2008 _ ____ 1 __ _ __ …
Codice di esempio: (pc.cr <- princomp(USArrests)) summary(pc.cr) loadings(pc.cr) ## note that blank entries are small but not zero Ricevo output diversi da ciascuno e non sono sicuro di capire quale sia la differenza. Ecco l'output: > summary(pc.cr) Importance of components: Comp.1 Comp.2 Comp.3 Comp.4 Standard deviation 82.8908472 14.06956001 6.424204055 2.4578367034 …
Sto cercando correlazioni tra le risposte a diverse domande in un sondaggio ("umm, vediamo se le risposte alla domanda 11 sono correlate a quelle della domanda 78"). Tutte le risposte sono categoriche (la maggior parte di esse va da "molto infelice" a "molto felice"), ma alcune hanno un diverso insieme …
Chiuso. Questa domanda è fuori tema . Al momento non accetta risposte. Vuoi migliorare questa domanda? Aggiorna la domanda in modo che sia in argomento per Cross Validated. Chiuso 2 anni fa . Voglio calcolare un riepilogo di una variabile in un data.frame per ogni combinazione unica di fattori nel …
A rischio di rendere la domanda specifica per il software, e con la scusa della sua ubiquità e idiosincrasie, voglio fare una domanda sulla funzione biplot()in R e, più specificamente, sul calcolo e sulla stampa delle sue frecce rosse sovrapposte predefinite, corrispondenti alle variabili sottostanti. [Per dare un senso ad …
Il mgcvpacchetto per Rha due funzioni per adattare le interazioni del prodotto tensore: te()e ti(). Comprendo la divisione di base del lavoro tra i due (adattamento di un'interazione non lineare rispetto alla scomposizione di questa interazione in effetti principali e un'interazione). Quello che non capisco è perché te(x1, x2)e ti(x1) …
Mi chiedo come posso simulare un campione di n vite di distribuzione Weibull che includano osservazioni censurate a destra di tipo I. Ad esempio, consente di avere n = 3, forma = 3, scala = 1 e il tasso di censura = .15 e il tempo di censura = .88. …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.