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`.
Esiste un modo in R (una funzione integrata) per calcolare la matrice di transizione per una catena di Markov da una serie di osservazioni? Ad esempio, prendere un set di dati come il seguente e calcolare la matrice di transizione del primo ordine? dat<-data.frame(replicate(20,sample(c("A", "B", "C","D"), size = 100, replace=TRUE)))
Sono abbastanza nuovo nelle statistiche e ho bisogno del tuo aiuto. Ho un piccolo campione, come segue: H4U 0.269 0.357 0.2 0.221 0.275 0.277 0.253 0.127 0.246 Ho eseguito il test Shapiro-Wilk usando R: shapiro.test(precisionH4U$H4U) e ho ottenuto il seguente risultato: W = 0.9502, p-value = 0.6921 Ora, se presumo …
Uso sempre lm()in R per eseguire la regressione lineare di su . Tale funzione restituisce un coefficiente tale chex β y = β x .yyyxxxββ\betay=βx.y=βx.y = \beta x. Oggi ho imparato a conoscere i minimi quadrati totali e quella princomp()funzione (analisi dei componenti principali, PCA) può essere utilizzata per eseguirlo. …
Uso la funzione auto.arima () nel pacchetto di previsione per adattarsi ai modelli ARMAX con una varietà di covariate. Tuttavia, ho spesso un gran numero di variabili tra cui scegliere e di solito finisco con un modello finale che funziona con un sottoinsieme di esse. Non mi piacciono le tecniche …
Possiamo usare lm()per prevedere un valore, ma in alcuni casi abbiamo ancora bisogno dell'equazione della formula del risultato. Ad esempio, aggiungere l'equazione ai grafici.
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 . Sto usando il cursore per eseguire una foresta casuale convalidata in modo incrociato su un set di …
Sto usando il pacchetto robustbase per eseguire una stima glm. Tuttavia quando lo faccio, ottengo il seguente errore: Error in solve.default(crossprod(X, DiagB * X)/nobs, EEq) : system is computationally singular: reciprocal condition number = 1.66807e-16 Cosa significa / indica? E come posso eseguire il debug? PS. Se hai bisogno di …
Sto usando i modelli R (3.1.1) e ARIMA per le previsioni. Vorrei sapere quale dovrebbe essere il parametro "frequenza", che è assegnato nella ts()funzione , se sto usando dati di serie temporali che sono: separato da minuti e si sviluppa su 180 giorni (1440 minuti / giorno) separato da secondi …
Sto usando glmnet per calcolare le stime di regressione della cresta. Ho ottenuto alcuni risultati che mi hanno fatto sospettare che glmnet stia davvero facendo quello che penso faccia. Per verificare questo ho scritto un semplice script R in cui comparo il risultato della regressione della cresta effettuata da risolvere …
[Il titolo iniziale "Misurazione della somiglianza per gli alberi del cluster gerarchico" è stato successivamente modificato da @ttnphns per riflettere meglio l'argomento] Sto eseguendo una serie di analisi di gruppi gerarchici su un frame di dati dei record dei pazienti (ad esempio simile a http://www.biomedcentral.com/1471-2105/5/126/figure/F1?highres=y ) Sto sperimentando misure di …
Vorrei ottenere un valore p e una dimensione dell'effetto di una variabile categoriale indipendente (con più livelli) - che è "complessiva" e non per ogni livello separatamente, come è l'output normale di lme4in R. È proprio come la cosa che la gente riferisce quando esegue un ANOVA. Come posso ottenerlo?
Bloccato . Questa domanda e le sue risposte sono bloccate perché la domanda è fuori tema ma ha un significato storico. Al momento non accetta nuove risposte o interazioni. Capisco che questa domanda sia piuttosto ampia, ma mi chiedo quali dovrebbero essere i punti decisivi nel decidere di creare (o …
Spero che a tutti voi non dispiaccia questa domanda, ma ho bisogno di aiuto per interpretare l'output per un output del modello a effetti misti lineari. Ho cercato di imparare a fare in R. Sono nuovo nell'analisi dei dati longitudinali e nella regressione lineare degli effetti misti. Ho un modello …
Ad esempio, considera il ChickWeightset di dati in R. La varianza ovviamente aumenta nel tempo, quindi se uso una semplice regressione lineare come: m <- lm(weight ~ Time*Diet, data=ChickWeight) Le mie domande: Quali aspetti del modello saranno discutibili? I problemi si limitano all'estrapolazione al di fuori Timedell'intervallo? Quanto è tollerante …
Mi sono appena imbattuto in questo documento , che descrive come calcolare la ripetibilità ( nota anche come affidabilità, nota anche come correlazione intraclasse) di una misurazione tramite la modellazione di effetti misti. Il codice R sarebbe: #fit the model fit = lmer(dv~(1|unit),data=my_data) #obtain the variance estimates vc = VarCorr(fit) …
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.