Eu tenho os seguintes dados e gostaria de ajustar um modelo de crescimento exponencial negativo a ele: Days <- c( 1,5,12,16,22,27,36,43) Emissions <- c( 936.76, 1458.68, 1787.23, 1840.04, 1928.97, 1963.63, 1965.37, 1985.71) plot(Days, Emissions) fit <- nls(Emissions ~ a* (1-exp(-b*Days)),...
Estou usando o pacote nnet no R para tentar construir uma ANN para prever preços de imóveis para condomínios (projeto pessoal). Eu sou novo nisso e não tenho formação em matemática, então, por favor, fique comigo. Eu tenho variáveis de entrada que são binárias e contínuas. Por exemplo, algumas...
Eu estava lendo sobre o cálculo da estimativa imparcial do desvio padrão e a fonte que li declarada (...) exceto em algumas situações importantes, a tarefa tem pouca relevância para aplicações de estatística, uma vez que sua necessidade é evitada por procedimentos padrão, como o uso de testes...
Eu queria saber qual pacote estatístico de software vocês recomendam para executar a inferência bayesiana. Por exemplo, eu sei que você pode executar o openBUGS ou o winBUGS como autônomo ou também pode chamá-los de R. Mas o R também possui vários de seus próprios pacotes (MCMCPack, BACCO), que...
mgcv::gamPágina de ajuda da Reading : intervalos de confiança / credibilidade estão prontamente disponíveis para qualquer quantidade prevista usando um modelo ajustado No entanto, não consigo descobrir uma maneira de realmente conseguir uma. Eu pensei predict.gamque teria um type=confidencee...
Preciso estimar a função densidade com base em um conjunto de observações usando o estimador de densidade do kernel. Com base no mesmo conjunto de observações, também preciso estimar a primeira e a segunda derivada da densidade usando as derivadas do estimador de densidade do núcleo. A largura de...
Estou tentando plotar um QQ-plot com dois conjuntos de dados de cerca de 1,2 milhão de pontos, em R (usando qqplot e alimentando os dados no ggplot2). O cálculo é fácil, mas o gráfico resultante é muito lento para carregar, porque há muitos pontos. Eu tentei a aproximação linear para reduzir o...
Eu tenho um monte de dados de duas amostras (controle e tratamento), cada uma contendo vários milhares de valores que serão submetidos a testes de significância em R. Teoricamente, os valores devem ser contínuos, mas devido ao arredondamento feito pelo software de medição, eles não são ' te eles...
Eu recentemente implementei um filtro Kalman no exemplo simples de medir uma posição de partículas com velocidade e aceleração aleatórias. Descobri que o filtro Kalman funcionava bem, mas então me perguntei qual era a diferença entre isso e apenas fazer uma média móvel. Descobri que, se eu usasse...
Eu apenas comecei a aprender a usar Stan e rstan. A menos que eu sempre tenha ficado confuso sobre como o JAGS / BUGS funcionou, pensei que você sempre tivesse que definir uma distribuição anterior de algum tipo para cada parâmetro no modelo a ser extraído. Parece que você não precisa fazer isso em...
Eu tenho um vetor de números que eu enviei aqui (... / code / MyData.Rdata) usando dput. Gostaria de obter o bca ci, então eu escrevi este código: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean, R=1000) Mas...
Eu usei randomForest para classificar 6 comportamentos de animais (por exemplo, em pé, andando, nadando etc.) com base em 8 variáveis (diferentes posturas e movimentos corporais). O MDSplot no pacote randomForest me fornece essa saída e tenho problemas em interpretar o resultado. Eu fiz um PCA...
Estou fazendo uma validação cruzada usando o método deixar de fora. Eu tenho uma resposta binária e estou usando o pacote de inicialização para R e a função cv.glm . Meu problema é que não entendo completamente a parte "custo" dessa função. Pelo que entendi, é a função que decide se um valor...
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á 2 anos . Eu tenho um conjunto de dados bastante...
Atualmente, estou trabalhando em um modelo de regressão em que tenho apenas variáveis categóricas / fatoriais como variáveis independentes. Minha variável dependente é uma proporção transformada em logit. É bastante fácil apenas executar uma regressão normal em R, pois R sabe automaticamente...
Recebi dados x = c(21,34,6,47,10,49,23,32,12,16,29,49,28,8,57,9,31,10,21,26,31,52,21,8,18,5,18,26,27,26,32,2,59,58,19,14,16,9,23,28,34,70,69,54,39,9,21,54,26) y = c(47,76,33,78,62,78,33,64,83,67,61,85,46,53,55,71,59,41,82,56,39,89,31,43,29,55,
Os dados que tenho são um valor de declive de regressão de y ~ time, um erro padrão, um valor n e um valor ap, para uma espécie específica em duas áreas diferentes. Quero verificar se a inclinação de regressão para uma área é significativamente diferente da inclinação de regressão para a outra área...
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 fique no tópico de Validação cruzada. Fechado há 6 meses . Ao testar a suposição de regressão paralela...
Eu sou novo na linguagem R. Eu gostaria de saber como simular a partir de um modelo de regressão linear múltipla que atenda às quatro suposições da regressão. ok .. obrigado. Digamos que eu queira simular os dados com base neste conjunto de
A fórmula a seguir está correta se eu quiser medir o erro padrão da mediana no caso de uma amostra pequena com distribuição não normal (estou usando python)? sigma=np.std(data) n=len(data)