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`.
Il modello Bradley – Terry – Luce (BTL) afferma che , dove è la probabilità che l'oggetto sia giudicato "migliore", più pesanti, ecc., dell'oggetto , e e sono parametri.p i j j i δ i δ jpji=logit−1(δj−δi)pji=logit−1(δj−δi)p_{ji} = logit^{-1}(\delta_j - \delta_i)pijpijp_{ij}jjjiiiδiδi\delta_iδjδj\delta_j Questo sembra essere un candidato per la funzione glm, …
I miei dati sono una serie temporale di popolazione occupata, L e periodo di tempo, anno. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e. NaN NaN NaN NaN 1.6058 0.0008 sigma^2 estimated as 1.503e-06: log likelihood=107.55 …
So che questa è probabilmente una domanda di base ... Ma non sembra trovare la risposta. Sto adattando un GLM con una famiglia Poisson, quindi ho cercato di dare un'occhiata alle previsioni, tuttavia l'offset sembra essere preso in considerazione: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict (model_glm, type="response") Ottengo casi non …
@whuber ha dimostrato come simulare esiti multivariati ( , y 2 e y 3 ) per un punto temporale.y1y1y_1y2y2y_2y3y3y_3 Come sappiamo, i dati longitudinali si verificano spesso negli studi medici. La mia domanda è come simulare esiti multivariati di misure ripetute in R? Ad esempio, misuriamo ripetutamente , y 2 …
Ho bisogno di un pacchetto che possa darmi l'equazione per un modello SVM lineare. Attualmente sto usando e1071 in questo modo: library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho #Offset Tuttavia, non sono sicuro di come e1071::svm()selezioni le …
Ho problemi nell'uso delle funzioni cor()e cor.test(). Ho solo due matrici (solo valori numerici e lo stesso numero di righe e colonne) e voglio avere il numero di correlazione e il corrispondente valore p. Quando uso cor(matrix1, matrix2)ottengo i coefficienti di correlazione per tutte le celle. Voglio solo un singolo …
Sto usando il pacchetto kernlab in R per creare un SVM per classificare alcuni dati. SVM sta funzionando bene in quanto fornisce "previsioni" di un'accuratezza decente, tuttavia il mio elenco di variabili di input è più grande di quello che vorrei e non sono sicuro dell'importanza relativa delle diverse variabili. …
Supponi di avere dati di sopravvivenza come questo: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Per eseguire un test standard del rango di registro, è possibile utilizzare survdiff(Surv(time, status) ~ group, data = obs, …
Come hobby secondario, ho esplorato le serie temporali di previsione (in particolare, usando R). Per i miei dati, ho il numero di visite al giorno, per ogni giorno che risale a quasi 4 anni fa. In questi dati ci sono alcuni schemi distinti: Dal lunedì al venerdì ci sono molte …
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 3 anni fa . Sono molto interessato a imparare come creare report ricorrenti dal mio codice R e dalla visualizzazione ggplot2. …
Sto cercando un semplice esempio di codice su come eseguire un filtro antiparticolato in R. Il pacchetto pomp sembra supportare il bit matematico dello spazio degli stati, ma gli esempi sono un po 'complicati da seguire programmaticamente per un semplice sviluppatore OO come me, in particolare come caricare i dati …
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 9 mesi fa . Il problema menzionato in questa domanda è stato risolto nella versione 1.7.3 del pacchetto R glmnet. Sto …
Ho creato la mia versione leggermente migliorata del termplot che uso in questo esempio, puoi trovarla qui . Ho precedentemente pubblicato su SO, ma più ci penso e credo che ciò sia probabilmente più correlato all'interpretazione del modello di rischi proporzionali di Cox che alla codifica effettiva. Il problema Quando …
Ho alcuni dati di serie temporali in cui la variabile misurata è numeri interi positivi discreti (conteggi). Voglio testare se c'è una tendenza al rialzo nel tempo (o meno). La variabile indipendente (x) è nell'intervallo 0-500 e la variabile dipendente (y) è nell'intervallo 0-8. Ho pensato di rispondere a questo …
Se ho una variabile con 4 livelli, in teoria devo usare 3 variabili fittizie. In pratica, come viene effettivamente realizzato? Uso 0-3, utilizzo 1-3 e lascio vuoto il 4? Eventuali suggerimenti? NOTA: lavorerò in R. AGGIORNAMENTO: Cosa accadrebbe se uso solo una colonna che utilizza 1-4 corrispondente a AD? Funzionerà …
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.