Re: [R-br] Anova fatorial com medidas repetidas (e aninhamento?)

2018-10-14 Por tôpico Cesar Rabak por (R-br)
Thallita, O excerto dos seus dados mostra que (pelo menos para o R) o fator sensibilidade tem dois níveis: *sensível* e *tolerante*. Escrevendo de outra forma um caso ou é sensível ou é tolerante. Então: o quê você quiz dizer com a expressão "subconjunto de"? On Thu, Oct 11, 2018 at 10:22 PM

Re: [R-br] Anova fatorial com medidas repetidas (e aninhamento?)

2018-10-15 Por tôpico Cesar Rabak por (R-br)
; > > > - > *Thallita Grande, M.Sc.* > Doutorado em Ecologia > Instituto de Ciencias Biologicas > Universidade de Brasilia - Campus Darcy Ribeiro > Skype: thallita.grande > CRBio 080256/04-D > -- > *De:* R-br em nome

Re: [R-br] Análise de experimentos em faixas (medidas repetidas no tempo)

2018-10-30 Por tôpico Cesar Rabak por (R-br)
a função dic do ExpDes ou a design.rcbd do agricolae não lhe servem? On Sun, Oct 28, 2018 at 8:13 AM Andre Oliveira por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > bom dia, > de todos os métodos análises que busquei em livros, no google e aqui na > lista r-br, sempre achei análises de experimento

Re: [R-br] GLM model

2018-11-09 Por tôpico Cesar Rabak por (R-br)
A resposta está no 1º parágrafo da *vignette* para o pacote quasi: « Computing “quasi-AIC” (QAIC), in R is a minor pain, because the R Core team (or at least the ones who wrote glm, glmmPQL, etc.) are purists and don’t believe that quasi- models should report a likelihood. » HTH -- Cesar Rabak

Re: [R-br] usar funções em C dentro do objeto DLL

2018-11-19 Por tôpico Cesar Rabak por (R-br)
Se a lógica de chamar uma função C se mantém com a API do R que vc está usando, necessitar-se-ia que na chamada a FDwfGetVersion você passasse uma posição de memória para a função preencher a cadeia (a.k.a. *ponteiro* para memória no jargão C). O que a documentação da biblioteca R para ligação co

Re: [R-br] usar funções em C dentro do objeto DLL

