Skip to main content

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 ...
Marcus Nunes's user avatar
  • 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 ...
Carlos Eduardo Lagosta's user avatar
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 ...
Daniel Falbel's user avatar
  • 12,5mil
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 ...
Marcus Nunes's user avatar
  • 19mil
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 ...
Rui Barradas's user avatar
  • 19mil
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 ...
Rui Barradas's user avatar
  • 19mil
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 ...
Rui Barradas's user avatar
  • 19mil
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 ...
Carlos Eduardo Lagosta's user avatar
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) # ...
Robert's user avatar
  • 940
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 ...
Rui Barradas's user avatar
  • 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 ...
Rui Barradas's user avatar
  • 19mil
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 ...
Ruben O. Chiavone's user avatar
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 ...
brow-joe's user avatar
  • 2.214
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 ...
Rui Barradas's user avatar
  • 19mil
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....
Marcus Nunes's user avatar
  • 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 ...
Marcus Nunes's user avatar
  • 19mil
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 ...
Jorge Mendes's user avatar
  • 1.634
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 ...
Rui Barradas's user avatar
  • 19mil
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 ...
Carlos Eduardo Lagosta's user avatar
3 votos
Aceito

Como aplicar uma função a vários objetos de uma vez?

A forma como eu gosto de fazer isso é criando uma lista de modelos no próprio objeto para fácil controle e organização com o comando nest do pacote tidyr. Primeiro eu crio uma função que descreva o ...
M Aurélio's user avatar
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 ...
Daniel Falbel's user avatar
  • 12,5mil
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....
Marcus Nunes's user avatar
  • 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 ...
Carlos Cinelli's user avatar
  • 17,2mil
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....
AlexCiuffa's user avatar
  • 2.459
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 ...
AlexCiuffa's user avatar
  • 2.459
2 votos
Aceito

ValueError: Input contains NaN, infinity or a value too large for dtype('float64') ( fiz a verificação e nenhuma dessas procede)

Sem acesso aos dados não consigo testar para ver se é isso, mas note que as linhas data.drop(['budget'],axis = 1) data.dropna() não são operações in-place. Tanto DataFrame.drop quanto DataFrame....
jfaccioni's user avatar
  • 1.931
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 ...
Rui Barradas's user avatar
  • 19mil
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 =...
bbiasi's user avatar
  • 814
2 votos
Aceito

Reta de regressão ajustada considerando diferentes fatores no ggplot2

Em primeiro lugar, para ler os dados é melhor, mais simples, usar a função help("read.csv"), que é a versão de read.table para ficheiros CSV com sep = ";" e "dec = ",&...
Rui Barradas's user avatar
  • 19mil
2 votos

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

Como complemento à excelente resposta do Jorge Mendes, esta resposta separa os resíduos por níveis de DummyVariable no gráfico p1b. Isso é feito com o aes(group = DummyVariable). O segundo gráfico, p2,...
Rui Barradas's user avatar
  • 19mil

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