Simulação de novo de dados a partir de um quadro de dados de projeto experimental. Com foco em R (embora outra solução de idioma seja ótima). Ao projetar um experimento ou uma pesquisa, a simulação de dados e a realização de uma análise desses dados simulados podem fornecer informações excelentes...
Estou trabalhando na minha tese de mestrado no momento e planejei executar as estatísticas com o SigmaPlot. No entanto, depois de passar algum tempo com meus dados, cheguei à conclusão de que o SigmaPlot pode não ser adequado para o meu problema (posso estar enganado) e iniciei minhas primeiras...
Eu tenho usado o MCMCglmmpacote recentemente. Estou confuso com o que é referido na documentação como estrutura R e estrutura G. Eles parecem estar relacionados aos efeitos aleatórios - em particular especificando os parâmetros para a distribuição anterior sobre eles, mas a discussão na...
Estou lidando com um problema de classificação binária supervisionada. Eu gostaria de usar o pacote GBM para classificar os indivíduos como não infectados / infectados. Eu tenho 15 vezes mais não infectado do que indivíduos infectados. Eu queria saber se os modelos GBM sofrem no caso de tamanhos...
Eu tenho um conjunto completo de seqüências (432 observações para ser mais preciso) de 4 estados : por exemploA−DA−DA-D Y=⎛⎝⎜⎜⎜⎜AB⋮BCA⋮CDA⋮ADC⋮DBA⋮AA−⋮BC−⋮A⎞⎠⎟⎟⎟⎟Y=(ACDDBACBAACA−−⋮⋮⋮⋮⋮⋮⋮BCADABA)Y=\left(\begin{array}{c c c c c c c} A& C& D&D & B & A &C\\ B& A& A&C & A&-...
Na rotina rpart () para criar modelos CART, você especifica o parâmetro de complexidade para o qual deseja remover sua árvore. Eu vi duas recomendações diferentes para escolher o parâmetro de complexidade: Escolha o parâmetro de complexidade associado ao erro mínimo possível de validação cruzada....
Em seu artigo amplamente citado, distribuições anteriores para parâmetros de variância em modelos hierárquicos (916 citação até agora no Google Scholar) Gelman propõe que boas distribuições prévias não informativas para a variação em um modelo bayesiano hierárquico são a distribuição uniforme e a...
Usando Rpara executar a decomposição do STL, s.windowcontrola a rapidez com que o componente sazonal pode mudar. Valores pequenos permitem mudanças mais rápidas. Definir a janela sazonal como infinita equivale a forçar o componente sazonal a ser periódico (ou seja, idêntico ao longo dos...
Versão curta: Eu tenho uma série temporal de dados climáticos que estou testando quanto à estacionariedade. Com base em pesquisas anteriores, espero que o modelo subjacente (ou "gerador", por assim dizer) dos dados tenha um termo de interceptação e uma tendência temporal linear positiva. Para...
Acabei de assistir novamente a uma palestra do curso Machine Learning no Coursera. Na seção em que o professor discute o PCA para pré-processamento de dados em aplicativos de aprendizado supervisionado, ele diz que o PCA deve ser executado apenas nos dados de treinamento e, em seguida, o mapeamento...
Estou impressionado com o forecastpacote R , bem como, por exemplo, o zoopacote para séries temporais irregulares e interpolação de valores ausentes. Meu aplicativo é na área de previsão de tráfego de call center; portanto, os dados nos finais de semana estão (quase) sempre ausentes, o que pode...
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 6 anos . Qual é a diferença entre as funções rnorme...
Estou trabalhando em conjunto de dados electricitydisponíveis no pacote R TSA. Meu objetivo é descobrir se um arimamodelo será apropriado para esses dados e, eventualmente, se ajustará a eles. Portanto, procedi da seguinte maneira: 1º: Plote a série temporal que resultou se o seguinte gráfico: 2º:...
Acabei de ler um artigo em que os autores realizaram uma regressão múltipla com dois preditores. O valor geral ao quadrado de r foi 0,65. Eles forneceram uma tabela que dividiu o quadrado do r entre os dois preditores. A tabela ficou assim: rsquared beta df pvalue whole model 0.65 NA 2, 9...
Não estou acostumado a usar variáveis no formato de data em R. Estou apenas imaginando se é possível adicionar uma variável de data como uma variável explicativa em um modelo de regressão linear. Se for possível, como podemos interpretar o coeficiente? É o efeito de um dia na variável...
Estou tendo problemas para entender a saída do meu lmer()modelo. É um modelo simples de uma variável de resultado (Suporte) com interceptações de estado variáveis / efeitos aleatórios de estado: mlm1 <- lmer(Support ~ (1 | State)) Os resultados de summary(mlm1)são: Linear mixed model fit by...
Eu tenho dados para uma rede de estações meteorológicas nos Estados Unidos. Isso me fornece um quadro de dados que contém data, latitude, longitude e algum valor medido. Suponha que os dados sejam coletados uma vez por dia e conduzidos pelo clima em escala regional (não, não entraremos nessa...
Eu conheço três métodos para fazer estimativa de parâmetros, ML, MAP e Bayes. E para a abordagem MAP e Bayes, precisamos escolher os anteriores para os parâmetros, certo? Digamos que eu tenho esse modelo , no qual são parâmetros, para fazer a estimativa usando MAP ou Bayes, li no livro que é...
Em R, eu uso a ldafunção da biblioteca MASSpara fazer a classificação. Pelo que entendi LDA, a entrada xxx receberá o rótulo yyy , que maximiza p(y|x)p(y|x)p(y|x) , certo? Mas quando ajustar o modelo, no qual x=(Lag1,Lag2)x=(Lag1,Lag2)x=(Lag1,Lag2)y=Direction,y=Direction,y=Direction, eu não...
Estou usando a glmerfunção do lme4pacote em R e usando o bobyqaotimizador (ou seja, o padrão no meu caso). Estou recebendo um aviso e estou curioso para saber o que isso significa. Warning message: In optwrap(optimizer, devfun, start, rho$lower, control = control, : convergence code 3 from...