Mande um email para o autor
paulo...@leg.ufpr.br (mailto:paulo...@leg.ufpr.br)
Ele é da ufpr. Acredito que ele seja membro dessa lista .
att
On Oct 4 2018, at 6:00 pm, Cesar Rabak via R-br
wrote:
>
> Para esmagar esse bug, pode-se setar uma semente específica em ambos os SO e
> ver se aí há
Edson,
Segue uma possível solução. Deixo a conferência dos resultados com você.
Tive que fazer algumas correções na tabela devido a forma como os dados estão
armazenados. Evite utilizar pontuação nas palavras (ex: Discordância, prefira
Discordancia ) e acredito que o nome das colunas podem ser
Mande um CMR porque montar uma tabela baseada em sua descrição ou criar
soluções de cabeça é isso. Bom mínimo de um dput dos dados ou parte
representativa deles, copie e cole aqui
Em qua, 26 de set de 2018 12:52, Edson Lira via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:
> Bom dia a todos! Tenho
Thallita
Embora as figuras que enviou não mostrem o ajuste dos modelos que você está
estudando. A escolha do modelo é um pouco mais complexa e não deve ser baseada
em palpite.
Se você procurar na literatura sobre métodos de comparação de modelos, você
encontrará uma vasta literatura.
Na minha
Cara pelo que parece você possui um banco de dados muito grande sendo carregado
e pouca memória RAM e Swap para carregar todo esse dado
On Aug 1 2018, at 1:20 pm, Andre Oliveira via R-br
wrote:
>
> Ubuntu 64, estou apenas repetindo uma leitura que funcionava normalmente.
>
> obg
> Em
Acho que isso resolve
a <- c(1,2,3,4,5)
b <- c(2,8,6,4,10)
dat <- data.frame(a,b)
dat[-c(2,4),]
ou
out<-c(2,4) # sem aspas
dat[-out,]
On Jul 28 2018, at 11:08 am, Antonio Silva via R-br
wrote:
>
> Olá,
>
> Algo que deve ser simples mas no qual empaquei.
>
> Tenho um vetor com o nome de
Caros colegas
Desejo fazer um gráfico, semelhante ao retornado pela função augPred{nlme},
contendo os valores observados, os valores preditos pelo modelo fixo, e os
valores individuais preditos. como esse exemplo (Só que para modelo não linear)
https://rdrr.io/cran/nlme/man/plot.augPred.html
Usando apply fica assim
variavel=apply(DADOS,1,sum)
Sent from Mailspring
(https://link.getmailspring.com/link/1531333852.local-7c031b20-c273-v1.2.2-96fb3...@getmailspring.com/0?redirect=https%3A%2F%2Fgetmailspring.com%2F=ci1ickBsaXN0YXMuYzNzbC51ZnByLmJy),
the best free email app for work
On
Uma solução simples
BANCO<- data.frame(cartao=c(1,2,3,4,5),emprestimo=c(6,7,8,9,10))
variavel<-DADOS$cartao +DADOS$emprestimo
variavel
Sent from Mailspring
Não sou usuário do expedes. Geralmente faço a análise usando a função
anova() . Seus dados estão normais. Você que tranformá-los porque? É devido
a variância? se for antes de buscar transformação faça uma análise mais
profunda dos seus dados para tentar compreender a causa. Se ainda assim não
Seja mais específica em sua definição de transformar de preferência envie
um CMR (código mínimo reproduzível) com seu problema para que possa obter
ajuda.
Acredito que você se refira a transformações dos dados para obter algum
pressuposto estatístico. Existem vários que dependem das
Faça um dput(data_estacoes) copie e cole a saída aqui para que possamos lhe
ajudar.
att
Em 7 de maio de 2018 16:12, Marcele de Jesus Correa via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:
> Olá, boa tarde!
>
> Sou iniciante no R e preciso de muita ajuda. Tenho recorrido aos
> tutoriais, livros,
No livro Singer, Nobre e Rocha (2017) disponível em:
https://www.ime.usp.br/~jmsinger/MAE0610/Singer
os autores citam como estratégia para identificação da estrutura da matriz
de variãncia covariãncia entre outros a construção do seguinte gŕafico:
" No caso de estudos longitudinais sem dados
Coloque um código reproduzível.
De um dput() de seus dados ou parte dele ,copie a saída e cole no e-mail.
Fica mais fácil te ajudar assim
Em seg, 26 de mar de 2018 11:59, Edmar Caldas via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:
> Eu estou usando RStudio Version 1.0.136.
>
> eu só fiz isso
Danielle, se erro está na forma como escreveu a condição. Se seu objetivo é
realizar uma sequência de comandos para valores maiores q a variável
limiteinferior então o teste e
For(i in 1:lenght(Aux){
If(Aux[i]>=limiteinferior){
Comandos para executar caso o teste lógico seja true
}
}
Supondo q Aux
Pessoal
Estou lendo uma tese onde a autora utilizou o modelo misto não linear para
avaliar a curva de degradação in vivo, utilizando o modelo de orskov (não
linear)
Ela fez toda a análise utilizando o modelo misto, para fazer a análise de
influência ela montou um modelo OLS utilizando a função
Caros amigos
Como é realizado análise de influência em modelos mistos não lineares?
Há alguma implementação em R para realizar tal análise ?
Vocês teriam alguma material para me passar ou indicar para que eu possar
ler e me inteirar sobre o tema?
Obrigado
--
Tem o livro do Douglas Bates e José Pinheiro, criadores do pacote nlme.
Para esse livro é o melhor livro para atender a sua demanda. É melhor q o
do bates
https://goo.gl/images/iULt4R
Em 09/03/2018 17:52, "Manasses Nóbrega via R-br"
escreveu:
> Olá a todos e todas,
Use FactoMineR
Seus desenvolvedores estão dando um curso Mooca q iniciou esta semana. Caso
interesse é possível se inscrever gratuitamente. Veja lá
https://googleweblight.com/i?u=https://www.r-bloggers.com/exploratory-multivariate-data-analysis-with-r-enroll-now-in-the-mooc/amp/=pt-BR
Em
http://www.leg.ufpr.br/~walmes/cursoR/cursoR2.pdf
Em 15/02/2018 15:13, "ana paula coelho madeira via R-br" <
r-br@listas.c3sl.ufpr.br> escreveu:
> Prezados colegas,
>
> boa tarde!
>
>
> Estou com um experimento em DIC no esquema de parcelas subdividas no
> tempo. No pacote ExpDes, temos
Espero que ajude.
library(agricolae)
data(sweetpotato)
model<-aov(yield~virus,data=sweetpotato)
comparison<- LSD.test(model,"virus",alpha=0.01,group=TRUE)
dados <- comparison$means
library(tibble)
dados <- rownames_to_column(dados, "TRAT")
dados <- rownames_to_column(dados, "ID")
A função revalue do pacote plyr é ótima para isso
dados <-
expand.grid(Resposta=c("otimo","bom","regular","ruim","NaoUtilizou","pessimo"),Questao=1:10)
head(dados)
library(plyr)
dados$novaclass <- revalue(dados$Resposta, c("otimo"="positiva","bom" =
"positiva","regular"="","NaoUtilizou" =
O pacote parece estar em desenvolvimento no gitHub. Veja se baixando por lá
e instalando vc não consegue instalar
No github não fala nada dele ser exclusivo para linux.
https://github.com/gabrielrvsc/HDeconometrics
Em 24 de janeiro de 2018 19:54, Mauro Sznelwar via R-br <
Um simples google
https://pt.wikipedia.org/wiki/Logaritmo_natural
Em 24 de janeiro de 2018 19:34, Mauro Sznelwar via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:
> Logaritmos neperiano é na base 10? É na base 'e' que é aprox. 2,7.
>
> Boa tarde,
>
> desculpe a demora em responder.
> @Cesar
Yury,
segue solução,
Basicamente consite em aumentar espaços na margem com a função par() e
então modificar o rótulo do eixo com a função title(), você alter a
distancia modificando o argumento "line", e o tamanho da fonte em "cex.lab"
Att
a.box = matrix(sample(1:500, 100, replace = T), ncol =
Yuri,
Seu código não é reproduzível. Envie um código reproduzível para que possa
obter ajuda de forma eficiente.
abcs
Em 24 de janeiro de 2018 10:43, Yury Duarte via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:
> Bom dia colegas listeiros!
>
> Estou encontrando um pequeno problema e gostaria de
Caros amigos
Como faço para criar uma coluna ID no data frame abaixo que identifique as
observações de um indivíduos com um único numero.
Cada linha representa uma observação em um indíviduo que recebeu um
tratamento. Cada indivíduo teve 6 observações (SEMANA).
Desejo identificar as observações
Caros amigos, no exemplo abaixo, como faço para colocar as identificações
do indivíduo em ordem crescente?
att
library(nlme)
library(lattice)
install.packages("tidyr")
library(tidyr)
DADOS.CURTO <- read.csv("https://www.dropbox.com/s/b4cckybgrcg86me/
galinhas.csv?raw=1
Caros colegas,
Nos comandos abaixo, realizo a análise de um modelo não linear misto. Após
selecionar o modelo, realizei um teste outliers do modelo. Identifiquei os
outliers e os removi. Ajustei novamente o modelo, porém a nova análise
gráfica dos resíduos desses modelo indica novos outliers.
Caros Amigos
o que pode significar essa mensagem em modelos não linear mixto (nlme)?
Error in MEEM(object, conLin, control$niterEM) :
Singularity in backsolve at level 0, block 1
--
=
Fernando Souza
Zootecnista, DSc. Produção e Alimentação Animal
Muito obrigado pelo retorno Walmes!
abcs
Em 12 de dezembro de 2017 10:56, Walmes Zeviani via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:
> A nlme usa a lattice. Então você muda da mesma forma que mudaria um objeto
> que veio da xyplot, por exemplo.
>
> GALINHA$ID <- factor(GALINHA$ID,
>
Caros colegas,
Devido ao elevado número de individuos avaliados há uma sobreposição dos
identificados no eixo y do gráfico, o que o torna incompreensível.
Como posso organizar esses indicadores? Ou como posso removê-los.
library(nlme)
library(lattice)
GALINHA<-structure(list(ID =
t;paisagem" e faça o ajuste da figura (*cropp*) dela.
> . .
>
> HTH
> --
> Cesar Rabak
>
>
> On Mon, Dec 4, 2017 at 8:34 PM, Fernando Souza via R-br <
> r-br@listas.c3sl.ufpr.br> wrote:
>
>>
>> Caros colegas
>>
>> Criei esses gráfico
install.packages("tidyr")
library(tidyr)
dados<-structure(list(id = 1:6, UF = structure(1:6, .Label = c("Acre",
"Alagoas", "Amapá", "Amazonas", "Bahia", "Ceará"), class = "factor"),
x2002 = c(192L, 1251L, 204L, 958L, 6885L, 3705L), x2003 = c(189L,
1269L, 198L, 1050L, 6901L, 3871L), x2004
Talvez possa ajudar
https://stackoverflow.com/questions/23547929/determine-if-a-vector-is-ordered
Em 25/10/2017 18:31, "Elias Carvalho via R-br"
escreveu:
> Como posso identificar se uma variável é 'ordered factor nas seguintes
> condições:
>
>- o read.csv leu ela
Eu faria assim
datanasc<-as.Date("02/06/1982","%d/%m/%Y")
idade<-difftime(Sys.Date(),datanasc) /365.25
Em 24 de outubro de 2017 18:57, Jônatan via R-br
escreveu:
> 'dttm' é a abreviatura para 'datetime' usada no print de data_frame ou
> tibble (tbl) do pacote
Uma solução.
O comando retorna um dataframe contendo uma coluna com as médias e outra
com os desvios. A partir daí fazer os graficos fica a seu sabor.
populacao <- rchisq(10, df = 10)
tamanhoAmostra <- c(6,14,16,18,20)
##x:dados a serem amostrados
##n= tamanho da amostra
##r =número de
Acredito que isso deva te ajudar. Não entendi muito bem o que quis dizer
com " para cada tam repetir nam (vezes)" , por isso o código abaixo não
faz nenhuma repetição dos valores amostrados , porém acredito que é fácil
de ajustar.
populacao <- rchisq(10, df = 10)
tamanho <-
Amigos como faço para colocar a linha dos valores preditos e intervalo de
confiança da predição sobre esses gráfico?
Fiz uma tentativa, o gráfico não ficou agradável
GALINHA<-structure(list(ID = structure(c(128L, 48L, 109L, 52L, 42L, 16L,
210L, 37L, 281L, 96L, 255L, 1L, 189L, 323L, 77L, 184L,
Caros colegas ao rodar a função augPred {nlme} recebo a seguinte mensagem:
"Error in tapply(as.character(object[[nm]]), groups, FUN[[dClass]]) :
arguments must have same length"
Meus dados possuem diferentes números de individuos por tratamento e há
alguns NA's. Acredito que é por isso que nao
Os comando so Fernando Toledos estão corretos, porém necessitava de algumas
correções na função melt
Os pacotes necessários para o código são: {plyr} e {reshape}
setwd("~/Documentos/Biblioteca/EstudoR" )
#raiz =
'C:\\Users\\Yury\\Desktop\\Mestrado\\1_TESE\\Model_Data\\Maize\\PREVISAO\\PA\\'
Caros colegas
Possuo R instalado em uma máquina Debian 8. Recentemente ao utilizar o
pacote nome recebi algumas mesgs estranhas de erro. Fui pesquisar e parece
estar relacionado a falhas de compilação devido ao período de congelamento
de debian8 para passaram para debian9. Já testei e o mesmo
install.packages(car)
library(car)
Anova(modelo, type="lll")
Em 03/07/2017 23:20, "Fabiana Goncalves Barbosa via R-br" <
r-br@listas.c3sl.ufpr.br> escreveu:
> Olá Pessoal
>
> Como faço uma ANOVA para amostras de tamanhos diferentes?
>
>
>
> ___
> R-br
Se você possui 4 GL para o fator A e 4 GL para o fator B, como podem haver
somente 1 GL para a interação? Acredito que por isso o R retorna "there are
aliased coefficients in the model". Uma vez que você não possui informações
em todas as combinações dos fatores A e B.
Vc pode armazenar em um array ou trabalhar com listas. Veja documentação
Em 19/06/2017 13:12, "Clodoaldo José Figueredo via R-br" <
r-br@listas.c3sl.ufpr.br> escreveu:
> Caros(as)
>
> Estou replicando dados e calculando a matriz de covariância em cada
> reamostragem.
> Preciso armazenar essas
Caros colegas,
Estou realizando um exercício de parcela subdividida utilizando a
função lme{nlme}. Acontece que quando rodo a função summary.glht as
comparações pairwise são feitas e eu recebo as seguintes mensagens
warnings:
Warning messages:
1: In RET$pfunction("adjusted", ...) :
Em Sat, 21 Jan 2017 10:15:10 -0300
Marcos Vital via R-br escreveu:
> Olá, colegas de lista!
>
> Peço desculpas para aqueles que estão lá na comunidade de R no
> facebook, para quem isso será um pouco de cross-posting...
>
> Estou construindo uma lista de blogs
Eu trabalho assim
#-Analise de Cook's distance
---
install.packages(sfsmisc); library(sfsmisc)
analise<-lm(CONSUMO~factor(GEST)*factor(MANEJO),data=agua)
n<-length(agua$CONSUMO) # número de observações
48 matches
Mail list logo