Realmente não sei o que é possível e gostaria de um ponteiro na direção certa. Tenho medições de tempo e posição que podem ser qualquer coisa de uma pessoa andando, um veículo na estrada ou em um estacionamento ou uma impressora em um escritório. Preciso calcular os tempos de viagem dos veículos...
Eu tenho o modelo que preciso estimar, com e .Y= π0 0+ π1X1+ π2X2+ π3X3+ ε ,Y=π0+π1X1+π2X2+π3X3+ε, Y = \pi_0 + \pi_1 X_1 + \pi_2 X_2 + \pi_3 X_3 + \varepsilon, π k ≥ 0 para k ≥ 1∑kπk= 1 para k ≥ 1∑kπk=1 for k≥1\sum_k \pi_k = 1 \text{ for }k \geq 1πk≥ 0 para k ≥ 1πk≥0 for k≥1\pi_k\ge0 \text{...
Uma rápida olhada na página de ajuda do gerador de números aleatórios de shows R que você pode escolher entre 7 geradores predefinidos ( Wichmann-Hill, Marsaglia-Multicarry, Super-Duper, Mersenne-Twister, Knuth-TAOCP-2002, Knuth-TAOCP, L'Ecuyer-CMRG). ?Random O padrão é Mersenne-Twister , que...
Estou implementando um classificador SVM não linear com o kernel RBF. Disseram-me que a única diferença em relação a um SVM normal era que eu tinha que simplesmente substituir o produto do ponto por uma função do kernel: Sei como funciona um SVM linear normal, ou seja, depois de resolver o...
Eu tenho dois grupos de 10 participantes que foram avaliados três vezes durante um experimento. Para testar as diferenças entre os grupos e nas três avaliações, executei um ANOVA de desenho misto 2x3 com group(controle, experimental), time(primeiro, segundo, três) e group x time. Ambos timee...
Consegui criar uma rede neural dos meus dados. Mas não tenho tanta certeza sobre a interpretação da saída R. Eu usei o seguinte comando para criar rede neural: > net=nnet(formula = category~iplen+date_time, size=0,skip=T,lineout=T) # weights: 3 initial value 136242.000000 final value...
Eu medi dois parâmetros (carbono orgânico dissolvido DOC = y, e descarga = x). Quando essas duas variáveis são plotadas uma contra a outra, obtemos um loop de histerese (veja o exemplo de código e a figura). Agora, para uma análise mais aprofundada, eu gostaria de determinar a área desse loop...
Desejo calcular o parâmetro da distribuição exponencial de uma população de amostra retirada dessa distribuição sob condições tendenciosas. Até onde eu sei, para uma amostra de n valores, o estimador usual é . No entanto, minha amostra é tendenciosa da seguinte forma:e - λ x λ = nλλ\lambdae- λ...
Meu conjunto de dados possui as seguintes variáveis: Tratamento (4 tipos - fixo) Localização (8 localizações - fixas) Posição no local (3 posições por local - fixo) As amostras são coletadas em cada posição (3 amostras por posição aleatória) Tempo (dois tempos de amostragem - fixos) Taxa de...
Pesquisei o pacote de mouses e ainda não descobri uma maneira de usar as várias imputações para criar um modelo de Cox e validá-lo com a função do pacote rmsvalidate() . Aqui está um código de amostra do que tenho até agora, usando o conjunto de dados
Como acompanhamento da sugestão de modelo para uma regressão de Cox com covariáveis dependentes do tempo, aqui está o gráfico de Kaplan Meier, responsável pela natureza dependente do tempo das gestações. Em outras palavras, o conjunto de dados agora é dividido em um longo conjunto de dados com...
Alguém tem um breve exemplo educacional rápido de como usar redes neurais ( nnetem R, por exemplo) para fins de previsão? Aqui está um exemplo, em R, de uma série temporal T <- seq(0,20,length=200) Y <- 1 + 3*cos(4*T+2) +.2*T^2 + rnorm(200) plot(T,Y,type="l") Este é apenas um exemplo, mas...
Atualmente, estou usando R para prever uma série temporal com estas instruções: X <- ts(datas, frequency=24) X.arima <- Arima(X, order=c(2,1,0), seasonal=c(1,1,1)) pred <- predict(X.arima, n.ahead=24) plot.ts(pred$pred) Como você pode ver, tenho dados a cada hora e escolhi o período...
Estou tentando executar uma análise de crescimento de classe latente (LCGA) e / ou modelos de mistura de crescimento (GMMs) em R. Os dados que estou usando são um número crescente de bifurcações de repositórios git (variável discreta, não categórica), como você pode veja neste conjunto de dados...
Os dados consistem em espectros ópticos (intensidade da luz em relação à frequência) gravados em momentos variados. Os pontos foram adquiridos em uma grade regular em x (tempo), y (frequência). Para analisar a evolução do tempo em frequências específicas (um aumento rápido, seguido de uma...
Tenho problemas para encontrar uma solução sobre como executar um teste post-hoc (Tukey HSD) após uma ANOVA de medidas repetidas de dois fatores (ambos em indivíduos) em R. Para a ANOVA, usei a função aov: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1)) Depois de ler as respostas...
Ao lidar com dados com fatores R, pode ser usado para calcular as médias de cada grupo com a função lm (). Isso também fornece os erros padrão para as médias estimadas. Mas esse erro padrão difere do que recebo de um cálculo manualmente. Aqui está um exemplo (retirado daqui Prevendo a diferença...
Estou usando o lme4pacote no R para fazer alguma modelagem logística de efeitos mistos. Meu entendimento era que a soma de cada efeito aleatório deveria ser zero. Quando eu uso modelos mistos lineares de brinquedo lmer, os efeitos aleatórios geralmente são < confirmando minha crença de que...
Eu tenho dados que se parecem um pouco com isso: ID Status 01 A 02 G 03 E ... ... 100 G Você entendeu a ideia, eu acho. Eu tenho esses dados de duas populações separadas (coortes) e quero comparar a distribuição da variável status em uma população com a distribuição em outra. A pergunta que estou...
Eu tenho uma pergunta sobre como um estatístico normalmente interpretaria uma saída anova. Digamos que tenho saída anova de R. > summary(fitted_data) Call: lm(formula = V1 ~ V2) Residuals: Min 1Q Median 3Q Max -2.74004 -0.33827 0.04062 0.44064 1.22737 Coefficients: Estimate Std. Error t...