Re: [R-br] Dados ausaentes (NA)

2018-08-02 Por tôpico Cesar Rabak via R-br
Edson, Estou "de passagem" e lendo de forma perfunctória todos os posts, por isso *pode ser* que o quê vou sugerir seja apenas outro idioma R para o que você postou, mas para eliminar registros de um dataframe que contenha NA o que estou acostumado a ver, já usando o nome do seu, seria: > b <- b[

Re: [R-br] Curvas não lineares - qual o melhor modelo?

2018-08-10 Por tôpico Cesar Rabak via R-br
Thallita, Uma visão complementar no que tange à escolha de modelos, é ver qual explicaria o fenômeno melhor e aí ver qual o ajuste dos dados a ele. SE você fizer dessa forma o modelo poderia ser usado para extrapolar dados ("fazer predições") enquanto um ajuste (por mais "adequado" que possa pare

Re: [R-br] Diferenças entre sistemas operacionais (geoR)

2018-10-04 Por tôpico Cesar Rabak via R-br
Para esmagar *esse* bug, pode-se setar uma semente específica em ambos os SO e ver se aí há alguma diferença nos resultados. On Thu, Oct 4, 2018 at 5:45 PM Thiago V. dos Santos via R-br < r-br@listas.c3sl.ufpr.br> wrote: > Da maneira que eu vejo, a unica explicacao plausivel seria se a krige.conv

Re: [R-br] Tamanho de Gráficos

2016-07-10 Por tôpico Cesar Rabak via R-br
Pedro, Via de regra o *publisher* estabelece isso nas suas regras para publicação, seja num manual seja num documento exemplo ou no modelo. Tem a ver com a permissão da impressão ser em uma ou duas colunas e (e às vezes se é possível certas figuras para ocupar a página inteira [eventualmente até "

Re: [R-br] Ajustes do comando pdf

2016-07-20 Por tôpico Cesar Rabak via R-br
Elias, Acho que a forma mais fácil seria você ver qual dos comandos da família dev.<> força um "nova página" no seu dispositivo PDF, ou alternativamente qual opção em semPlot::semPath() faria essa função. HTH 2016-07-13 0:17 GMT-03:00 Elias Carvalho via R-br : > Pessoal > > Com o código abaixo,

Re: [R-br] Levene Test em experimentos fatoriais

2016-08-07 Por tôpico Cesar Rabak via R-br
Maurício, Você deve proceder *não fazendo* esse teste a priori de qualquer outro de inferência estatística: 1) Além dos problemas de esses testes serem ou pouco sensíveis se a amostra for pequena e demasiado sensíveis (sempre acabam rejeitando a hipótese nula) se a amostra for "grande"; 2) eles eq

Re: [R-br] Matriz distâncias rodoviárias

2016-08-09 Por tôpico Cesar Rabak via R-br
Rafael, Esse problema (de obter alguma base de dados com as distâncias "rodoviárias") é recorrente em vários *fora*, geralmente de programação ou ligados a otimização, etc. O busilis está em obter esses dados já em uma base, pois hodiernamente todo mundo que "disponibiliza" essa informa (leia-se

Re: [R-br] Matriz distâncias rodoviárias

2016-08-10 Por tôpico Cesar Rabak via R-br
lto:r-br-boun...@listas.c3sl.ufpr.br] *Em nome de *Cesar > Rabak via R-br > *Enviada em:* terça-feira, 9 de agosto de 2016 15:25 > *Para:* Raphael Saldanha ; a lista Brasileira > oficial de discussão do programa R. > *Assunto:* Re: [R-br] Matriz distâncias rodoviárias > > >

Re: [R-br] ISO8601 time duration

2016-08-10 Por tôpico Cesar Rabak via R-br
Raphael, Essa funcionalidade está "estacionada" em um pedido de "feature request" para um pacote R que cuida exatamente dessas coisas lubridate desde dezembro de 2015! O *parsing* de cadeias ISO 8601 não é trivial e é bem documentado em *fora* de

Re: [R-br] Digest R-br, volume 68, assunto 17

2016-08-19 Por tôpico Cesar Rabak via R-br
Ana, No R para ver a "aderência" a uma distribuição qualquer usa-se o teste de Kolmogorov-Smirnov. no seu terminal do R tecle: > ks.test? Procure nesta lista as razões porque você *não *deve fazer isso. HTH -- Cesar Rabak 2016-08-19 9:21 GMT-03:00 ana paula coelho madeira via R-br < r-br@lista

Re: [R-br] Teste T-pareado

2016-08-22 Por tôpico Cesar Rabak via R-br
Olá Simone! Como outros listeiros já resolveram o problema "profilático" do seu exemplo, e dado a responsabilidade dele por ser empregado em uma aula de estatística, deixo aqui para você o seguinte "besouro atrás da orelha": Por que esse exemplo (desse tipo de dados) viola uma das premissas do te

Re: [R-br] Teste T-pareado

2016-08-22 Por tôpico Cesar Rabak via R-br
NÃO, Emerso, a diferença entre os "paired" nas duas chamadas deve-se que no primeiro caso se usa um par de vetores de dados e no segundo um único vetor, e portanto o "paired" não tem sentido nenhum... Veja as respostas mais acima para ver a razão real da "discrepância". 2016-08-22 14:15 GMT-03:

Re: [R-br] Teste de Aderência

2016-08-22 Por tôpico Cesar Rabak via R-br
Ana, Antes de mais, nada, recomendo (e eu não estou sozinho nessa) que você não faça testes de aderência em distribuições de dados obtidas de experimentos, especialmente se elas forem servir para algum teste inferencial posterior. Com isso posto, o teste ks.test() do R pode verificar a aderência

Re: [R-br] Avaliação da evolução do peso

2016-08-31 Por tôpico Cesar Rabak via R-br
Já tentou o pacote dplyr? 2016-08-31 10:43 GMT-03:00 Edson Lira via R-br : > Eu gostaria de fazer que a estrutura do banco ficasse assim: > > id sexo dt_nasc dt1_acom dt2acom dt3_acom peso1 peso2 peso3... > > [ ]'s > Edson Lira > Estatístico > Manaus-Amazonas > > > Em Terça-feira, 30 de Agosto

Re: [R-br] Parcela Subdividida

2016-09-05 Por tôpico Cesar Rabak via R-br
Angélica, Você entendeu a mensagem de erro? Não é possível analisar uma matriz nula 2016-09-04 1:22 GMT-03:00 Angélica Ricarte via R-br < r-br@listas.c3sl.ufpr.br>: > Prezados, > > Estou tentando realizar uma análise de variância para um experimentos em > parcelas subdivididas em DIC, mas

Re: [R-br] Parcela Subdividida

2016-09-06 Por tôpico Cesar Rabak via R-br
Angélica, Como você deu nomes "crípticos" para seu dataframe, e não postou um CMR, dá só para fazer uma 'interpretação' das mensagens. Aparentemente você quer a resposta PS[,4] para a interação entre as variáveis (grupos[?]) PS[,1] e PS[,2] e considera a combinação dos fatores PS[,1] e PS[,3] com

Re: [R-br] (sem assunto)

2016-09-07 Por tôpico Cesar Rabak via R-br
Elias, Se você exportou para docx como propôs o Leonard, a tabela é um objeto Word e portanto as larguras e outros detalhes de formatação são manipuláveis no docx, não? Se não é esse o caso, por favor especifique um pouco melhor o quê exatamente o "assim" significa! Melhor ainda, usando um dos c

Re: [R-br] Parâmetros do plot

2016-09-07 Por tôpico Cesar Rabak via R-br
Você viu a documentação do plot.bpca2d no help do pacote? 2016-09-07 8:34 GMT-03:00 Andre Oliveira via R-br : > Bom dia, > pesquisei e conversei com algumas pessoas e não resolvi o meu problema. > Como tirar estes vetores desta imagem feita com o plot? > require(bpca) > fit=bpca > plot(bpca(ir

Re: [R-br] Passando um objeto para "character"

2016-09-09 Por tôpico Cesar Rabak via R-br
Diego, A ideia parece boa, mas apenas para uma função para uso "interno" ou "pessoal". . . Uma função de uso geral com um *attach* escondido dentro pode gerar uma série de erros inesperados cuja elucidação por parte do emissor do pacote (Pedro Emmanuel) será um pesadelo. A opção do with() é a ma

Re: [R-br] Criação de um Fluxograma

2016-09-14 Por tôpico Cesar Rabak via R-br
Caro Abel, Sem muita elaboração, me parece que o quê você deseja se aproxima da proposta do pacote data.tree. Veja a documentação e os exemplos e se você consegue: ou adequá-lo às suas necessidades ou você se adequar às características dele!! HTH 2016-09-14 11:28 GMT-03:00 Abel Brasil Ramos da

Re: [R-br] Distribuição exponencial mista

2016-09-19 Por tôpico Cesar Rabak via R-br
. . . e sua dúvida seria? On Mon, Sep 19, 2016 at 10:54 AM, Alisson Lucrécio via R-br < r-br@listas.c3sl.ufpr.br> wrote: > Caro Colegas, > > Bom dia. > > Algum colega da poderia me ajudar na solução de uma distribução > exponencial mista. > > # Script > > library(MASS) > > # Distribuição Gamma >

Re: [R-br] Usar R para Salvar Referencias

2016-09-19 Por tôpico Cesar Rabak via R-br
Teoricamente possível, mas dado que o R não é por excelência uma linguagem para scripting para esse tipo de coisas (Internet), ao menos que a geração das referências a serem buscadas fosse originada em algum procedimento do R a coisa mais lúcida IMNSHO seria fazer isso via alguma outra linguagem de

Re: [R-br] Salvar uma tabela/matriz como imagem

2016-09-20 Por tôpico Cesar Rabak via R-br
Elias, Você mencionou explicitamente a palavra *imagem*, por isso acho que a resposta do Jônatan lhe será mais apropriada do que este meu comentário, mas mesmo assim talvez lhe ajude. Considerando que uma tabela é por excelência um objeto mais voltado para texto que para imagem (inclusive do pont

Re: [R-br] Precipitação diária baseada em uma rede neural

2016-09-25 Por tôpico Cesar Rabak via R-br
Alisson, O uso de redes neurais brilha num quadro onde se tenha uma massa enorme de dados e algum tipo de "resposta" frente a variáveis explicativas e não se tenha nenhuma ideia de uma modelagem mais descritiva e o mais importante após o "treinamento" da rede neural a explicação seja menos importa

Re: [R-br] Regexp no grep para linhas que não iniciam com sequência

2016-10-04 Por tôpico Cesar Rabak via R-br
Marcus, De fato existe o PowerShell da Microsoft, mas ele não inclui a funcionalidade dos comandos Unix like, como, por exemplo, o grep que o Rafael usou. No Windows 10 (IIUC somente no dez) você pode instalar o "ambiente bash" que inclui em adição ao próprio bash, o grep, sed, awk (pelo menos se

Re: [R-br] Regressão linear - Excluindo outliers

2016-10-14 Por tôpico Cesar Rabak via R-br
FCosta, Você viu o posto do Valmes sobre o influence.measures()? 2016-10-14 17:45 GMT-03:00 Felinto COSTA via R-br : > Antes do estabelecimento de critérios empíricos de corte para os valores > da distância de Cook, > há esse trabalho sobre sua distribuição exata: > > https://www.researchgate.

Re: [R-br] boxplot comparativos

2016-10-15 Por tôpico Cesar Rabak via R-br
Elisa, Veja se estes exemplos a "inspiram". no seu ambiente R, tente example(boxplot) e veja os exemplos um por um, mais de um deles mostram boxplots comparativos. Pelo que você escreveu, parece-me que o exemplo que mais se aproxima do seu objetivo é do crescimento de dentes dos porquinhos da Índi

Re: [R-br] comparação de médias - problema com os pressupostos

2016-11-02 Por tôpico Cesar Rabak via R-br
Maurício, Toda vez que você tem "comparação de médias" e "pressupostos" é necessário pensar no objetivo da comparação: - Você está usando um teste de médias para testar parametricamente se as distribuições podem ser consideradas iguais ou não, ou; - você está querendo testar as médias e

Re: [R-br] Teste de independência dos resíduos

2016-11-02 Por tôpico Cesar Rabak via R-br
Mauricio, Em adição às indicações do Walmes, pesquise nesta lista outras trédis sobre esse assunto onde eu comento a precariedade desses testes e necessidade de você ter mais segurança quanto ao processo que gerou os dados. HTH 2016-10-27 0:36 GMT-02:00 Maurício Sangiogo via R-br < r-br@listas.c

Re: [R-br] RES: mutate vs within

2016-11-02 Por tôpico Cesar Rabak via R-br
Leonard, Seu problema me parece que seria muito bem atendido ("código limpo", "claro", etc.) pela função recode do pacote car. Como parece-me que você tem bastante coisa para converter talvez a curva de aprendizado dessa nova função possa ser um investimento frente ao possível resultado. HTH --

Re: [R-br] Sorteios aleatórios de um cadastro (AAS ou amostra estratificada)

2016-11-02 Por tôpico Cesar Rabak via R-br
André, Como você passou pouca informação, posso sugerir-lhe o pacote survey a função stratsample() para a estratificada, enquanto por óbvio, a AAS você pode realizá-la com facilidade seja com sample() (como o Walmes sugeriu) ou gerando um conjunto de números aleatórios e pós-processando-os caso a

Re: [R-br] Heterocedasticidade e teste de média

2016-11-08 Por tôpico Cesar Rabak via R-br
Luiz, Malgrado você tenha endereçado a pergunta, dou meu "pitaco" enquanto os outros colegas não se manifestam. Todas as regressões necessitam que em última análise os erros aleatórios possam ser modelados por um processo cuja melhor representação é a distribuição gaussiana com média zero e algu

Re: [R-br] PCA

2016-11-09 Por tôpico Cesar Rabak via R-br
Ana, Além das funções citadas, recomendo que você veja os pacotes FactoMineR http://factominer.free.fr/ e ade4 http://pbil.univ-lyon1.fr/ade4/home.php?lang=eng. Em adição aos pacotes propriamente ditos (que você pode baixar diretamente dos sites CRAN no ambiente R) nos sites deles há uma extensa

Re: [R-br] Eliminar observações iguais em variáveis distintas

2016-11-11 Por tôpico Cesar Rabak via R-br
Se eu entendi corretamente, você precisa "forçar" os níveis nas duas variáveis a serem "os mesmos", por isso você precisa fazer: levels(expand$ido) <- levels(expand$idd) NÃO TESTADO! Verifique com cópia dos seus dados. . . HTH -- Cesar Rabak 2016-11-11 12:34 GMT-02:00 Alexandre Loures via R-br

Re: [R-br] Tamanho de amostras

2016-11-13 Por tôpico Cesar Rabak via R-br
Existem vários. Para ver qual se adéqua melhor às suas necessidades dê uma pesquisada em CRAN Task Views . HTH -- Cesar Rabak 2016-11-13 18:23 GMT-02:00 Emerson Cotta Bodevan via R-br < r-br@listas.c3sl.ufpr.br>: > Pessoal, boa noite. > > Existe algum pac

Re: [R-br] Residuos do Teste T de Student

2016-11-19 Por tôpico Cesar Rabak via R-br
Gracialiano, Por que você esperaria resíduos de um teste t de Student? Você conhece alguma referência ou obra que trate desse teste e fale sobre resíduos calculados devido a esse teste? 2016-11-19 15:55 GMT-02:00 Graciliano Galdino via R-br < r-br@listas.c3sl.ufpr.br>: > Olá a todos, > > Quand

Re: [R-br] sobre a utilização do paste para formar nomes de objetos

2016-11-24 Por tôpico Cesar Rabak via R-br
Antônio, *Acho* que não entendo sua dúvida e sua descrição do que você deseja versus o que você tentou me parece confusa. para fazer a cadeia "matriz" ser concatenada com, digamos, 1 até 12, (o seu "a in 1 to ..." onde o "to" seria em 12), basta apenas fazer o seguinte: paste0("matriz",1:12) HT

Re: [R-br] logica condicional e criar arquivo netCDF

2016-11-25 Por tôpico Cesar Rabak via R-br
Mateus, A lógica que você propõe no seu post é a correta e funcionará se em " casos <- UV[UV>=resultado]" resultado for um escalar. No seu script a linha que parece ser similar ao que você quer é "acima_95 = as.matrix(Con(UV>=Percentile95,UV,Con(UV>=Percentile95, UV, 0))) " onde 'Percentile95' nã

Re: [R-br] Cálculo do coeficiente de assimetria (skewness coefficient) - ecologia de populações

2016-11-30 Por tôpico Cesar Rabak via R-br
Juliana, Como é definido esse coeficiente de assimetria que você usa para prever « a direção de mudança de uma população »? 2016-11-29 10:48 GMT-02:00 Juliana Loureiro via R-br < r-br@listas.c3sl.ufpr.br>: > Olá pessoal! > > Estou com uma dúvida relacionada ao cálculo do coeficiente de assimet

Re: [R-br] Dúvida sobre o gls

2016-11-30 Por tôpico Cesar Rabak via R-br
Luiz, Dê uma olhada em: Modelos DE Regressão com apoio computacional Modelos Lineares Generalizados e Extensões Modelos Lineares Generalizados - da teoria à prática -

Re: [R-br] trabalhando com números pequenos

2016-12-10 Por tôpico Cesar Rabak via R-br
Luiz, Esse problema volta e meia é discutido aqui na lista, você poderá encontrar respostas para essa sua dúvida em posts anteriores. Para uma visão mais completa ("definitiva" segundo alguns) remeto-lhe ao texto do Goldberg: http://dl.acm.org/citation.cfm?id=103163 HTH -- Cesar Rabak 2016-12-

Re: [R-br] trabalhando com números pequenos

2016-12-12 Por tôpico Cesar Rabak via R-br
A solução seria você escrever uma função de teste de Shapiro usando a biblioteca Rmprf! E se você o fizer colocá-la no namespace do pacote Rmprf de tal sorte que o uso dela a ativasse. HTH -- Cesar Rabak PS.: Não entro na seara, neste post, se faz sentido fazer um teste de Shapiro para números d

Re: [R-br] Gráfico

2016-12-14 Por tôpico Cesar Rabak via R-br
Além do excelente pacote do Minato Nakazawa já apontado pelo Maciel, o qual já usei, e posso dizer que o único senão na minha opinião é a necessidade de "massagear" os dados para o formato que a função radarchart desse pacote deseja, meus associados têm usado o pacote "radarchart" cuja versão atual

Re: [R-br] classes de colunas: iris como exemplo

2016-12-19 Por tôpico Cesar Rabak via R-br
Cleber, Minha preferência para ver esse tipo de informação é aplicar a função str() ao data.frame. HTH -- Cesar Rabak On Sun, Dec 18, 2016 at 10:50 AM, Cleber N.Borges via R-br < r-br@listas.c3sl.ufpr.br> wrote: > Prezados, > > tentando detectar quais as classes de colunas de dados me deparei

Re: [R-br] dose-resposta

2016-12-21 Por tôpico Cesar Rabak via R-br
Veja se na Task View de farmacocinética algum dos pacotes indicados não lhe serve. HTH -- Cesar Rabak 2016-12-20 22:26 GMT-02:00 Maurício Lordêlo via R-br < r-br@listas.c3sl.ufpr.br>: > Caros, > Foi realizado um experimento no qual deseja-se saber qual composto > extraído de uma determinada esp

Re: [R-br] Modelo linear com autocorrelação

2016-12-21 Por tôpico Cesar Rabak via R-br
Paulo, Considerando que você está com duas modelagens de processos completamente diferentes, qual é o seu racional para aferir a « . . . diferença relativamente grande nos coeficientes estimados do modelo a partir da inclusão. »? Além da questão mais matemática da análise estatística que você est

Re: [R-br] dose-resposta

2016-12-22 Por tôpico Cesar Rabak via R-br
Walmes, Apenas um observação: não seria mais interessante aproveitar mais a informação da dose se a tratássemos como variável contínua fazendo uma ANCOVA? Minha interpretação apressada do seu código me faz crer que a dose virou fator também. . . []s -- Cesar Rabak On Thu, Dec 22, 2016 at 9:46

Re: [R-br] dose-resposta

2016-12-28 Por tôpico Cesar Rabak via R-br
Achei a solução do Walmes muito interessante. Por > isso, > gostaria de testar a sua abordagem, também! > > Marcelo > > On 22/12/16 at 01:21, Cesar Rabak via R-br wrote: > > Walmes, > > > > Apenas um observação: não seria mais interessante aproveitar mais a > &g

Re: [R-br] R companies

2016-12-29 Por tôpico Cesar Rabak via R-br
Enquanto a lista é alvissareira e nos dá munição para evitar críticas ao tentar usá-lo em algum ambiente que seja mais avesso a software livre, etc., é importante ver que algumas empresas são tão usuárias do R como de todos os outros SW estatísticos, e portanto a lista acaba tendo apenas valor mora

Re: [R-br] R companies

2016-12-30 Por tôpico Cesar Rabak via R-br
onardo Aguirre : > Pra complementar, segue artigo com as principais linguagens de programação > de 2016 (R em 5º): > > http://spectrum.ieee.org/computing/software/the-2016- > top-programming-languages > > Abraços, > > Leonardo > > Em 29 de dezembro de 2016 16:26, C

Re: [R-br] Modelo linear com autocorrelação

2017-01-02 Por tôpico Cesar Rabak via R-br
OK Paulo, Fico feliz que você se "reconciliou" com a diferença dos coeficientes entre as modelagens. Sua abordagem, seria então preparar um modelo retendo uma parte dos dados para teste e usar o restante para um "treinamento" como se faz em certas abordagens mais exploratórias? SE eu entendi cor

Re: [R-br] Extrair Período de datas

2017-01-23 Por tôpico Cesar Rabak via R-br
Salah, Suspeito que o uso da função month seja insuficiente para o Amiko (OP), posto que ele não especifica a faixa de anos que o banco de dados dele abarca, e por outro lado as datas do Carnaval são "feriados móveis". Ademais, como ele fala em "período de Carnaval", cuja definição é bastante fro

Re: [R-br] Analise fatorial

2017-01-23 Por tôpico Cesar Rabak via R-br
Bruce, A mensagem do factanal não permite dizer que "...não pode fazer,", apenas que com os dados e os pontos de partida para o algoritmo ele não conseguiu otimizar. . . veja que este exemplo no SO o indivíduo tinha 10 fatores!! http://stackoverflow.com/questions/30938752/factor-analysis-in-r-err

Re: [R-br] Extrair Período de datas

2017-01-24 Por tôpico Cesar Rabak via R-br
Excelente exemplo Karina! Eu tinha esquecido do lubridate. . . não uso há anos, literalmente. Quanto ao exemplo do salah, a linha: > subset(DF, data > "2016-02-05" & DF$data < "2016-02-10") pode ser (IMHNHO melhor): > subset(DF, data > "2016-02-05" & data < "2016-02-10") Ademais, *acho* que o

Re: [R-br] Importar arquivo para Objeto - Erro

2017-01-24 Por tôpico Cesar Rabak via R-br
Marcelo, Poste aqui as primeiras cinco linhas desse arquivo "pratica1.txt" 2017-01-24 19:43 GMT-02:00 Marcelo Menezes via R-br < r-br@listas.c3sl.ufpr.br>: > Boa Noite pessoal! Sou Marcelo, iniciante no uso do R, ainda estou > aprendendo os comandos básicos a partir do ótimo material disponibil

Re: [R-br] Kruskal-Wallis com comparações entre tratamentos

2017-02-04 Por tôpico Cesar Rabak via R-br
Maurício, Você não nos contou um dos aspectos que importam para estabelecer qual das duas funções seria mais apropriada: você está comparando três tratamentos ou dois tratamentos contra um controle? -- Cesar Rabak 2017-02-04 19:06 GMT-02:00 Maurício Sangiogo via R-br < r-br@listas.c3sl.ufpr.br>

Re: [R-br] boxplot com outliers, como colocar os valores?

2017-02-06 Por tôpico Cesar Rabak via R-br
Cristina, tentei rodar seu CMR, mas ele falhou com: Error in rect(x$breaks[-nB], 0, x$breaks[-1L], y, col = col, border = border, : objeto 'cores' não encontrado Quanto a plotar os valores que você deseja, recomenda-se, em geral, evitar ao máximo esse expediente. Se o gráfico não apresenta as

Re: [R-br] Trabalhando com datas

2017-02-06 Por tôpico Cesar Rabak via R-br
Edson, Não sabendo exatamente para quê você necessita a diferença, talvez lhe convenha ter este aditamento: a diferença após a conversão das cadeias de caracteres com strptime() pode ser feita diretamente: > hht[2]- hht[1] Time difference of 30 mins se você quiser apenas o valor numérico, pode-se

Re: [R-br] Modelo de Clutter - método dos mínimos quadrados em dois estágios

2017-02-10 Por tôpico Cesar Rabak via R-br
Gustavo, Veja no guia de postagem, http://www.leg.ufpr.br/r-br-guia como você deve apresentar essa pergunta para ter mais chances de obter uma orientação efetiva, em particular, forneça o CMR. Uma sugestão, enquanto você não atende a essa diretriz, é você rodar o exemplo das funções dentro do R c

Re: [R-br] Distribuição para regressão de resposta binária

2017-02-10 Por tôpico Cesar Rabak via R-br
Marcos, Dei uma olhada "perfunctória" nas respostas e me fixei nos resultados de duas regressões que você enviou, na mensagem de nove do corrente 10h34, e a primeira mensagem dia sete. O que me salta os olhos nas duas regressões é o valor-p do intercepto: Na regressão do dia sete: > Mod1 <- glm

Re: [R-br] Fatorial incompleto

2017-02-13 Por tôpico Cesar Rabak via R-br
Fernando, O conceito de "fatorial incompleto" vem do desenho do experimento não do acidente, como é o seu caso. Outrossim, você tem tão poucos casos para a IDADE==1, que comparar esses resultados com as outras idades pode ser pouco producente. Qual seria a variável de interesse na comparação e q

Re: [R-br] Fatorial incompleto

2017-02-14 Por tôpico Cesar Rabak via R-br
Assis via R-br < > r-br@listas.c3sl.ufpr.br> escreveu: > > Lembrando que planejamento de experimento é uma coisa planejada, você > desenha como vai ser, não espera acontecer para ver onde encaixa. > > O mesmo vale para amostragem > > Em seg, 13 de fev de

Re: [R-br] Referências para justificar o uso de teste de médias baseados em preditores lineares [OFF-TOPIC]

2017-02-28 Por tôpico Cesar Rabak via R-br
Prez@dos, Como a resposta do Walmes está completa, no meu [humilde] ponto de vista (e bate com minha experiência), gostaria de tentar aditar oferecendo uma abordagem complementar mais intuitiva, mas que pode abrir a famosa lata de minhocas (*can of worms*) como gostam de colocar os americanos: An

Re: [R-br] Simulação com Monte Carlo

2017-02-28 Por tôpico Cesar Rabak via R-br
Rodolfo, Hodiernamente o "método Monte Carlo" tem sido usado como um apelido para quase qualquer forma de calcular um resultado empregando geração de números pseudo aleatórios, portanto a resposta pode ser dependente do ambiente, instituição e/ou instrutor a qual você deve os resultados. O uso de

Re: [R-br] Teste de Normalidade Multivariada

2017-03-09 Por tôpico Cesar Rabak via R-br
Pollyanna, Veja se você realmente precisa fazer esses testes apesar de tudo! Eu já fiz vários posts criticando essa prática, incorreta, de testar normalidade, que geralmente serve como um "seguro", na verdade uma muleta psicológica, para efetuar análises que têm como "pressuposto" a normalidade d

Re: [R-br] Teste de Normalidade Multivariada

2017-03-12 Por tôpico Cesar Rabak via R-br
Marcelo Laia, "Extremamente robusto" não uma medida que se possa usar para determinar se algo é adequado ou não, no atacado. Numa situação onde seu orientador/patrão e você conseguem negociar essa métrica qualitativamente, ou talvez a inclusão dos "anonymous referees" ou editores daquela famosa pu

Re: [R-br] Dúvida de como usar o split com o Anova do pacote Car

2017-03-13 Por tôpico Cesar Rabak via R-br
Olympio, O que diz a documentação da função car::Anova()? 2017-03-13 15:17 GMT-03:00 Olympio Neto via R-br : > Prezados, > > Quero fazer um desdobramento interação em experimento fatorial muito > semelhante ao que foi apresentado na página R-íduculas (Desdobramento de > interação em experiment

Re: [R-br] grafico de distribuicao acumulada de uma variavel continua

2017-03-17 Por tôpico Cesar Rabak via R-br
http://bfy.tw/Ahi3 2017-03-16 14:49 GMT-03:00 Nilson Guiotoku via R-br < r-br@listas.c3sl.ufpr.br>: > Bom dia! > > > Eu preciso construir o gráfico da densidade acumulada de uma variável > continua (domínio dos Reais) que não tem um modelo identificado. Usei o > gráfico de pareto construído por t

Re: [R-br] Ajuda comando restruturação

2017-03-18 Por tôpico Cesar Rabak via R-br
Procure pela documentação do pacote *reshape* do R. HTH -- Cesar Rabak 2017-03-18 10:44 GMT-03:00 Edmar Caldas via R-br : > > > Bom Dia! > > Preciso de uma ajuda e não sei como começar. > > tenho uma arquivo que está na coluna e quero que ele fique na linha. segue > exemplo. Estou procurando um

Re: [R-br] Ajuda no calculo entre duas datas

2017-03-24 Por tôpico Cesar Rabak via R-br
Edmar, Se você quer fazer o cálculo "na unha" para o número de dias, o melhor é usar o divisor 365.25 para levar em conta os anos bissextos¹ Enquanto as propostas tanto de floor() como as.integer() façam o que você especificou, veja se não vale a pena você arrendondar levando em conta a fração, p

Re: [R-br] Selecionar linhas e colunas

2017-03-24 Por tôpico Cesar Rabak via R-br
Na verdade, é melhor: subset(raiox, Fratura=="Ausente *&&* id >=50) raiox[raiox$Fratura=="Ausente" *&&* raiox$id>=50,] 2017-03-24 10:49 GMT-03:00 Fernando Antonio de souza via R-br < r-br@listas.c3sl.ufpr.br>: > Acredito que > > subset(raiox, Fratura=="Ausente & id >=50) > > raiox[raiox$Fratura

Re: [R-br] Função aggregate.

2017-03-24 Por tôpico Cesar Rabak via R-br
o parâmetro FUN só pode aparecer uma vez em aggregate, se você precisa de uma função que "faça mais coisas" defina-a com mais parâmetros, como, p. ex. no s/caso, ...function(data_int, valor_int). 2017-03-24 13:36 GMT-03:00 Edmar Caldas via R-br : > > > Estou com um arquivo neste formato. > > [i

Re: [R-br] GLM

2017-03-29 Por tôpico Cesar Rabak via R-br
Ana Paula, Se você está interessada na teoria da extração, v. a obra de G. A. Paulo "Modelos de Regressão com apoio computacional", IME-USP, 2013, disponível no link da página dele: https://www.ime.usp.br/~giapaula/texto_2013.pdf, os exemplos do livro são na linguagem S-plus/R, portanto você terá

Re: [R-br] Calibratio Plot

2017-03-30 Por tôpico Cesar Rabak via R-br
Ary, As formas de fazer esse gráfico (especialmente se você abre mão da complexidade do histograma superposto) são muitas no R, veja plot() no pacote base (para colocar a linha de fé você pode acrescentá-la com lines(), xyplot() no pacote Lattice, scatterplot() do pacote car, enfim esta lista não

Re: [R-br] Data frame não reconhece uma variavel

2017-03-31 Por tôpico Cesar Rabak via R-br
Uma ou mais das idades gerou um NA, tente > mean(idade, na.rm = TRUE) _se_ funcionar, você deve investigar qual das idades gera o NA (provavelmente erro na data de nascimento) OBS.: como a data hoje é "muda" é uma "variável que não varia", talvez você devesse fazer: > Banco$idade <- difftime(Sy

Re: [R-br] Ajuda: If else

2017-04-07 Por tôpico Cesar Rabak via R-br
Edmar, No seu R, tecle o comando e leia o help sobre o assunto que você tem dúvida: >?Logic HTH -- Cesar Rabak 2017-04-03 17:07 GMT-03:00 Edmar Caldas via R-br : > Pessoal estou com outra dúvida, preciso criar uma coluna de > inconsistência que funciona assim. Zero quando a condição está ok

Re: [R-br] Duvida em estatística no R

2017-04-08 Por tôpico Cesar Rabak via R-br
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 .

Re: [R-br] Como escolher o melhor conjunto de dados imputador no MICE

2017-04-21 Por tôpico Cesar Rabak via R-br
Elias, A imputação é uma técnica que super simplificadamente objetiva colocar dados em variáveis cujos valores estão ausentes para não "perturbar" uma determinada estatística e permitir analisar as outras variáveis dos casos sob estudo. Um exemplo simples seria colocar a média da variável em todo

Re: [R-br] Como escolher o melhor conjunto de dados imputador no MICE

2017-04-24 Por tôpico Cesar Rabak via R-br
> > Prof. Edson Lira, Me > Estatístico > Manaus-Amazonas > > > Em Sexta-feira, 21 de Abril de 2017 10:16, Cesar Rabak via R-br < > r-br@listas.c3sl.ufpr.br> escreveu: > > > Elias, > > A imputação é uma técnica que super simplificadamente objetiva colocar >

Re: [R-br] SOMAR LINHAS

2017-05-03 Por tôpico Cesar Rabak via R-br
Veja a documentação para o comando agregate, especialmente a interface com uso de "fórmulas". Talvez você tenha que fazer a concatenação de cada soma. YMMV -- Cesar Rabak 2017-05-03 14:03 GMT-03:00 jose gomes de souza filho Gomes via R-br < r-br@listas.c3sl.ufpr.br>: > Boa Tarde, > > Tenho a s

Re: [R-br] SOMAR LINHAS

2017-05-04 Por tôpico Cesar Rabak via R-br
opie a saída e cole no email, para que possamos refazer > e lhe ajudar aqui. > > Att > > Em 3 de maio de 2017 22:36, Cesar Rabak via R-br > escreveu: > >> Veja a documentação para o comando agregate, especialmente a interface >> com uso de "fórmulas". >> >

Re: [R-br] Ajuda com Erro em DataFrame

2017-05-15 Por tôpico Cesar Rabak via R-br
Pesso@ll, Geralmente os objetos retornados pelas funções mais complexas do R retornam o *dataframe* empregado também na regressão. No caso em tela, o *dataframe* já expurgado dos casos *missing* está no elemento da lista "model": No caso hipotético do OP, ele poderia fazer o seguinte: novo_banc

Re: [R-br] knitr - definindo permissões qundo estou Gerando PDF

2017-05-15 Por tôpico Cesar Rabak via R-br
Como esta é meio off, dirijo-me apenas a você, Leonard: Você coloca que sabe fazer isso no TeX, você sabe como fazer isso no LaTeX também? Poderia me passar algum link a respeito? Gto! 2017-05-15 9:05 GMT-03:00 Leonard de Assis via R-br < r-br@listas.c3sl.ufpr.br>: > Bom dia a todos, > > > >

Re: [R-br] Fwd: Des-transformar variáveis

2017-05-19 Por tôpico Cesar Rabak via R-br
Hugo, O nome dessa operação que deseja fazer em inglês é « backtranformation ». Numa das n respostas a essa questão nos *fora* do R a proposta de se usar uma função: > invBoxCox <- function(x, lambda) if (lambda == 0) exp(x) else (lambda*x + 1)^(1/lambda) você poderá então retrotransformar s

Re: [R-br] R + COM - Component Object Model

2017-06-02 Por tôpico Cesar Rabak via R-br
Cleber, Esses pacotes são mantidos pelos seus idealizadores/desenvolvedores e IIRC há uma versão para .NET tão recente como 2016. . . Você já olhou lá? 2017-06-01 20:54 GMT-03:00 Cleber N.Borges via R-br < r-br@listas.c3sl.ufpr.br>: > Olá a todos, > Alguém conhece algum pacote ou modos alterna

Re: [R-br] Integral de um modelo não linear para calculo da área abaixo da curva

2017-06-14 Por tôpico Cesar Rabak via R-br
Alexandre, Considerando que os dados da sua curva "empírica" já estão dados, é muito mais interessante do ponto de vista matemático e prático (inclusive no R) fazer a integral por meios numéricos, especialmente por você deixar claro que deseja a integral definida (« área abaixo da curva para meus

Re: [R-br] Escolha do valor que mais ocorre

2017-06-15 Por tôpico Cesar Rabak via R-br
Elias, Começando pelo mais fácil: A resposta à sua questão número 3 é que não se define moda para dados dicotômicos porque ou é um ou outro nível que teria a maior frequência (noves fora o "empate"), dependendo do tipo de análise pode-se se falar em "prevalência", etc. Com relação às duas outras

Re: [R-br] Resgate de um valor de um summary

2017-06-16 Por tôpico Cesar Rabak via R-br
Sim, você precisa dar uma olhada no objeto cfaOut com uma chamada a str: > str(cfaOut) e ver em que lista ou dataframe embutido econtram-se os índices que lhe interessam. . . 2017-06-15 20:21 GMT-03:00 Clodoaldo José Figueredo via R-br < r-br@listas.c3sl.ufpr.br>: > Obrigado, mas ao executar

Re: [R-br] Vetor de Matrizes

2017-06-19 Por tôpico Cesar Rabak via R-br
Clodoaldo, Sua questão ou é estranhamente primária e simples ou há algo que não se percebe que necessitaria que você seguisse o guia de postagem deste grupo e em particular o CMR: 1) Para "armazenar essas matrizes" você as atribui a variáveis, no seu caso uma variável indexada serviria; 2) Para "

Re: [R-br] Assunto: Criar novo dataframe a partir de uma estrutura

2017-06-19 Por tôpico Cesar Rabak via R-br
Tente (*caveat emptor*, código não testado): > df2 <- dataframe.original[FALSE,] teste: > str(dataframe.original) e > str(df2) HTH -- Cesar Rabak 2017-06-19 16:46 GMT-03:00 Edson Lira via R-br : > Não seria o caso de subset? > > Enviado do Yahoo Mail no Android >

Re: [R-br] Baixar imagens

2017-06-21 Por tôpico Cesar Rabak via R-br
Sem CMR. . . a msg que você recebe parece ser falta de "fechar" o arquivo. On Wed, Jun 21, 2017 at 11:29 AM, Nicolli Albuquerque via R-br < r-br@listas.c3sl.ufpr.br> wrote: > Estou tentando baixar fotos através da função download.file: > > > foto<-as.character(linksvalidosespecies$link[50]) > s

Re: [R-br] Somas de Quadrados Tipo III

2017-06-21 Por tôpico Cesar Rabak via R-br
Eduardo, Além das excelentes respostas que vc já recebeu, note que até a *terminologia* sobre "Somas de Quadrados Tipos II e III" é na verdade estranha à Estatística e sim um artifício inventado por um dos 'fabricantes' de pacote de SW estatístico cuja terminologia "pegou". Comecei o post com ess

Re: [R-br] Ajuda com organização de dados

2017-06-25 Por tôpico Cesar Rabak via R-br
Jackson, Com a informação passada por você podemos apenas dar sugestões. Se você quiser uma resposta mais assertiva você precisa seguir o guia de postagem(http://www.leg.ufpr.br/r-br-guia), e em particular fornecer o CMR -- código mínimo reproduzível. Minha sugestão é que você veja a documentação

Re: [R-br] RES: Dúvida sobre pesos amostrais

2017-06-30 Por tôpico Cesar Rabak via R-br
Carlos, Essa dúvida não exatamente pertinente ao R, mas sim uma dúvida de Estatística Básica. No seu caso você obteve uma proporção da preferência por uma variedade de serviço. Não existe o conceito de "extrapolação" em amostragem, mas sim a questão do intervalo de confiança -- IC dessa estimati

Re: [R-br] Duvida estatística com R

2017-07-04 Por tôpico Cesar Rabak via R-br
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 indi

Re: [R-br] GLM multivariado

2017-07-12 Por tôpico Cesar Rabak via R-br
Ana Paula, Um texto bem completo que inclui o assunto de seu interesse está em https://www.ime.usp.br/~giapaula/textoregressao.htm 2017-07-11 9:21 GMT-03:00 ana paula coelho madeira via R-br < r-br@listas.c3sl.ufpr.br>: > Bom dia! > > > Alguém tem alguma referência de estudos com GLM multivaria

Re: [R-br] tratamento de pontos discrepantes (pesos de bovinos)

2017-07-20 Por tôpico Cesar Rabak via R-br
Prezada Adriele, A questão de "detecção", e mais ainda, a remoção de dados aberrantes em bancos de dados tem várias dimensões, dentre elas a puramente matemática. Como você menciona uma fonte desses dados aberrantes como erros de balança,imagino que os outros pontos seriam escrutinizados antes de

Re: [R-br] geom_boxplot

2017-07-21 Por tôpico Cesar Rabak via R-br
http://bfy.tw/CxnM HTH -- Cesar Rabak 2017-07-21 8:36 GMT-03:00 Fabiani via R-br : > Bom dia listeiros, estou com uma dúvida, e talvez alguém possa me indicar > alguma maneira de resolver. > > Estou com uma série histórica de dados diários de chuva, e preciso fazer > um boxplot dos valores agr

Re: [R-br] Erro de memória importando objeto pequeno de SQL

2017-07-28 Por tôpico Cesar Rabak via R-br
Pedro, Mesmo sem poder fazer uma tentativa, me parece que o erro é estranho demais para não ser notado: a chamada à primitiva da ling. C "calloc" está pedindo 200 Gbytes, que falha e me parece algo grande demais para um "default" da chamada a sqlQuery(). HTH 2017-07-28 9:42 GMT-03:00 Pedro Emma

Re: [R-br] Erro de memória importando objeto pequeno de SQL

2017-07-29 Por tôpico Cesar Rabak via R-br
Pedro, Calloc é uma chamada "primitiva" para o R, uma função da biblioteca C e é uma função feita para "alocar" (reservar) memória e zerá-la. O estranho, para mim, é a tentativa do aplicativo de reservar *tanta* memória assim na chamada da função do R. Essa função sqlQuery() pode ser testada faz

  1   2   >