2018-11-26 Por tôpico Cesar Rabak por (R-br)
Cleber, Seu raciocínio no C está correto com o *proviso* que cadeia szVersion tem que estar pré-alocada antes da chamada à função FDwfGetVersion da API do Windows. Veja na página 7 do manual do SDK que você está usando a documentação para a função em questão diz: FDwfGetVersion(char szVersio

Re: [R-br] usar funções em C dentro do objeto DLL

2018-11-26 Por tôpico Cesar Rabak por (R-br)
Concordo que possa ser uma saída. On Mon, Nov 26, 2018 at 5:08 PM Rodrigo Ângelo wrote: > Cleber, > > Talvez dê pra você fazer uma outra dll em C que encapsula essa e adiciona > as bibliotecas que você precisa. > > Att, > Rodrigo > > On Mon, Nov 26, 2018 at 4:54 PM Cleber Borges por (R-br) < > r

Re: [R-br] COMPARAR TRATAMENTOS

2018-12-02 Por tôpico Cesar Rabak por (R-br)
Edson, Com seis casos por tipo de tratamento, não vejo como você conseguirá somente com um teste estatístico separar dos tratamentos o "mais eficiente". Com dez tratamentos você teria que fazer 45 testes o que exigiria que o limiar de rejeição da hipótese nula (que todos os tratamentos são equiva

Re: [R-br] análise de correlação entre variáveis qualitativas e quantitativas

2018-12-07 Por tôpico Cesar Rabak por (R-br)
Prezada Flávia, A distribuição de quaisquer forma de renda ou seus *proxies* como chamam os americanos, em quase a totalidade das populações humanas onde esse conceito é aplicável não é gaussiana ("Normal") mas sim tende a seguir a curva descoberta por Pareto, portanto nada de inesperado aqui. Qu

Re: [R-br] IC para previsao

2018-12-08 Por tôpico Cesar Rabak por (R-br)
Edmar, Por gentileza: Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. Também conhecido por "CMR". On Fri, Dec 7, 2018 at 10:37 PM Edmar Caldas por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Preciso de uma ajuda para criar o IC para ajuste do mode

Re: [R-br] Digest R-br, volume 96, assunto 8

2018-12-10 Por tôpico Cesar Rabak por (R-br)
Edmar, O CMR precisa atender a cláusula *reproduzível *se você espera que a gente possa entender seu problema: R version 3.5.1 (2018-07-02) -- "Feather Spray" Copyright (C) 2018 The R Foundation for Statistical Computing Platform: i386-w64-mingw32/i386 (32-bit) R é um software livre e vem sem GA

Re: [R-br] Fwd: Criação de um dataframe a partir das diferenças

2018-12-13 Por tôpico Cesar Rabak por (R-br)
Silvia, O suposto resultado não apresenta *um* sequer dos dados apresentados como exemplo 🤔. . . Daria para você mostrar como ficaria o resultado se os dados do exemplo fosse processado como você gostaria? On Thu, Dec 13, 2018 at 2:53 PM SILVIA PATRICIA BARRETO SANTANA por (R-br) < r-br@listas.

Re: [R-br] Distribuição gama bivariada no R

2018-12-20 Por tôpico Cesar Rabak por (R-br)
Faça uma pesquisa no site CRAN por *Task View* "gamma bivariate" On Thu, Dec 20, 2018 at 4:20 PM Ana Beatriz Souza Souza por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Boa tarde pessoal! > > Estou trabalhando com distribuições gama bivariada e outras distribuições > relacionadas. Alguma indicaçã

Re: [R-br] Particionar soma de quadrado

2019-01-09 Por tôpico Cesar Rabak por (R-br)
Andre, O seu CMR está cheio de problemas... A linha: fit=aov(indice ~ cultivar*factor(dose)+Error(bloco:(cultivar+dose)) dá erro por falta de um parêntese no final. Colocando-se o símbolo faltante caímos em: > fit=aov(indice ~ cultivar*factor(dose)+Error(bloco:(cultivar+dose))) Warning messag

Re: [R-br] Particionar soma de quadrado

2019-01-10 Por tôpico Cesar Rabak por (R-br)
Acho que não. . . a página indicada faz esta contribuição acaciana que IMNSHO não ajuda um iota OP na solução da ANOVA dele: « . . . This can happen due to a bug in your programming, a participant being noncompliant, data trimming after the fact, or a whole host of other reasons. . . . » O caso do

Re: [R-br] Particionar soma de quadrado

2019-01-10 Por tôpico Cesar Rabak por (R-br)
André, o fato de a soma dos quadrados estar correta não muda o fato de você estar fazendo o equivalente a encontrar uma reta com apenas dois pontos e querer tratar isso como uma generalização e estudar o valor-p da ANOVA, etc. A razão porque aquela função é denominada *Error* é que ela deveria aju

Re: [R-br] Modelo em aov semelhante ao do ExpDes.pt

2019-01-10 Por tôpico Cesar Rabak por (R-br)
Usando o código do próprio pacote: > fat3.dbc Não resolve seu problema? On Thu, Jan 10, 2019 at 11:55 AM Marcelo Laia por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Colegas, > > Minha aluna está rodando no ExpDes o seguinte modelo: > > attach(clorofila) > fat3.dbc(Isolado,Genótipo,Tempo,Bloc

Re: [R-br] Linha de tendência em xyplot (lattice): melhor ajuste

2019-01-11 Por tôpico Cesar Rabak por (R-br)
Marcelo, Tanto o LOESS como a Spline são curvas "heurísticas" portanto, embora nada seja impossível (que alguém tenha inventado algum tipo de *pseudo R²* para elas), não teriam muito sentido porque a parte que o R² supostamente é chamada a responder é o ajuste a um *modelo* dos seus dados que prec

Re: [R-br] Ajuda com a passagem incorreta de variável string para uma função do R

2019-01-11 Por tôpico Cesar Rabak por (R-br)
Bruno, apenas uma dúvida: como você espera fazer um *merge* em dois conjuntos de dados cujas "chaves" (no data.table, que um dataframe) as variáveis ou colunas... ou sendo mais simplório: como você espera que o SW estatístico *saiba* que Bolsonaro na data.table perc_validos é a *mesma coisa* que

Re: [R-br] Ajuda com a passagem incorreta de variável string para uma função do R

2019-01-15 Por tôpico Cesar Rabak por (R-br)
tratarei na sequência do meu código. Para mim o importante é > acertar a seleção para o cruzamento... > > Assim, como eu faço para selecionar uma das colunas com contagem de votos > além da variável chave nas funções? O que eu estou fazendo de errado? > > Grato mais uma vez pela at

Re: [R-br] Ajuda com a passagem incorreta de variável string para uma função do R

2019-01-18 Por tôpico Cesar Rabak por (R-br)
t;. >>> >>> Meu desafio é separar, dentre as várias colunas na tabela *votos_regiao*, >>> apenas a coluna chave (*Regiao*) e a coluna que eu apontar como >>> parametro na função (neste exemplo acima a "Demais Candidatos"). Da maneira >>> que

Re: [R-br] JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES

2019-01-29 Por tôpico Cesar Rabak por (R-br)
Edson, Veja uma trédi recente nesta lista sobre o uso da função merge ela parece servir exatamente para o que você deseja. Na trédi você poderá ver que podes escolher entre a merge da biblioteca base ou a do pacote data.table. HTH -- Cesar Rabak On Fri, Jan 25, 2019 at 6:40 PM Edson Lira por (

Re: [R-br] JUNTAR DOIS ARQUIVOS COM DIMENSÕES DIFERENTES

2019-01-30 Por tôpico Cesar Rabak por (R-br)
wrote: > Cesar, vc tem o link? Não consegui encontrar. > > > [ ]'s > Prof. Edson Lira, Me > Estatístico > Manaus-Amazonas > > > Em terça-feira, 29 de janeiro de 2019 10:21:39 AMT, Cesar Rabak por (R-br) > escreveu: > > > Edson, > > Veja

Re: [R-br] "probability heat map" no R

2019-02-14 Por tôpico Cesar Rabak por (R-br)
E como essa probabilidade seria calculada a partir de "friction"? Sendo que "distance" e "position" são dados como você que identificar seus resultados no *heat map*? On Thu, Feb 14, 2019 at 11:37 AM Luiz Leal por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Prezados, bom dia. Tenho o seguinte c

Re: [R-br] "probability heat map" no R

2019-02-18 Por tôpico Cesar Rabak por (R-br)
iz > > > [image: Inline image] > > > On Thursday, February 14, 2019, 2:54:10 PM GMT-2, Cesar Rabak por (R-br) < > r-br@listas.c3sl.ufpr.br> wrote: > > > E como essa probabilidade seria calculada a partir de "friction"? > > Sendo que "distance"

Re: [R-br] Heatmap plot

2019-02-21 Por tôpico Cesar Rabak por (R-br)
Sua pergunta está algo difusa ou vaga, mas admitindo como premissa que você já tenha a matriz de correlações, talvez obtida de outro cálculo ou oriunda de outro SW, a resposta seria sim salvo pelo detalhe que a matriz precisa ser representada no formato esperado para as funções do R que plotarão se

Re: [R-br] Teste de Qui-quadrado aplicado para a distribuição binomial

2019-03-06 Por tôpico Cesar Rabak por (R-br)
Na minha tentativa de rodar s/exemplo antes da linha > res As seguintes msgs de alerta são impressas: Warning messages: 1: In densfun(x, parm[1], parm[2], ...) : NaNs produzidos 2: In chisq.test(table(emp), table(est), correct = TRUE) : Aproximação do qui-quadrado pode estar incorreta 3: In dens

Re: [R-br] Teste de Qui-quadrado aplicado para a distribuição binomial

2019-03-06 Por tôpico Cesar Rabak por (R-br)
Alexandre: Por favor explique-no o seu raciocínio para fazer este teste aqui: *> result<-ifelse(testChi$p.**value>0.05,"Binomial Negative","Other distribution")* ## Classifica os testes em binomial negativa ou outra distribuição HTH -- Cesar Rabak On Tue, Mar 5, 2019 at 5:07 PM ASANTOS por (R-

Re: [R-br] two-way anova - violação das pressuposições

2019-03-28 Por tôpico Cesar Rabak por (R-br)
Embora chegando atrasado a esta discussão, gostaria de contribuir com algumas informações. Primeiramente, uma questão de conceito: há uma certa confusão na literatura, especialmente quando ela é ligada a algum manual de SW recente a respeito da "premissa de normalidade" numa ANOVA. Não é distribu

Re: [R-br] Redistribuição das contagens de uma estrutura de dados

2019-04-01 Por tôpico Cesar Rabak por (R-br)
Apenas aditando a elegante solução proposta pelo Alan, Este comando: > grep(",",colnames(df), invert=T) permitiria identificar todas as colunas com nomes "simples" no jargão do OP. Então, um primeiro data.frame poderia ser construído : df2 <- df[grep(",",colnames(df), invert=T)] Aí um laço f

Re: [R-br] rfm - retornando "0" linhas

2019-04-05 Por tôpico Cesar Rabak por (R-br)
View com "V" maiúsculo 🤔? On Fri, Apr 5, 2019 at 4:12 PM Edmar Caldas por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Boa Tarde! > Estou tentando criar um data.frame, uma View ou até mesmo salvar no > formato csv o resultado da RFV criada a partir do pacote RFM, mas ele > retorna com "0" linhas n

Re: [R-br] Redistribuição das contagens de uma estrutura de dados (Cesar Rabak / Alan Rodrigo Panosso)

2019-04-05 Por tôpico Cesar Rabak por (R-br)
👍 On Tue, Apr 2, 2019 at 2:09 PM Paulo Eduardo de Mesquita por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Boa tarde César e Alan > > Muito obrigado pela ajuda. > > Sobre a dúvida do César: > > " *. . . porque mesmo após as várias explicações, ainda não entendi se as > somas (como Alan fez) seria

Re: [R-br] Ler csv e enviar whatsapp

2019-04-20 Por tôpico Cesar Rabak por (R-br)
Como a pergunta é muito inespecífica, cabe aqui usar o método socrático, de forma a levantar o quê sabemos e o que não sabemos. *De plano* o qualificador "...usando o R." permite descartar a óbvia, mas provavelmente não desejada verdade que lendo-se no R os números de telefone pode-se enviar as me

Re: [R-br] Duvida sobre tipos de categoria e análise estatística aplicável

2019-05-17 Por tôpico Cesar Rabak por (R-br)
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. F

Re: [R-br] Duvida sobre tipos de categoria e análise estatística aplicável 2

2019-05-24 Por tôpico Cesar Rabak por (R-br)
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 quan

Re: [R-br] Converter dados

2019-05-24 Por tôpico Cesar Rabak por (R-br)
Pelo seu exemplo, você vai considerar um mês como trinta dias *sempre*. Correto? Esse erro não atrapalhará sua análise quando um período de tratamento inclua (ou não) o mês de fevereiro, por exemplo? On Fri, May 24, 2019 at 3:14 PM Emerson Cotta Bodevan por (R-br) < r-br@listas.c3sl.ufpr.br> wr

Re: [R-br] Como extrair ano de data

2019-05-24 Por tôpico Cesar Rabak por (R-br)
A resposta exata depende de como essa " 03/12/2013" está sendo vista na imagem do seu R. Se for uma cadeia de caracteres, veja a função gsub(), se for um tipo Date, user *format()* para obter apenas o ano. Frente às questões que você já postou, recomendaria a leitura da sua referência favorita o

Re: [R-br] [Loop em Análise de Tendência sobre série histórica]

2019-05-24 Por tôpico Cesar Rabak por (R-br)
Com relação à sua segunda questão, talvez valha a pena você considerar em vez de usar um *for loop*, o uso da função *apply()* com a margem escolhida adequadamente (pelo que entendi no seu caso seria as colunas). HTH -- Cesar Rabak PS.: Via de regra o uso de laços em vez do uso das capacidades ve

Re: [R-br] Aplicação de métodos de dados faltantes no ambiente R

2019-05-31 Por tôpico Cesar Rabak por (R-br)
Rapaz! Com 24 variáveis e apenas 40 objetos, que entendo aqui como número de dados, amostras ou "linhas" do banco de dados, você se arrisca a fazer um sobreajuste (*overfitting*)¹ e não uma análise estatística! Ademais, com dados faltantes em alguns casos você ainda vai uma nova questão do que es

Re: [R-br] Aplicação de métodos de dados faltantes no ambiente R

2019-05-31 Por tôpico Cesar Rabak por (R-br)
Caio, Neste parágrafo acho que é necessário esclarecer se você se expressou mal ou é problema de entendimento: - Modelo de regressão logística. A equação gerada permitiria estimar os > valores dos dados perdidos. Não assume uma distribuição normal das > variáveis de resposta. Ao contrário da regr

Re: [R-br] Método de estimação do glm

2019-06-04 Por tôpico Cesar Rabak por (R-br)
*Read The **Fine** Manual? **😎*? Especialmente a descrição da função no subitem *method*. HTH -- Cesar Rabak On Tue, Jun 4, 2019 at 4:11 PM Elias Carvalho por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Estou fazendo uma regressão logistica com o comando: > modelo.rl1 <- glm(V1 ~ V2, data = da

Re: [R-br] (sem assunto)

2019-06-12 Por tôpico Cesar Rabak por (R-br)
Não uma restrição teórica para o teste Kruskall Wallis que as pernas da amostra sejam balanceadas. Obviamente, o *bom senso* pediria que o desbalanço não seja tão grande a ponto de ter um grupo com uma ou duas amostras e outro na casa de milhares... HTH -- Cesar Rabak On Mon, Jun 10, 2019 at 9:

Re: [R-br] Experimento em faixas

2019-06-15 Por tôpico Cesar Rabak por (R-br)
Ana Paula, Essa *doença* de propor testes de normalidade para "testar" a premissa de normalidade e homoscedasticidade *a priori* de um teste estatístico é muito propalada por manuais, autores eviesados por fabricantes de pacotes de SW estatísticos, geralmente comerciais e caros, que para justifica

Re: [R-br] Problema para plotar intervalo de confiança no ggplot2 para modelo GLM de Poisson com junção de níveis

2019-06-19 Por tôpico Cesar Rabak por (R-br)
Já houve boas interações, mas seu exemplo me faz questionar: você escreveu ". . . uma variável resposta (Consumption) e duas explicativas, *sendo uma categórica e outra quali* . . . " se a variável Income é qualitativa, porque você deixa o programa tratá-la como numérica? HTH -- Cesar Rabak On

Re: [R-br] Boxplot no R

2019-06-21 Por tôpico Cesar Rabak por (R-br)
De fato o *help* da função boxplot() diz: > > range: this determines how far the plot whiskers extend out from the box. > If ‘range’ is positive, the whiskers extend to the most extreme data point > which is no more than ‘range’ times the interquartile range from the box. A > value of zero causes t

Re: [R-br] Contrastes - MLG

2019-08-07 Por tôpico Cesar Rabak por (R-br)
André, A função glht está usando a ordem "alfabética" dos seus fatores, como a função não deve ter meios de mudar isso (você já *olhou* no FM, não *olhou*?) 😉, a única forma seria fazer uma cópia do objeto retornado 'p' neste teu exemplo, e inverter os sinais de *todos* os contrastes e colocar os

Re: [R-br] Como reportar o valor do Odds Ratio

2019-08-28 Por tôpico Cesar Rabak por (R-br)
Elias, Você já recebeu respostas em relação ao aspecto *matemático* da questão, e portanto vou pular essa e partir para observações sobre a segunda parte da sua questão: O que reportar na tabela de resultados? Há diferentes respostas "corretas" para essa questão em função do "público alvo" e/ou

Re: [R-br] Problema criar colunas: socorro!

2019-09-01 Por tôpico Cesar Rabak por (R-br)
Ellerson, SE a amostra de dados que você passou é um bom exemplo, você tem um problema *anterior* para resolver que é o de multiplicidade de entradas com grafia levemente diferentes. . . Veja as linhas cinco e seis do seu exemplo que o enquadramento, embora igual do ponto de vista legal, está esc

Re: [R-br] Problema criar colunas: socorro!

2019-09-02 Por tôpico Cesar Rabak por (R-br)
u BD pela > cadeia de caracteres que contenha o "ART. 14 - CP" e a partir dessa lista > auxiliar ir eliminando as linhas e incluir numa nova coluna nos número de > processo correspondentes” > > Estou sem saber que passo a passo considerar aqui... > > > > > Em

Re: [R-br] Ajuda: Escolha de distribuição de probabilidade (GLM)

2019-09-30 Por tôpico Cesar Rabak por (R-br)
Chiara, A escolha da "melhor" distribuição deve ser sempre feita considerando o *processo* de geração dos seus dados e menos usando informações obtidas dos próprios dados « ... dados não são normais e também não apresentam homoscedasticidade, mesmo após realizar transformação (log) ... ». Outros

Re: [R-br] Organizar data.frame e sumarizar dados em transposição

2019-10-09 Por tôpico Cesar Rabak por (R-br)
Sem um CMR é difícil ser assertivo, mas o delineamento do que você deseja seria obter uma lista sem repetições da coluna dos CPF (*uniq*), e usando os itens dessa lista como parâmetro para *which* ou *subset* fazer a soma e contar o número de parcelas adicionando no novo dataframe esses resultados.

Re: [R-br] Organizar data.frame e sumarizar dados em transposição

2019-10-10 Por tôpico Cesar Rabak por (R-br)
t; > Obrigado > > 1. http://www.rpubs.com/Kassio_Ferreira/tutorialdplyr > > Marcelo > > > On 09/10/19 at 04:32, Cesar Rabak por (R-br) wrote: > >Sem um CMR é difícil ser assertivo, mas o delineamento do que você > >deseja seria obter uma lista sem repe

Re: [R-br] Correlação de Pearson com escala Likert

2019-10-17 Por tôpico Cesar Rabak por (R-br)
A análise estatística de resultados obtidos por meio de uma escala Likert tende a ser um assunto controverso, pois há bastante "espaço" para violações de quase quaisquer premissas das análises e especialmente os testes estatísticos. Com isso posto, e [minha] recomendação que você deve pesquisar ma

Re: [R-br] Como identificar colinearidade?

2019-10-17 Por tôpico Cesar Rabak por (R-br)
Chiara, Colinearidade é quando você pode expressar alguma das variáveis que você pensava que eram independentes (entre si) como combinação linear de outras duas (ou mais...). O gráfico anexo mostra uma porção de coisas que chamam a atenção: o número de pontos para TP parece ser único para o caso "

Re: [R-br] derivada da densidade de probabilidade

2019-11-02 Por tôpico Cesar Rabak por (R-br)
Nilson, O que você deseja fazer tem um nome que nossos antepassados que já enfrentaram esse problema deram: "construção da ROC pelo método bi-normal". Por transitividade matemática, você pode descobrir a curva logística que identifica o valor da sua variável que discrimina melhor as duas condiçõe

Re: [R-br] Design of Experiments: replicatas no ponto central e intervalo de confiança

2019-11-03 Por tôpico Cesar Rabak por (R-br)
Esse assunto está menos para tutorial e mais para livro texto. Qual obra de referência você está usando para consultar sobre Planejamento de Experimentos? On Sat, Nov 2, 2019 at 10:33 PM Nei por (R-br) wrote: > Prezados e Prezadas, > > Estou estudando Planejamento Fatorial Completo de 2 níveis

Re: [R-br] Dúvida de 'warning message'

2019-11-06 Por tôpico Cesar Rabak por (R-br)
Você entende o mecanismo de criação de "postos" (*ranks* em inglês)? Quando dois ou mais deles têm a mesma posição devido ao cálculo como é feito, você tem o *empate*. Caso você precise aprofundar, um livro texto de Estatística tem o material para esclarecer esse assunto. HTH -- Cesar Rabak On

Re: [R-br] densidade de uma variável

2019-11-11 Por tôpico Cesar Rabak por (R-br)
Nilson, Qual é a razão de você ter expectativa de uma variável normal ter intervalo [0,1] e não 0±d.p.? On Mon, Nov 11, 2019 at 6:33 PM Nilson Guiotoku por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Pessoal, > > Estou com uma dúvida sobre a aplicação da função density, em um banco de > dados d

Re: [R-br] Última Versão do R em ambiente Linux

2019-11-11 Por tôpico Cesar Rabak por (R-br)
Malgrado meu Linux seja Ubuntu, as instruções acima citadas valem para ambos, e posso lhe assegurar que *funcionam*: csrabak@salsicha:~$ R --version R version 3.6.1 (2019-07-05) -- "Action of the Toes" Copyright (C) 2019 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-

Re: [R-br] Tamanho do efeito para Wilcoxon

2019-11-14 Por tôpico Cesar Rabak por (R-br)
Elias, Tua pergunta enseja a oportunidade de discutir duas coisas: 1. O quê é efeito, e como pode ser medido; e 2. a confusão das pessoas sobre esse assunto devido à, desastrosa IMNSHO, introdução do conceito de [*tamanho de] efeito padronizado*, que um dos próprios autores da ideia

Re: [R-br] Desdobramento de um fatorial triplo

2019-11-19 Por tôpico Cesar Rabak por (R-br)
Prez@dos, Quero dar um pitaco nessa conversa sobre CV e *pressupostos* (que são um eufemismo para *normalidade*): Considerando que o importante são os resíduos, que quando são hígidos têm média zero, pela definição de CV seria indeterminado (tendendo ao infinito)! Por transitividade, se a média

Re: [R-br] Validação cruzada

2019-12-03 Por tôpico Cesar Rabak por (R-br)
Gente 🤔. . . Os competentes *scripts* apresentados *não* me parecem validação cruzada. . . apenas uma validação de um modelo com um conjunto seperado de dados dos de treino. Na validação cruzada deve-se fazer a quebra em dados de treino e teste em n "dobras" (*folds*) e os resultados das n valida

Re: [R-br] Validação cruzada

2019-12-05 Por tôpico Cesar Rabak por (R-br)
Ok ! On Wed, Dec 4, 2019 at 11:27 AM Walmes Zeviani wrote: > Cesar, > > De fato não é uma validação cruzada, mas um *hold out*. > Mas entendi que aquilo é apenas um ciclo do loop como CMR. > Já que apontou, dado o pouco que sei, o mesmo modelo poderia ser ajustado > com a lme4::lmer() ou a nlme:

Re: [R-br] Comparar elementos em 2 dataframes

2019-12-06 Por tôpico Cesar Rabak por (R-br)
Sérgio, Como você já recebeu a informação mais relevante (vc falhou em enviar um CMR [veja o rodapé automaticamente inserido nas msg deste grupo com instruções de como fazê-lo]), gostaria de tentar ajudar de forma paralela: A descrição do seu problema parece-me centrada num domínio de problema pa

Re: [R-br] JUNÇÃO DE GRÁFICOS DE LINHA

2020-01-19 Por tôpico Cesar Rabak por (R-br)
Uma forma "mais automática" de fazer usando o pacote car, supondo seus dados no dataframe dat. > library(car) > sp(quantidade ~ taxadecoleta, dat, groups=dat$especie) HTH -- Cesar Rabak On Sun, Jan 19, 2020 at 4:22 PM Chiara Lubich por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Oi, gente > Est

Re: [R-br] Criar numeração única para dados iguais

2020-02-03 Por tôpico Cesar Rabak por (R-br)
Admitindo que você gerou a classificação na variável grupo por meios adequados, ou seja *você já disponha *dessa coluna, minha sugestão é que você crie a variável GC pela conversão em inteiro do grupo : Código não testado: seu_dataframe$GC <- as.integer(seu_dataframe$grupo) Dependendo como a fo

Re: [R-br] knitr::kable() para pdf como aumentar o tamanho da fonte do caption?

2020-02-08 Por tôpico Cesar Rabak por (R-br)
*Qual* dos modelos você usou? Veja o modelo do artigo (eventualmente no *site* dele) e observe onde o *caption* é estabelecido. SE vc vai fazer um número pequeno de vezes isso, então outra alternativa seria pegar o .TeX intermediário e compilar com outra classe de documento. HTH On Fri, Feb 7,

Re: [R-br] Identificar valores fora da elipse

2020-02-10 Por tôpico Cesar Rabak por (R-br)
Veja a *vignette* do pacote SIBER. On Sat, Feb 8, 2020 at 9:58 PM Luiz Leal por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > > Prezados, como faço para criar um script que identifique que os > laboratórios 7042 e 9541 estão fora da elipse. > Os dados seguem em anexo. > Desde já agradeço. > Luiz >

Re: [R-br] ajuda com o R

2020-02-14 Por tôpico Cesar Rabak por (R-br)
Os seus dados geram um " NA/NaN/Inf " quando a chamada a qr( D %*% ss[[nt]] %*% D, tol = tol) é executada. . . Veja se vc consegue ver o que há de errado com eles. On Fri, Feb 14, 2020 at 4:26 PM Dayana Alves da Silva por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Estou tentando rodar uma anal

Re: [R-br] Editar dados em arquivo *txt em um data frame

2020-03-25 Por tôpico Cesar Rabak por (R-br)
O arquivo que o OP postou o *link* não lhe parece adequado? On Wed, Mar 25, 2020 at 6:33 PM Mauro Sznelwar por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Tem um arquivo para rodar e ver como fica? > > > Olá, > > Segue uma sugestão: > > > library(tidyverse) > myfile <- read_lines("log.txt") > >

Re: [R-br] Editar dados em arquivo *txt em um data frame

2020-03-25 Por tôpico Cesar Rabak por (R-br)
Veja esta linha no texto do OP: > #Arquivo original > https://www.dropbox.com/s/pncmjwl3camap6d/log.txt?dl=0 > Fazendo a entrada conforme a sugestão do Max: > str(myfile) chr [1:1445] "obj" "Learning Rate: 0.001, Momentum: 0.9, Decay: 0.0005" ... HTH -- Cesar Rabak On Wed, Mar 25, 2020 at 9:2

Re: [R-br] Editar dados em arquivo *txt em um data frame

2020-03-25 Por tôpico Cesar Rabak por (R-br)
Só gostaria de comentar, que o resultado mydf finda sendo um *tibble* e não mais "apenas" um *data.frame*. Se isso é vantajoso ou requer que o tiblle seja demovido para data.frame só Alexandre pode aquilatar... HTH -- Cesar Rabak On Wed, Mar 25, 2020 at 1:28 PM Max Richard Verginio por (R-br) <

Re: [R-br] Análise de dados microbiológicos - Estatística não paramétrica

2020-03-27 Por tôpico Cesar Rabak por (R-br)
Olá Ana Paula, O problema com seus dados tem um nome procure nos pacotes do R por *interval-censored data*. Os dados ficariam assim: (0, 1000), ... , 2780, ... (certo valor, determinado valor), etc. HTH On Fri, Mar 27, 2020 at 2:30 PM ana paula coelho madeira por (R-br) < r-br@listas.c3sl.ufpr.br

Re: [R-br] ggplot2 (Mapa/UDHs/IVS)

2020-03-29 Por tôpico Cesar Rabak por (R-br)
Marcelo, Seu código não atende à diretriz deste grupo em Guia de Postagem no R-br , em particular, a solicitação de fornecer o CMR (*código mínimo reproduzível*). As chances de você obter respostas são maiores quando essa recomendação é seguida. HTH On Sun, Mar

Re: [R-br] Regressão Logística pode ser influenciada pelo balanceamento de resultados de variável dependente

2020-04-21 Por tôpico Cesar Rabak por (R-br)
Elias, Como cada VI « apresenta uma prevalência pequena da doença », que eu entendo como "cada VI tem uma contingência cruzada pequena para o caso de doença e exposição" você realmente está com um conjunto de variáveis que não vão conseguir explicar a VD como os testes que foram descritos mostram

Re: [R-br] Duvidas Análise de Componentes Principais

2020-04-22 Por tôpico Cesar Rabak por (R-br)
Fernando, Há dois aspectos a ser vistos na sua ACP. O número de variáveis que você tem antes da ACP e número de casos para você fazer uma análise multivariada (nos dois casos). Um gráfico muito interessante para auxiliá-lo a decidir sobre seus resultados é o "mapa de indivíduos" que plota cada c

Re: [R-br] Duvidas Análise de Componentes Principais

2020-04-23 Por tôpico Cesar Rabak por (R-br)
0DkFQBpLztwXks > > Aproveito para pedir dicas de boas referencias para estudo, principalmente > como tirar as conclusões do PCA para o modelo. > > Atenciosamente > > > > > > Em qua., 22 de abr. de 2020 às 22:03, Cesar Rabak por (R-br) < > r-br@listas.c3sl.uf

Re: [R-br] Duvidas Análise de Componentes Principais

2020-04-23 Por tôpico Cesar Rabak por (R-br)
o > 5 componentes: Eu explico 78% da variação > > Existe uma variação explicada mínima a ser considerada na escolha do > número de componentes? > > Obrigado! > > > > On Apr 23 2020, at 9:10 pm, Cesar Rabak por (R-br) < > r-br@listas.c3sl.ufpr.br> wrote: > &g

Re: [R-br] Problema com a instalação de pacotes no R!

2020-04-27 Por tôpico Cesar Rabak por (R-br)
André, Uma forma sutil de esse problema aparecer na instalação de pacotes no Windows é ele estar carregado pelo R na sessão que você está usando para fazer a instalação. Certifique-se que você começa de uma imagem sem carga de nenhum pacote que deseja atualizar ou reinstalar para evitar esse tipo

Re: [R-br] IC em tabelas de frequência

2020-05-06 Por tôpico Cesar Rabak por (R-br)
1) Você deve contatar os autores do artigo; 2) Análise de Dados Categorizados, Paulino & Singer, Ed. Edgar Blücher, 2006; 3) Sim; há diversos: https://pt.lmgtfy.com/?q=confidence+interval+for+proportion+in+r&pp=1 On Wed, May 6, 2020 at 8:21 PM Emerson Cotta Bodevan por (R-br) < r-br@listas.c3sl.u

Re: [R-br] Duvida teste de correlação

2020-05-27 Por tôpico Cesar Rabak por (R-br)
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 ser

Re: [R-br] Variancia

2020-06-03 Por tôpico Cesar Rabak por (R-br)
Elias, O "enunciado" da sua questão traz vários pontos que precisariam ser melhor elucidados: Com um desbalanço na resposta como o indicado 42 variáveis independentes pode levar a análise (que entimemicamente aponta para a Regressão Logística ou uma da família log-linear) ao que se denomina sobre

Re: [R-br] Ratio

2020-06-06 Por tôpico Cesar Rabak por (R-br)
Eugenio, Como outros já apontaram, e o FAQ deste grupo também solicita, para se poder ajudar precisa-se, em geral, de um CMR (código mínimo reproduzível). Ademais, você faz uma pergunta num grupo de apoio a um SW estatístico usando o jargão da sua área (*Ratio* de ações) que reduz ainda mais o nú

Re: [R-br] Manipulando banco de dados no R

2020-06-06 Por tôpico Cesar Rabak por (R-br)
Wenceslau, SE um BD com as informações já existe, por que você não o consulta diretamente, com, por exemplo, subset() etc..? Do contrário, você vai ter que fazer *loops* e calcular as três porcentagens, usando-as para gerar o comando subset, e no caso de o resultado ter mais de u'a linha, tratá-l

Re: [R-br] help Pacote rsm, superfície de resposta

2020-06-08 Por tôpico Cesar Rabak por (R-br)
Você precisa estudar os manuais do R (ñ o RStudio) para entender a assim denominada "*formula interface*" que é praticamente uma sublinguagem dentro do R. De fato: > modelo<-lm(y~x1*x2) > modelo1<-lm(y~x1*x2) Produzem regressões idênticas. Como seria idêntico: > modelo2<-lm(y~x1+x2+x1:x2) Talvez v

Re: [R-br] GLM com alta proporção de zeros

2020-07-23 Por tôpico Cesar Rabak por (R-br)
Emerson: Algum dos pares de VIs do seu problema apresenta correlação muito alta (pelos nomes vejo dois pares altamente candidatos)? Qtos casos você para cada desfecho para calcular a RP? -- Cesar Rabak On Wed, Jul 22, 2020 at 5:03 PM Emerson Cotta Bodevan por (R-br) < r-br@listas.c3sl.ufpr.br>

Re: [R-br] Explicação gráfico e resultados da PCA

2020-07-24 Por tôpico Cesar Rabak por (R-br)
Chiara, Já fiz uma sugestão para Ana Paula em nove de março de 2018 sobre esse mesmo assunto nesta lista, repito aqui para evitar a pesquisa: Veja estas obras: http://math.agrocampus-ouest.fr/infoglueDeliverLive/enseignement/support2cours/livres/ana.do.R Livre - Analyse de données avec R - math.a

Re: [R-br] Diferença de todos os valores

2020-08-05 Por tôpico Cesar Rabak por (R-br)
Diogo, Na comunidade *R* temos um ditado que diz mais ou menos o seguinte (tradução livre deste escrivinhador) : « No *R* a pergunta nunca é *se* mas sim como *como* algo pode ser feito.» Pelo teor da sua pergunta você não domina os rudimentos da linguagem R (a minha evidência é o jeito como escr

Re: [R-br] Diferença de todos os valores

2020-08-05 Por tôpico Cesar Rabak por (R-br)
OK Jerônimo... Para mostrar como o ditado dos erristas é apropriado veja esta solução¹: > res<-NULL > vetor<-c(19,23,27,29) > l <- length(vetor) > for (i in 1:l) {ifelse(length(vetor)>1,res<-c(res,abs(car(vetor)-cdr(vetor))),print(res)); vetor <- cdr(vetor)} [1] 4 8 10 4 6 2 Numa implemen

Re: [R-br] Diferença de todos os valores

2020-08-06 Por tôpico Cesar Rabak por (R-br)
Você leu a nota de rodapé? On Wed, Aug 5, 2020 at 8:59 PM Mauro Sznelwar por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Não consegui rodar esta rotina: > > res<-NULL > > vetor<-c(19,23,27,29) > > l <- length(vetor) > > for (i in 1:l) > {ifelse(length(vetor)>1,res<-c(res,abs(car(vetor)-cdr(vet

Re: [R-br] Assunto: Re: Diferença de todos os valores

2020-08-06 Por tôpico Cesar Rabak por (R-br)
👍 😎 On Thu, Aug 6, 2020 at 12:03 AM Diogo Jerônimo por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > Mauro, acho que tem de baixar a package lisp, como indicou o César > > Enviado do Yahoo Mail no Android >

Re: [R-br] Ponto de corte Cook's Distance

2020-08-07 Por tôpico Cesar Rabak por (R-br)
Esses diagnósticos "sintéticos" sobre um objeto matemático são *avisos* para o pesquisador, mas não substituem a compreensão do assunto no domínio do problema. Usando o MIto da Caverna como ideia, sua regressão seria a *sombra* e o problema real não totalmente visível para você. A distância de Co

Re: [R-br] Artigo sobre p.value

2020-08-07 Por tôpico Cesar Rabak por (R-br)
Já procurou no servidor da lista? On Fri, Aug 7, 2020 at 9:38 AM Luiz Roberto Martins Pinto por (R-br) < r-br@listas.c3sl.ufpr.br> wrote: > > Bom dia. > > Há algum tempo alguém encaminhou para a lista um artigo sobre p.value. O > autor se referiu ao artigo como "o artigo que ele sempre quiz escre

Re: [R-br] Artigo sobre p.value

2020-08-07 Por tôpico Cesar Rabak por (R-br)
.cnpq.br/2732314327604831 > > "*The s**cience exists because there are patterns. * > * The patterns exist because God created them*. > * The statistic exists to research the patterns that God created.*" > > > > > Em sex., 7 de ago. de 2020 às 11:54, Cesar Rabak po

Re: [R-br] Mudar eixo numérico para labels

2020-08-07 Por tôpico Cesar Rabak por (R-br)
Diogo, Mude a chamada a plot para : > plot(c(min(dados$CR),max(dados$CR)),c(minesc,maxesc),type="n",main=Título,xlab="",ylab=""*, xaxt='n'*) Veja o acréscimo em negrito:___ *↥* e acrescente a seguinte chamada: > axis(1,at=1:5,labels=I

Re: [R-br] Artigo sobre p.value

2020-08-09 Por tôpico Cesar Rabak por (R-br)
Seguramente, e *mais* o uso do limiar de 0,05 é arbitrário, desde o tempo que o Fisher explicou porque, para os testes dele, que não usavam o conceito do Teste de Hipótese Nula "moderna" baseado no conceito de Erros Tipo I e Tipo II (que Fisher discordava). MAS, para usar a teoria com os dois tipo

Re: [R-br] Diferença entre o teste de normalidade normalidade sobre a variável resposta e sobre a aov() dos dados em função da variável explicativa

2020-08-22 Por tôpico Cesar Rabak por (R-br)
Pesso@ll, Aproveitando o excelente "gancho" do Rodrigo, e abusando da nossa língua pátria que permite aproveitar a ambiguidade dela para a gente usar numa discussão, gostaria de contribuir nesta *trédi* com os seguintes argumentos: 1) A pergunta foi « … Gostaria de saber a diferença entre realizar

Re: [R-br] Assimetria

2020-09-24 Por tôpico Cesar Rabak por (R-br)
Sem mais informações, a sugestão geral seria: plote usando escalas que não comprimam demais os dados na parte mais à direita do eixo das abscissas, ou se o domínio do problema permite usar uma transformação nessa escala. HTH -- Cesar Rabak On Tue, Sep 22, 2020 at 9:07 AM ana paula coelho madeira

Re: [R-br] Iniciar no Pyton

2020-10-22 Por tôpico Cesar Rabak por (R-br)
Caro Mauro, Para Pyton não precisa procurar mais: "iço non ekziste" 😀 ‼ Para Python a já madura linguagem de programação há tantos meios de discutir desde Google groups até os canais do Discord, que você vai ter que ver qual o meio que mais lhe apetece... Discutir vantagens de A × B invariavelme

Re: [R-br] Modelo Misto - Análise Experimento Switchback

2020-11-04 Por tôpico Cesar Rabak por (R-br)
Se a única variável manipulada ("independente") foram os tratamentos, e a resposta é unicamente uma única variável também, matematicamente você não acaba com a mesma análise, apenas por caminhos diferentes? On Wed, Nov 4, 2020 at 12:28 PM Fernando Souza por (R-br) < r-br@listas.c3sl.ufpr.br> wrot

  1   2   >