Respostas interessantes marcadas com a tag

8 votos
Aceito

Valores iniciais na função nls no R

A função nls é utilizada para fazer uma regressão não-linear. Ela utiliza processos iterativos para obter as estimativas dos parâmetros do teu modelo. Em teoria estes processos iterativos vão pouco a ...
  • 19mil
6 votos
Aceito

Data Frame e Regressão Linear

Para fazer o log você pode colocar diretamente na fórmula. Exemplo: # sem log lm(mpg ~ cyl, mtcars) # com log lm(log(mpg) ~ cyl, mtcars) Para colocar o termo quadrático você vai utilizar a função ...
  • 17,1mil
5 votos

Como colocar a linha de regressão em um gráfico?

Uma outra maneira é usar a função curve. curve(coef(modelo)[1]+coef(modelo)[2]*x,add=TRUE,col = "blue",lwd=2) Uma vantagem dessa função é que ela também serve caso o modelo envolva componentes não ...
  • 2.373
5 votos
Aceito

Problema com forecast no R

O erro está dizendo o problema: as variáveis t e t2 possuem comprimentos diferentes da sua variável resposta. Veja o erro que aparece nesse código bem simples em que as variáveis possuem ...
  • 12,4mil
5 votos
Aceito

Variaveis selecionadas no GLM sendo utilizadas na GLMM

nos modelos que tenho testado em meu trabalho, o que apresenta o melhor resultado (AIC BIC etc) não é o GLMM que partiu de todas as variáveis que possuo e sim o GLMM que parte das variáveis mais ...
  • 19mil
5 votos
Aceito

Gráfico de Regressão no R ou Python

As linhas no gráfico que quer reproduzir não correspondem ao resultado de uma regressão linear, mas simplesmente à porcentagens do valor de x. Pode adicioná-las com stat_function. Se for usar em ...
5 votos
Aceito

gráfico de resíduos half-normal plots no ggplot2

