Você nunca conseguirá *testar* a normalidade em uma correlação. . .
E o uso do método de Spearman *não* é para dados não normais (*apenas*) mas
sim para dados cuja medida não seja intervalar e a única forma de
interpretar as diferenças seja via postos (*ranks* na língua de
Shakespeare).
O que
Irmão, como quer avaliar correlação, acho que o teste de normalidade deve ser
por variável...
Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia
- PUC-Rio/PósMQICONRE: 8514 - SÉRIE AEmail:
diogojos...@yahoo.com.brhttp://lattes.cnpq.br/8996149312896520
Em
Oi Diogo, boa tarde. Entendi sim, mas ai tem que ser uma por uma, ou tem
alguma maneira de testar com as variaveis juntas?
Abraços,
*Rodrigo Zucaratto*
Bacharel em Ciências Biológicas
Mestre em Ciências - Conservação da Natureza
Doutor em Ecologia e Evolução
Laboratório de Ecologia de Mamíferos
Boa tarde a todos!!!
Rodrigo, você testa normalidade nos dados que vão gerar a correlação, não na
correlação em si.
Você deve testar se suas variáveis possuem distribuição normal (aí tem vários,
Kolmogorov-Smirnov com a correção de Lilliefors, Anderson-Darling,
Shapiro-Wilks...). Se uma dessas
Prezados,
Estou com uma duvida em relação ao teste de correlação. Quando faço uma
correlação, não ha uma relação de causa e efeito (até aqui entendo). No
entanto, estou querendo avaliar a inluência de variáveis
ambientais (temperatura e pluviosidade) na fenologia de uma palmeira. To
fazendo a
Prezado Alexandre,
De fato, frente à idade do texto, é prudente observar se algo de mais novo
pode aplicar-se à sua área de conhecimento.
SE sua obra tiver apenas quatro *escalas* (nominal, ordinal, intervalar, e
razão [racional]), pode ser que os questionamentos modernos lhe venham a
ajudar
Alexandre,
A teoria sobre tipos ou categorias de dados não tem uma fundamentação
teórica muito sólida, e é frequentemente desafiada por textos fora do
*mainstream* dos livros texto de estatística, que noves fora uma ou outra
minúcia ainda seguem muito do proposto por Stevens há mais de 50 anos.
Prezados
Tenho formatado pesquisas na empresa. Elaboro questionários e formulários,
aprendi a usar R, R-Studio e tenho alguma noção de banco de dados e SPSS.
Tenho alguma familiaridade no emprego de testes estatísticos.
Entretanto, persiste uma dúvida sobre os tipos de testes que podem ser
Perdão por não ter sido tão claro, mas já consegui resolver o problema.
Obrigado pessoal :)
Em 4 de julho de 2017 18:13, Cesar Rabak via R-br
escreveu:
> De fato, Fernando,
>
> Está difícil de entender o quê você deseja. . .
>
> Você fala que deseja saber "a
De fato, Fernando,
Está difícil de entender o quê você deseja. . .
Você fala que deseja saber "a proporção de acerto considerando a amostra B"
mas usa o valor da amostra A. . .
A razão da subtração na sua segunda questão é ainda mais obscura, pois
algebricamente ela não pode dar o resultado
Não entendi o que você quer fazer.
Em 3 de jul de 2017 12:03 PM, "Fernando Gama via R-br" <
r-br@listas.c3sl.ufpr.br> escreveu:
> Pessoal estou com uma duvida, estou com dificuldades de recordar, tenho
> duas amostras A e B, para ilustrar:
>
> A <- length(216) # o tamanho é variável para cada
Pessoal estou com uma duvida, estou com dificuldades de recordar, tenho
duas amostras A e B, para ilustrar:
A <- length(216) # o tamanho é variável para cada laço for
B <- length(87)
Suponhamos que eu tenha um sample
sample <- length(15)
Gostaria de obter a proporção de acerto considerando o
Veja a função HSD.test do pacote agricolae.
HTH
--
Cesar Rabak
2017-04-08 14:44 GMT-03:00 Tiago G. Santos via R-br <
r-br@listas.c3sl.ufpr.br>:
> Boa tarde,
> Alguém sabe como faço para apresentar os resultados dos testes tukey ou
> dunnet como as letras que indicam significância estatística
Seu desbalanço é por diferentes números de repetições? se for o caso você
pode analisar normalmente utilizando a função aov e realize a anova tipo
III
install.packages("car")
library(car)
library(lsmeans)
Anova(modelo.aov, type="III")
Desdobramentos deverá ser feito, caso haja efeito
Boa tarde Eu estou tentando rodar os dados de um experimento com dois fatores,
os dados estão desbalanceadosEu preciso fazer o desdobramentos desses
dados?Outra duvida é sobre os teste de comparação de médias para DIC com dados
desbalanceados, no teste de Tukey não houve diferenças para p<0.05,
Leonardo,
Eu já toquei nesse assunto em resposta a outra questão parecida com essa
sua, por isso serei mais breve.
A linguagem R é por excelência uma linguagem voltada para o tratamento de
dados em grupos que *grosso modo* podemos denominar de matriciais.
Por a maioria das coisas que se pode
Leonardo, não sei se entendi o que quer.
Tente ou
ifelse(x 4, a, b)
ou
ifelse(x4,print(a),print(b))
Em 23 de agosto de 2015 23:27, Leonardo Monteiro
monteiroleonar...@gmail.com escreveu:
Prezados, eu sei que existem outras maneiras (funcoes, tipo ifelse) de se
fazer isso, porem eu gostaria
Prezados, eu sei que existem outras maneiras (funcoes, tipo ifelse) de se
fazer isso, porem eu gostaria de fazer pelo if else , pq fica mais claro no
meu script.
bom.. eu venho tentando aplicar a condicao em um vetor, por exemplo
x- 1:10
if(x 4) {print (a)} else {print (b)}
porem isso nao
Leonardo,
você pode passar bold() dentro de uma expression:
plot(x,y, main = expression(bold(paste(dados
D
On Fri, Jul 17, 2015 at 9:51 PM, Paulo Abreu abre...@me.com wrote:
Não serve usar o comando title()?
Assim:
x- 1:10
y = 1/x
plot (x,y)
title(dados, font=2)
Dê uma olhada no
Pessoal, tenho o seguinte:
x- 1:10
y = 1/x
plot (x,y, main = expression (dados))
preciso que o título (no caso, dados) seja em negrito.
alguma sugestão?
Obrigado,
Leonardo.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
Não serve usar o comando title()?
Assim:
x- 1:10
y = 1/x
plot (x,y)
title(dados, font=2)
Dê uma olhada no help de par
Em 17/07/2015, à(s) 21:42, Leonardo Monteiro monteiroleonar...@gmail.com
escreveu:
Pessoal, tenho o seguinte:
x- 1:10
y = 1/x
plot (x,y, main = expression
Estou aprendendo a ajustar modelos com distribuição não normal. Tenho lido
alguns livros para entender a metodologia, entretanto alguns pontos não
estão claro para mim. Acredito que porque minha mente está viciada em
modelos normais apenas. Por isso gostaria de uma esclarecimento e se
tiverem
Boa tarde, estou com uma dúvida de como usar a função paste ou se tem como
eu fazer de outro modo.
Eu tentei primeiramente assim:
matrizAA-matrix(nrow=395, ncol=31)
Symbols-c(AA,AIG,AXP)
for(j in 2:3)
{
empresa-Symbols[j]
for (i in 2:4){
matrizAA[i-1,j]-cor(AA[,i],Symbols[j][,i],
Um começo pode ser a leitura desses dois trabalhos:
http://www.leg.ufpr.br/lib/exe/fetch.php/projetos:artigo_silvia-rs.pdf
Ambos são trabalhos que mapeiam o risco de um evento de interesse
(processo espacial pontual) através de modelos aditivos generalizados. O
primeiro usa a função Kernel
Ola.
Gostaria de saber se alguém poderia me sugerir algo para avaliar o risco,
ou aonde uma doença ocorre espacialmente, mais comumente.
Meus dados são constituídos de coordenadas geográficas de cidades, na forma
de latitude e longitude. Dai para cada cidade eu tenho o número de casos de
uma
Veja o pacote spatstat para mapear padrões. Supondo q a doenca pode
acontecer em toda area de estudo.
Em 25/02/2015 18:21, Augusto Ribas ribas@gmail.com escreveu:
Ola.
Gostaria de saber se alguém poderia me sugerir algo para avaliar o risco,
ou aonde uma doença ocorre espacialmente, mais
Prezados, boa tarde.
Sendo as matrix x e y:
x- matrix(rnorm(25, 10, 4),nrow = 5, ncol = 5); x #matrix 1
y- matrix(rnorm(45, 50, 10), nrow = 9, ncol = 5); y # matrix 2
alguem poderia me ajudar em como pegar as linhas que estão sobrando da
matrix y (6 até 9) e colar na matrix x (abaixo da linha
rbind(x,y[6:9,])
2014-08-20 9:41 GMT-03:00 Leonardo Monteiro monteiroleonar...@gmail.com:
Prezados, boa tarde.
Sendo as matrix x e y:
x- matrix(rnorm(25, 10, 4),nrow = 5, ncol = 5); x #matrix 1
y- matrix(rnorm(45, 50, 10), nrow = 9, ncol = 5); y # matrix 2
alguem poderia me ajudar em
Ae pessoal, eu qro saber se alguem pode me dar uma dica na seguinte:
##
anos.observado- c(1991:2000)
anos.estimado- c(1994:2000)
x.estimado-c(20,10,13,12, 21,24, 20)
x.observado- c(1,2,3,NA,5,6,NA,8,NA,10)
##
Eu quero preencher os valores observados 'NA', com os estimados, porem a
seq. de anos
Crie um data.frame para as variáveis observadas e outro para as estimadas.
Depois, faça um merge dos dois e sobrescreva a coluna de interesse após o
merge.
On Jul 4, 2014 4:04 PM, Leonardo Monteiro monteiroleonar...@gmail.com
wrote:
Ae pessoal, eu qro saber se alguem pode me dar uma dica na
Não conhecia esse testador de expressões regulares.
Muito bom!
Abraços,
Paulo Nogueira Starzynski
Em 3 de janeiro de 2014 16:01, Rodrigo Coster rcos...@gmail.com escreveu:
Não sei se são as soluções mais elegantes, mas seguem as minhas sugestões:
# Restringindo o conteúdo do parenteses:
Ola a todos, e feliz ano novo.
Suponha que eu tenho o seguinte exemplo:
#
exemplo-c(Bursey et al. (2001), Bursey et al. (2001), Bursey et al.
(2001),
Bursey et al. (2001), Martins and Fabio (2005), Martins and Fabio
(2005),
Martins and Fabio (2005), Vicente et al. (1991), Vicente et al.
(1991);
Não sei se são as soluções mais elegantes, mas seguem as minhas sugestões:
# Restringindo o conteúdo do parenteses:
grep('\\([0-9]{4}.?; ?[0-9]{4}.?\\)', exemplo, value=T) # Explicação:
http://regex101.com/r/nB6jH4
# Não restringindo o conteúdo do parenteses:
grep('\\([^\\)]*;.*\\)', exemplo,
Se tivesse a posteriori de um parâmetro dava para fazer um
histograma/densityplot e pegar intervalos com a função quantile. Como o
caso é bidimensional, tem-se que fazer uma densidade kernel em duas
dimensões. Para isso veja o help da função::kde2d() (
Ou, para considerar correcao de borda, usar kernel2d() do pacote
'splacs'. Particularmente util para espaco parametrico limitado com
estimativas proximas 'a borda.
On 11/04/2013 05:19 PM, walmes . wrote:
Se tivesse a posteriori de um parâmetro dava para fazer um
histograma/densityplot e pegar
Windows XP
@Mauro Sznelwar
Qual sistema operacional você esta usando?
Parece que o erro tem haver de estar faltando alguma coisa, o compilador de
c++?
Eu uso o ubunto aqui, e ja veio junto o compilador, o compilador está
instalado ae certinho?
Na verdade eu não entendo muito bem de
Pelo que eu entendo, o Rcpp compila o codigo que você escreve em c++ numa
função, então tem que ter o compilador, no linux ele usa o gcc, agora no
windows não sei, mas imagino que esteja faltando o compilador. De uma
olhada nesse post:
@Mauro Sznelwar
Qual sistema operacional você esta usando?
Parece que o erro tem haver de estar faltando alguma coisa, o compilador de
c++?
Eu uso o ubunto aqui, e ja veio junto o compilador, o compilador está
instalado ae certinho?
Na verdade eu não entendo muito bem de rcpp, estou apenas
Como eu faço para rodar este código? Veja a mensagem que recebi:
Error in sourceCpp(code = code, env = env, rebuild = rebuild, showOutput =
showOutput, :
Error 1 occurred building shared library.
WARNING: The tools required to build C++ code for R were not found.
Please download and
vetor.erase(n)?
On Oct 17, 2013 4:39 PM, Augusto Ribas ribas@gmail.com wrote:
Ola pessoal.
Eu estou com uma dificuldade em como fazer uma chamada recursiva usando o
Rcpp.
Por curiosidade, eu estava fazendo alguns códigos aqui usando o Rcpp, só
por experiencia.
Ai eu estava
Prezados
Estou tentando concatenar uma variável com a seguinte condição:
Se ela tiver tamanho 1, gostaria de concatenar 000, se ela tiver tamanho 2
concateno 00 e se ela tiver tamanho 3 concateno 0.
Eu pensei na seguinte função:
dados$len - nchar(dados$var)
func - function(){
if (dados$len ==
Se a variável for inteiro
sprintf(%04d, c(1,3,30,50,900))
On 07/22/2013 08:57 AM, Sérgio Henrique almeida da silva ju wrote:
Prezados
Estou tentando concatenar uma variável com a seguinte condição:
Se ela tiver tamanho 1, gostaria de concatenar 000, se ela tiver
tamanho 2 concateno 00 e se
Não entendi Elias!
Em 22 de julho de 2013 09:00, Elias Krainski
eliaskrain...@yahoo.com.brescreveu:
Se a variável for inteiro
sprintf(%04d, c(1,3,30,50,900))
On 07/22/2013 08:57 AM, Sérgio Henrique almeida da silva ju wrote:
Prezados
Estou tentando concatenar uma variável com a
sprintf(), usando '%04d', faz o que vc quer para inteiros.
Se dados$var é 'character', use
aaa - sapply(3:0, function(n) paste(rep(0, n), collapse=''))
dados$var1 - paste(aaa[nchar(dados$var)], dados$var, collapse='')
___
R-br mailing list
if (dados$len == 1){
A sua solucao nao funciona porque essa condicao tem comprimento igual ao
tamanho de dados$len, mas if() precisa de TRUE ou FALSE (tamanho 1).
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
Elias
Não funcionou!
O comando criou o resultado abaixo, tudo dentro da mesma variável.
000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
1000
Em 22 de julho de 2013 09:22, Elias Krainski
eliaskrain...@yahoo.com.brescreveu:
if (dados$len == 1){
A sua solucao nao
Sergio,
da um dput(head(dados$var)) e nos cola o que aparecer na tela
2013/7/22 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com
Elias
Não funcionou!
O comando criou o resultado abaixo, tudo dentro da mesma variável.
000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
dput(head(dados$Filial))
c(1L, 1L, 1L, 1L, 1L, 1L)
Exemplo a variável pode ter os seguintes valores: 0, 1, 01, 44, 456, 0001
e eu quero que ela fique assim: , 0001, 0001, 0044, 0456, 0001.
Abraços
Em 22 de julho de 2013 09:39, Rodrigo Coster rcos...@gmail.com escreveu:
Sergio,
da um
Entao o que tu precisa fazer é dados$Filial - sprintf('%04d', dados$Filial)
2013/7/22 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com
dput(head(dados$Filial))
c(1L, 1L, 1L, 1L, 1L, 1L)
Exemplo a variável pode ter os seguintes valores: 0, 1, 01, 44, 456, 0001
e eu quero que
Funcionou amigo!
Obrigado.
Em 22 de julho de 2013 09:48, Rodrigo Coster rcos...@gmail.com escreveu:
Entao o que tu precisa fazer é dados$Filial - sprintf('%04d',
dados$Filial)
2013/7/22 Sérgio Henrique almeida da silva ju sergio.edfis...@gmail.com
dput(head(dados$Filial))
c(1L, 1L,
Caro Colegas,
Bom dia.
Eu estou tentando construir uma PCA no base, mas estou entendo alguns problema.
Tenho uma data.frame com uma célula com valor #F7F7F7 (dm.x[,1+i]) e gostaria
de substitui-la por outro valor de acordo com o valor de uma célula especifica,
assim tenho:
dificil, muito dificil, ser preciso na ausencia de um exemplo reproduzivel...
mas, acho que o problema resume-se 'a falta de aspas em torno do valor...
blah - '#D9D9D9'
b
Em 22 de maio de 2013 10:16, Alisson Lucrecio
alisson...@yahoo.com.br escreveu:
Caro Colegas,
Bom dia.
Eu estou
: [R-br] Duvida plot base1
dificil, muito dificil, ser preciso na ausencia de um exemplo reproduzivel...
mas, acho que o problema resume-se 'a falta de aspas em torno do valor...
blah - '#D9D9D9'
b
Em 22 de maio de 2013 10:16, Alisson Lucrecio
alisson...@yahoo.com.br escreveu:
Caro Colegas
@listas.c3sl.ufpr.br; Alisson Lucrecio
alisson...@yahoo.com.br
Sent: Wednesday, May 22, 2013 10:26 AM
Subject: Re: [R-br] Duvida plot base1
dificil, muito dificil, ser preciso na ausencia de um exemplo
reproduzivel...
mas, acho que o problema resume-se 'a falta de aspas em torno do valor...
blah
Costa
From: Benilton Carvalho beniltoncarva...@gmail.com
To: Alisson Lucrecio alisson...@yahoo.com.br
Cc: r-br r-br@listas.c3sl.ufpr.br
Sent: Wednesday, May 22, 2013 12:36 PM
Subject: Re: [R-br] Duvida plot base1
Alisson,
o que exatamente vc quer fazer e
@listas.c3sl.ufpr.br
Sent: Wednesday, May 22, 2013 12:36 PM
Subject: Re: [R-br] Duvida plot base1
Alisson,
o que exatamente vc quer fazer e' substituir a coluna 'col' de dfm
pela coluna 'right' de var? E' isso?
pq nao usar merge?
saida = merge(dfm, var, by.x='Solo', by.y='wrong')
saida[[2]] = NULL
: [R-br] Duvida plot base1
novamente exemplo reproduzivel.
bom trabalho
Em 22 de maio de 2013 12:49, Alisson Lucrecio
alisson...@yahoo.com.br escreveu:
O exemplo mostrei coincidiu de ser 2 e 2, mas o meu problema é com um numero
diferente (ex. dfm com 36 lina e na var ex, 4). Assim teria que
Cc: r-br r-br@listas.c3sl.ufpr.br
Sent: Wednesday, May 22, 2013 1:27 PM
Subject: Re: [R-br] Duvida plot base1
novamente exemplo reproduzivel.
bom trabalho
Em 22 de maio de 2013 12:49, Alisson Lucrecio
alisson...@yahoo.com.br escreveu:
O exemplo mostrei coincidiu de ser 2 e 2, mas o
Muito obrigado.
Alisson Lucrécio da Costa
From: Benilton Carvalho beniltoncarva...@gmail.com
To: Alisson Lucrecio alisson...@yahoo.com.br
Cc: r-br r-br@listas.c3sl.ufpr.br
Sent: Wednesday, May 22, 2013 2:35 PM
Subject: Re: [R-br] Duvida plot base1
jura q
Ola, eu gostaria de saber, se tem como eu mudar de forma simples a
mensagem de erro que o R envia ao tentar usar uma função e receber um
erro.
Por exemplo:
#Imagine que eu tenho a seguinte lista, e precise aplicar uma função
#Exemplo de lista de dados
teste-list(2, a, 1 , 10 , b)
teste
Da para fazer com ifelse:
ifelse(sapply(saida, class) == 'try-error', 'Verificar', saida)...
Se a ideia é só verificar se deu certo, tu pode fazer direto com o class()
dentro de um if/for...
2013/5/9 Augusto Ribas ribas@gmail.com
Ola, eu gostaria de saber, se tem como eu mudar de forma
Augusto,
Complicado ein, no seu exemplo você deu uma entrada de um tipo diferente.
log(-1) não existe mas o R ainda retorna um valor indicando erro (NaN).
A sua do taxize retorna qual valor em caso de falha?
Abraços
Em 10 de maio de 2013 14:39, Rodrigo Coster rcos...@gmail.com escreveu:
Da
Tente assim:
sapply(teste, function(x)tryCatch(log(x), error = function(...)verificar))
2013/5/9 Augusto Ribas ribas@gmail.com
Ola, eu gostaria de saber, se tem como eu mudar de forma simples a
mensagem de erro que o R envia ao tentar usar uma função e receber um
erro.
Por exemplo:
#
sapply(teste, function(x)tryCatch(log(x), error = function(...)verificar))
#
Isso era o que eu queria :), funcionou perfeitamente.
David, o taxisize retorna uma mensagem de erro que é uma String. Por
exemplo eu guardava os resultados numa matriz de characters, onde uma
coluna é o nome que usei
Bom que o problema se resolveu com os exemplos acima dos colegas da lista.
Checagem de erro no R fora da biblioteca padrão é um caos mesmo. Poderia
muito bem retornar algo diferente de string (NA por exemplo) e você ia
checar e substituir como foi sugerido acima.
Abraços
Em 10 de maio de 2013
Sim, se saise um NA, um Null seria perfeito.
Eu estava sugerindo no fórum oficial dos pacotes uma forma mais amigável de
testar um vetor ou lista de nomes de forma mais simples.
Agora com exemplos acho que consigo propor algo mais amigável e razoável.
Mas na pagina do try tinha uma menção ao
Eles tem um github (taxize) envie um Pull request na função.
Em 10 de maio de 2013 17:12, Augusto Ribas ribas@gmail.com escreveu:
Sim, se saise um NA, um Null seria perfeito.
Eu estava sugerindo no fórum oficial dos pacotes uma forma
mais amigável de testar um vetor ou lista de nomes de
Esse tipo de atribuição não vai dar certo mesmo
options(contrasts)$contrasts[2]-contr.sum
Você deve fazer assim
options(contrasts=c(unordered=contr.sum, ordered=contr.poly))
Com isso você faz uma alteração global para toda sua sessão R, ou seja,
toda lm(), aov(), glm(), lme(), etc, todas
Por favor, perdoe minha falta de atenção, confundi tudo no script anterior.
#exemplo
set.seed(123)
tratamento-factor(rep(letters[1:3],each=10))
resposta-rnorm(30,rep(c(1,1,4),each=10))
#modelo usando opções default
modelo1-lm(resposta~tratamento)
summary(modelo1)
modelo2-
O problema do seu código foi só a ordem das colunas, veja
set.seed(123)
tratamento - factor(rep(letters[1:3],each=10))
resposta - rnorm(30,rep(c(1,1,4),each=10))
modelo1 - lm(resposta~tratamento)
summary(modelo1)
modelo2 - lm(resposta~tratamento,
Muito obrigado pela explicação, eu que não estou sabendo montar as
comparações então, com esses exemplos eu acho que entendi melhor.
2013/2/7 Walmes Zeviani walmeszevi...@gmail.com:
O problema do seu código foi só a ordem das colunas, veja
set.seed(123)
tratamento -
Ola caros colegas da lista, eu estava tentando estabelecer a matriz de
contrastes para usar na função lm, mas não estou conseguindo entender
como as coisas funcionam no R.
Eu li um tutorial antigo, esse aqui
Pessoal,
estou tentando ler um txt com
teste=read.table.ffdf(filename, first.rows=2, next.rows=10,
VERBOSE=TRUE,colClasses=c(character, character, character,character,
character,character,character,character, character,POSIXct,
POSIXct,character)
)
o filename possui o caminho com o nome do
Pessoal mais uma vez procurando ajuda de vc's, tenho um gigantesco banco de
dados, e minhas variáveis são: pmc, te, tof e quero ver como elas se
comportaram durante o tempo, ou seja por mês e por ano e por mat. So que no
momento de fazer um barplot, os meses das barras estão do mesmo jeito que
nos
Eu não entendi exatamente o problema, mas deve ser algo assim...
Aqui eu tentei fazer um exemplo pra exemplificar. Veja se é isso.
Mas leia o comentário do Walmes antes daqui :)
#criando uma sequencia de datas no R
datas-seq(as.Date(2000/1/1), as.Date(2005/1/1), months)
#colocando elas em um
Bom dia prezados,
Estou finalizando minha dissertação e estou com um problema se alguem puder
ajudar fico muito grata.
Tenho uma serie temporal de temperatura para essa ajustei um modelo que tem
uma diferença sazonal
fiz as previsoes 6 passos a frente, no entanto minhas previsoes saem na
escala
Augusto,
Para você customizar a aparência do seu R veja o post:
http://ridiculas.wordpress.com/2012/05/10/emacs-personalizado/
att,
FH
2012/11/1 Augusto Ribas ribas@gmail.com
Muito obrigados a todos, funcionou perfeitamente :)
Eu achei diferente, por que no windows o help em html é o
Muito obrigados a todos, funcionou perfeitamente :)
Eu achei diferente, por que no windows o help em html é o default, e no
linux é diferente.
Mas criei o Rprofile e coloquei os comandos e agora esta como eu gostaria.
Em 1 de novembro de 2012 01:03, Fernando Mayer
Ola colegas, recentemente começei a usar o ubunto e estou tendo algumas
dificuldades básicas.
Estou tentando usar o emacs como editor para o R e minha duvida é de como
fazer para as paginas de documentação, quando se usa o comando ?função
aparecam no navegador.
Eu não entendi se isso tem que ser
help(plot, help_type=html)
Você pode controlar globalmente alterando o options(), algo assim
options(help_type=html)
help(lm)
À disposição.
Walmes.
==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação,
Outra sugestão é iniciar o help do browser e depois fazer as pesquisas por
lá.
Use:
help.start()
E depois clique em: Search Engine
Keywordshttp://127.0.0.1:28741/doc/html/Search.html
Abs,
D
2012/10/31 Walmes Zeviani walmeszevi...@gmail.com
help(plot, help_type=html)
Você pode
O .Rprofile não é criado automaticamente. Você deve criá-lo se quiser
customizar algumas opções do R, como carregar sempre os mesmos
pacotes, ajustar opções, etc.
O .Rprofile criado no seu $HOME será global, ou seja, em qualquer
diretório que vc iniciar o R, ele irá ler o $HOME/.Rprofile.
Mas vc
Walmes,
não conhecia este uso pro match, vc acabou de me poupar algumas dezenas
de if em uma recodificação que eu estou fazendo aqui
valeu, hehehe
[]s
Leonard de Assis
assis dot leonard at gmail dot com
Em 31/08/2012 23:12, Walmes Zeviani escreveu:
Aqui vai uma dica para eliminar os for() que
Leonardo, eu fazia vários if() para recodificar e fazer outras operações.
Quando descobri o match() fique de cara. Augusto, esse seu conjunto de
dados e problema exposto é muito interessante. Qual a possibilidade de você
disponibilizar total ou parcialmente os dados para a gente documentar as
Pode usar os dados como exemplo sim, vou mandar um e-mail com os dados
organizados de forma organizada pra vc. É um prazer poder contribuir.
Em 1 de setembro de 2012 14:19, Leandro Marino
leandromar...@leandromarino.com.br escreveu:
*Este é o motivo pelo qual leio todos os tópicos desta
Ola pessoal.
Eu estava olhando aqui uns dados com um amigo e me surgiu uma duvida.
O problema é o seguinte:
Os dados são referentes a ocorrência de uma doença em passarinhos. Pegamos
muitos passarinhos em redes em um lugar e testamos se ele esta doente ou
não.
Durante a captura, registramos
Augusto,
Tente colocar NA (missing value) onde você não tem a informação.
Att,
Heloise
Em 31 de agosto de 2012 16:18, Augusto Ribas ribas@gmail.com escreveu:
Ola pessoal.
Eu estava olhando aqui uns dados com um amigo e me surgiu uma duvida.
O problema é o seguinte:
Os dados são
Então, mas ai:
na.action
a function which indicates what should happen when the data contain NAs.
The default is set by the na.action setting of
optionshttp://127.0.0.1:30387/library/stats/help/options,
and is na.fail http://127.0.0.1:30387/library/stats/help/na.fail if that
is unset. The
Arrumando o exemplo: (eu não tinha colocado pra usar o dataframe que criei
nas analises)
#dados, exemplo
set.seed(666)
n-90
idade-sample(c(Adulto,Sub-Adulto,Jovem),n,replace=T)
idade
efeito-ifelse(idade==Jovem,1,0)
efeito
diagnostico-rbinom(n,1,0.75+efeito*(-0.25))
sexo-NA
for(i in
Tente:
modelo01-glm(as.factor(diagnostico)~factor(sexo)+factor(idade),family=
binomial,data=dados,na.action=na.exclude)
Em 31 de agosto de 2012 17:32, Augusto Ribas ribas@gmail.com escreveu:
Arrumando o exemplo: (eu não tinha colocado pra usar o dataframe que
criei nas analises)
Então:
Ai eu obtenho:
Call:
glm(formula = as.factor(diagnostico) ~ factor(sexo) + factor(idade),
family = binomial, data = dados, na.action = na.exclude)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.8176 0.4892 0.6524 0.7731 0.7833
Coefficients:
Augusto, se o que te interessa é idade, como você disse no primeiro e-mail,
minha sugestão é que uses apenas idade no modelo. O sexo sozinho tem efeito
significativo?
Em 31 de agosto de 2012 18:01, Augusto Ribas ribas@gmail.com escreveu:
Então:
Ai eu obtenho:
Call:
glm(formula =
Então nos dados não. Mas por exemplo, machos podem investir em plumagem,
entre outras coisas, deixando o sistema imuno mais fraco(menos
investimento). O que propiciaria uma maior chance de infecção.
Então eu queria testar as duas coisas (sexo e idade), mas me parece que
dado esse problema de
Ainda é possível você assumir que o efeito da classe não sei seja uma
média das outras duas. Isso é razoável se os dois sexos ocorrem na mesma
proporção e se sua inabilidade de identificar o sexo não se altere com a
idade de forma diferencial. Ou seja, na idade 1 mês existe relação 1:1 dos
sexos
Nossa, parece que funciona muito bem. Eu fiz um exemplo pra ver se estou
entendendo bem:
#gerando dados de exemplo
#variaveis resposta
set.seed(3)
n-300
sexo-sample(c(macho,femea),n,replace=T)
idade-sample(c(jovem,subadulto,adulto),n,replace=T)
#como é o efeito delas
Aqui vai uma dica para eliminar os for() que você usou, pode-se usar o
match() para fazer essas recodificações. Deve ter uma forma ainda melhor,
mas atualmente eu uso essa.
sexo - as.factor(sample(c(M,F,I), 10, repl=TRUE))
recod - c(-1,0,1)
sexo.rec - recod[match(sexo, c(F,I,M))]
data.frame(sexo,
Ficou muito bom assim :)
Exatamente o que eu precisava.
Obrigado Walmes e Heloise pela atenção.
Em 31 de agosto de 2012 22:12, Walmes Zeviani walmeszevi...@gmail.comescreveu:
Aqui vai uma dica para eliminar os for() que você usou, pode-se usar o
match() para fazer essas recodificações. Deve
Enviadas: Terça-feira, 14 de Agosto de 2012 15:03
Assunto: Re: [R-br] Duvida no TAPPLY
use addmargins
Jonas
Em 14 de agosto de 2012 10:38, Edson Lira edinhoes...@yahoo.com.br escreveu:
Estou usando esta rotina e obtenho
tapply(acj$Quantidade,list(Laboratório=acj$lab,Mês=acj$Mês),sum,na.rm
Estou usando esta rotina e obtenho
tapply(acj$Quantidade,list(Laboratório=acj$lab,Mês=acj$Mês),sum,na.rm=T)#Laboratórios
E obtenho esta saída abaixo, que me dá o quantitativo de exames mensal, para
obter a soma nas linhas e a soma nas colunas,o que devo acrescdentar a rotina?
?rowSums, ?colSums
\begin{signature}
=
Prof. Dr. Ivan Bezerra Allaman
Universidade Estadual de Santa Cruz
Departamento de Ciências Exatas e Tecnológicas
Ilhéus/BA - Brasil
Fone: +55 73 3680-5596
E-mail: ivanala...@yahoo.com.br/ivanala...@gmail.com
@
1 - 100 de 194 matches
Mail list logo