Já faz um tempo desde que eu olhei para particionar árvores. A última vez que fiz esse tipo de coisa, eu gosto de festa no R (criada por Hothorn). A idéia de inferência condicional via amostragem faz sentido para mim. Mas a rpart também teve apelo. No aplicativo atual (não posso fornecer detalhes,...
Houve várias boas perguntas e conjuntos de respostas em livros introdutórios ou abordagens para aprender R, por exemplo, aqui e aqui . Mas eu tenho um problema um pouco diferente - a melhor maneira de executar uma sessão de uma hora (ou várias dessas sessões) em um laboratório de informática que...
Estou analisando dados de um experimento fatorial desequilibrado com SASe R. Ambos SASe Rfornecer semelhantes Tipo I soma dos quadrados, mas a sua soma Tipo III de quadrados são diferentes uns dos outros. Abaixo estão SASe Rcódigos e saídas. DATA ASD; INPUT Y T B; DATALINES; 20 1 1 25 1 2 26 1...
O CARET usará automaticamente uma grade de ajuste pré-especificada para criar vários modelos antes de selecionar um modelo final e depois treinar o modelo final com os dados completos do treinamento. Eu posso fornecer à minha própria grade de sintonia apenas uma combinação de parâmetros. No...
Como calcular a estatística do quadrado R ( r2r2r^2 ) em R para loesse / ou predictsaída da função? Por exemplo para esses dados: cars.lo <- loess(dist ~ speed, cars) cars.lp <- predict(cars.lo, data.frame(speed = seq(5, 30, 1)), se = TRUE) cars.lppossui duas matrizes fitpara modelo e...
Bloqueado . Esta pergunta e suas respostas estão bloqueadas porque a questão está fora do tópico, mas tem um significado histórico. No momento, não está aceitando novas respostas ou interações. Estou tendo o seguinte problema ao fazer algumas análises com R. Eu tenho um...
Thru este site que eu descobri recentemente Sankey diagramas, uma ótima maneira de visualizar o que está acontecendo em um tradicional fluxograma. Aqui está um bom exemplo de um diagrama de Sankey de George M. Whitesides e George W. Crabtree , Source; Não se esqueça da pesquisa fundamental de...
Eu tenho uma pergunta sobre o uso de uma variável de agrupamento em um modelo não linear. Como a função nls () não permite variáveis de fator, tenho me esforçado para descobrir se é possível testar o efeito de um fator no ajuste do modelo. Incluí um exemplo abaixo, no qual desejo ajustar um...
Estou usando rlm no pacote R MASS para regredir um modelo linear multivariado. Funciona bem para várias amostras, mas estou obtendo coeficientes quase nulos para um modelo específico: Call: rlm(formula = Y ~ X1 + X2 + X3 + X4, data = mymodel, maxit = 50, na.action = na.omit) Residuals: Min 1Q...
Estou tentando encontrar o máximo local para uma função de densidade de probabilidade (encontrada usando o densitymétodo de R ). Não consigo fazer um método simples de "olhar ao redor dos vizinhos" (onde se olha em torno de um ponto para ver se é um máximo local em relação aos vizinhos), pois há um...
Suponha que tenhamos X1∼unif(n,0,1),X1∼unif(n,0,1),X_1 \sim \textrm{unif}(n,0,1), X2∼unif(n,0,1),X2∼unif(n,0,1),X_2 \sim \textrm{unif}(n,0,1), onde unif(n,0,1)unif(n,0,1)\textrm{unif}(n,0,1) é uma amostra aleatória uniforme de tamanho n,
Estou fazendo uma regressão de riscos proporcionais de Cox em R usando coxph, que inclui muitas variáveis. Os resíduos de Martingale parecem ótimos, e os resíduos de Schoenfeld são ótimos para QUASE todas as variáveis. Existem três variáveis cujos resíduos de Schoenfeld não são planos e a...
Vi métricas de erro diferentes usadas nas competições do Kaggle: RMS, média quadrada, AUC, entre outras. Qual é a regra geral da escolha de uma métrica de erro, ou seja, como você sabe qual métrica de erro usar para um determinado problema? Existem
Estou fazendo um curso de pós-graduação em Estatística Aplicada que usa o seguinte livro (para dar uma idéia do nível do material coberto): Conceitos e Métodos Estatísticos , de GK Bhattacharyya e RA Johnson. O professor exige que usemos o SAS para os trabalhos de casa. Minha pergunta é: existe...
Um dos problemas que sempre tive com modelos mistos é descobrir visualizações de dados - do tipo que pode acabar em um papel ou pôster - depois de obter os resultados. No momento, estou trabalhando em um modelo de efeitos mistos de Poisson com uma fórmula semelhante à seguinte: a <-...
MAD = Desvio Absoluto Médio MSE = Erro Quadrático Médio Vi sugestões de vários lugares em que o MSE é usado, apesar de algumas qualidades indesejáveis (por exemplo , http://www.stat.nus.edu.sg/~staxyc/T12.pdf , que afirma na p8 "Acredita-se que MAD é um critério melhor que o MSE. No entanto,...
Estou usando o pacote gbm para classificação. Como esperado, os resultados são bons. Mas estou tentando entender a saída do classificador. Existem cinco termos na saída. `Iter TrainDeviance ValidDeviance StepSize Improve` Alguém poderia explicar o significado de cada termo, especialmente o...
A abordagem comum para estimar os parâmetros de uma distribuição normal é usar a média e o desvio / variância padrão da amostra. No entanto, se houver alguns discrepantes, a mediana e o desvio médio da mediana devem ser muito mais robustos, certo? Em alguns conjuntos de dados que eu tentei, a...
Eu tenho um banco de dados de eventos (ou seja, uma variável de datas) e covariáveis associadas. Os eventos são gerados pelo processo Poisson não estacionário, com o parâmetro sendo uma função desconhecida (mas possivelmente linear) de algumas covariáveis. Eu acho que o pacote NHPoisson existe...
O pacote de plotagem R ggplot2 possui uma função impressionante chamada stat_smooth para plotar uma linha de regressão (ou curva) com a banda de confiança associada. No entanto, estou tendo dificuldade em descobrir exatamente como essa faixa de confiança é gerada, para cada momento da linha de...