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 ...
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
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 ...
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 ...
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 ...
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 ...
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 ...
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)
# ...
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 ...
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 ...
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 ...
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 ...
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....
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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....
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
...
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
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
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....
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 ...
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 =...
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 = ",&...
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,...
Only top scored, non community-wiki answers of a minimum length are eligible
Tags relacionadas
regressão × 74r × 55
plot × 13
ggplot2 × 11
python × 9
modelagem × 8
matlab × 4
lm × 4
python-3.x × 3
gráfico × 3
machine-learning × 3
sklearn × 3
glm × 3
banco-de-dados × 2
data.frame × 2
model × 2
lme4 × 2
postgresql × 1
funções × 1
ios × 1
loop × 1
variáveis × 1
swift × 1
for × 1
dataframe × 1