Segundo a secção Value de help("hnp), o objeto Graph16 é uma lista de classe "hnp". Estes objetos têm membros x, lower, upper, median e residuals, que são os vetores representados no ...
  • 18,4mil
5 votos
Aceito

Modelo heterocesdástico de efeitos mistos via função lmer

Numa resposta do R-sig-mixed-models, Douglas Bates (GoogleScholar, ResearchGate) dá a seguinte resposta a uma pergunta sobre a possibilidade de modelar o problema da heteroscedasticidade da ...
  • 18,4mil
5 votos

Aggregate dados no R com função própria

Eis uma solução. Mas com a função by, não com aggregate. A função by divide a tabela por Categoria e aplica o lm a cada sub-base x. Depois o ciclo lapply cria uma lista de sumários, onde se incluem os ...
  • 18,4mil
5 votos
Aceito

Aggregate dados no R com função própria

Se quer apenas os valores de R² por categoria, sem armazenar as listas com o resultado das regressões, pode se valer da lógica "separa/aplica" (split/apply). Pacotes dplyr e data.table ...
4 votos
Aceito

Como proceder regressão não linear pelo modelo mitscherlich?

O ajuste não é bom mas segue: fm <- nls(COM ~ Alfa*(1-10^(-Gama*(TEMPO+Tetha))), data = dados, start = list(Alfa = 900, Gama = .1, Tetha = 1)) summary(fm) # ...
  • 897
4 votos
Aceito

Quantidade de pontos fora do half-normal plot em ggplot2

Eis uma solução. São criadas duas funções, hnp_texto e hnp_cores. Só a primeira destas funções é pedida na pergunta. A função hnp_cores é opcional. É ainda criado uma função theme_* para simplificar ...
  • 18,4mil
3 votos
Aceito

Regressão linear com python

Pelo o que li no site do scipy roda no Windows sim. [Google Translate] Para a maioria dos usuários, especialmente no Windows, a maneira mais fácil de instalar os pacotes da pilha SciPy é baixar uma ...
3 votos

Regressão linear com python

Segundo a própria documentação do scipy ele não roda muito bem em Windows pois tem algumas dependências que funcionam apenas em linux e mac Como alternativa, recomendo utilizar o sklearn, ela é uma ...
  • 2.154
3 votos
Aceito

R regressão em dados da mesma linha

Não é assim tão difícil como isso. Na solução abaixo Calculo o logaritmo de cada coluna exceto a primeira logo no início. Aplico (apply) o modelo lm(x ~ y) a cada linha. E depois há várias instruções ...
  • 18,4mil
3 votos

Regressão linear ponderada utilizando o inverso da variância como fator de ponderação

Basta criar o vetor de pesos desejado para resolver este problema. No teu caso, chamei este vetor de pesos: variancias_condicionais <- aggregate(df$Y, list(df$X), var)$x quantidade_X <- as....
  • 19mil
3 votos

Regressões lineares em DIC com parcela subdividida

A primeira coisa a se fazer quando vamos analisar um experimento é a análise exploratória dos dados. Não ficou explícito na tua pergunta, mas estou assumindo que o Fator1 diz respeito às parcelas e o ...
  • 19mil
3 votos
Aceito

Por que a função "nls" do R não apresenta convergência mesmo com a curva inicial bem próxima dos pontos?

O problema está no excesso de parâmetros para tão poucos pontos. Se a função a ajustar está a usar uma translação do eixo dos x representada por x - C, pode-se eliminar um dos parâmetros, C, ajustando ...
  • 18,4mil
3 votos
Aceito

Gráficos dos resíduos do modelo de efeitos mistos utilizando a função ggplot2

Olha, pelos nomes (.fitted, .resid) os resultados parecem ser ligados ao pacote broom, que usa esse padrão pros nomes das colunas. (ou mais especificamente o broom.mixed para modelos lme) Com a ...
  • 1.633
3 votos
Aceito

Curvas de densidade estimada por um modelo em gamlss

A solução é traçar os gráficos das densidades dos valores previstos e não esses valores. Em relação à distribuição Gama, o gamlss, função GA, ajusta um modelo logarítmico para os parâmetros mu e sigma ...
  • 18,4mil
3 votos
Aceito

Densidades de probabilidade estimadas por um modelo de regressão com covariáveis

Pode calcular média e desvio dos valores preditos e usá-los para estimar os parâmetros das funções de densidade de probabilidade, plotando-as para a extensão dos dados. Não uso o pacote gamlss. Usarei ...
2 votos

Como colocar a linha de regressão em um gráfico?

É possível também trabalhar com gráfico interativo, caso necessário. Para tal, fica a dica do pacote highcharter. set.seed(1) x <- rnorm(100) y <- rnorm(100) + 2*x +10 dados <- data.frame(x =...
  • 814
2 votos

Como fazer Multiple Linear Regression?

Como você não forneceu dados de exemplo do seu domínio de problema (apesar de eu ter pedido duas vezes), eu tive que usar alguma fonte de dados publicamente disponível na Internet. Optei por usar os ...
  • 34,4mil
2 votos
Aceito

O que significa o argumento alpha da função ltsReg, do pacote robustbase?

As estimativas deverão ser equivalentes quando você usar alpha = 1. Veja o exemplo abaixo: > library(robustbase) > data(heart) > coef(ltsReg(clength ~ height + weight, data = heart, alpha = ...
  • 12,4mil
2 votos
Aceito

resultados diferentes usando rpart e caret

O caret por padrão faz tuning de alguns hiperparâmetros de cada modelo. Ele tenta fazer isso de uma forma inteligente, mas que nem sempre é o adequado para o seu problema. Já o rpart ajusta o modelo ...
  • 12,4mil
2 votos
Aceito

Regressão Linear Múltipla no R em grupos diferentes

Uma maneira de resolver este problema é através dos pacotes dplyr e broom: library(dplyr) library(broom) iris.regressao <- iris %>% group_by(Species) %>% do(regressao = lm(Sepal....
  • 19mil
2 votos

Modificar teste t para parâmetros de regressão linear no R

Com o R base você pode pedir os intervalos de confiança com a função confint: confint(reg) 2.5 % 97.5 % (Intercept) -0.005635243 0.7256352 x 0.849756826 1.0702432 ...
  • 17,1mil
2 votos

Curva de regressão com eixo x igual a 0

O regressor LinearRegression() do sklearn tem o atributo intercept_, que retorna o y onde o regressor intercepta o eixo Y, ou seja, em x = 0. Inclusive, no exemplo da documentação do sklearn....
  • 2.424
2 votos
Aceito

Como utilizar um modelo de regressão quadrática?

Utilizando apenas o statsmodels: Com o statsmodels é possível escrever a fórmula desejada, como por exemplo: target ~ np.power(X1, 2) + X2 Neste exemplo, significa que estamos buscando os ...
  • 2.424
2 votos

Como ajustar a reta de regressão de modo que 90% dos fiquem abaixo da reta?

Aqui estão duas maneiras de resolver o problema utilizando o pacote quantreg. A fórmula y = a*x^b pode ser transformada aplicando logaritmos e ajustar o modelo resultante, ou seja, uma reta de ...
  • 18,4mil

Only top scored, non community-wiki answers of a minimum length are eligible