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 ...
- 7.746
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 ...
- 7.746
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 ...
- 7.746
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
Tags relacionadas
regressão × 74r × 54
plot × 13
ggplot2 × 11
python × 9
modelagem × 8
matlab × 4
lm × 4
banco-de-dados × 3
python-3.x × 3
rstudio × 3
gráfico × 3
machine-learning × 3
sklearn × 3
glm × 3
data.frame × 2
model × 2
lme4 × 2
postgresql × 1
funções × 1
excel × 1
ios × 1
loop × 1
variáveis × 1
swift × 1