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

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

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 <-

Re: [R-br] Análise de amostra pequena

2018-07-17 Por tôpico Cesar Rabak via R-br
Elias, não consigo ver como você poderia usar qualquer técnica para obter mais informação a respeito do fenômeno em questão usando esses artifícios, como o que você citou e outro muito em voga o* bootstraping*, etc. Nesses casos, o melhor a fazer é tentar identificar apenas via a estatística

Re: [R-br] Modelos Mistos e Lattice (Grafico de predições)

2018-07-14 Por tôpico Cesar Rabak via R-br
> packages <- c("tidyr", "nlme""latticeExtra","lattice") Erro: unexpected string constant in "packages <- c("tidyr", "nlme""latticeExtra"" > > #Instala pacotes necessários, se ja tiver instalado ativa-o > for(p in packages){ + if(!require(p,character.only = TRUE)) install.packages(p) +

Re: [R-br] Análise de amostra pequena

2018-07-11 Por tôpico Cesar Rabak via R-br
Elias, Há vários possíveis comentários a fazer, embora, de antemão aviso-o, que pelo fato de você ter dezessete casos apenas, não há muito a melhorar. 1. Calcular os coeficientes de associação quando a hipótese de independência não pode ser descartada tende a ser um exercício equivalente

Re: [R-br] Dados DIC

2018-07-10 Por tôpico Cesar Rabak via R-br
Fúlvia, Sua pergunta extrapola dúvidas sobre o R e o pacote que você utilizou e vai direto para a Estatística, no seu caso Análise de Experimentos [DIC?]. Não nem foro nem do escopo deste grupo discutir esse tipo de questão, que lhe recomendo você consulte uma obra de referência sobre o assunto,

Re: [R-br] Mauro Sznelwar compartilhou The R Book.pdf com você

2018-06-17 Por tôpico Cesar Rabak via R-br
Mauro, Quando eu clico em "visualizar o arquivo" ele entra numa tela de login com o e-mail da lista, que obviamente não tenho senha... Se eu tento acessar via minha conta no Dropbox, ele reclama que estou tentando de um e-mail diferente (pois o registro do Dropbox é com outro e-mail)... Não dá

Re: [R-br] Criar um VETOR de um dataframe

2018-06-12 Por tôpico Cesar Rabak via R-br
Manasses, Já recebeste uma resposta com uma solução. Uma outra maneira seria você se familiarizar com a função *subset*. HTH -- Cesar Rabak 2018-06-12 0:36 GMT-03:00 Manasses Nóbrega via R-br < r-br@listas.c3sl.ufpr.br>: > Olá, pessoal! > > Este parece ser um problema fácil de resolver, mas

Re: [R-br] Erro na análise de regressão

2018-06-04 Por tôpico Cesar Rabak via R-br
OK, Angélica. Veja que seu plot indica uma varição não linear com relação à sobrevivência, ela é máxima para um certo valor de concentração do óleo, enquanto sua curva e uma regressão linear faz entender que o aumento do óleo diminui monotonicamente a taxa de sobrevivência. SE for feita uma

Re: [R-br] ggplot2_apresentar equação no gráfico

2018-06-04 Por tôpico Cesar Rabak via R-br
Henrique, Se você está pedindo para o SW, via geom_smooth(method = lm, formula = y ~ exp(-x), se = FALSE) fazer o alisamento de uma curva usando a equação y ~ exp(-x), como você esperaria que o SW "adivinhasse" que você quer calcular um coeficiente para o expoente? Por outro lado, se você está

Re: [R-br] Erro na análise de regressão

2018-06-04 Por tôpico Cesar Rabak via R-br
Angélica, O tratamento seria a variável *Oleo* e a resposta, a var. *S* ? Você já experimentou plotar esses dados? Qual é o seu pressuposto para chamar dic do pacote ExpDes.pt, dessa maneira? -- Cesar Rabak 2018-06-04 20:15 GMT-03:00 Angélica Ricarte via R-br < r-br@listas.c3sl.ufpr.br>: >

Re: [R-br] Ajuda com amostra pequena

2018-06-01 Por tôpico Cesar Rabak via R-br
Caro Elias, Antes de mergulhar no poço dos dados, vale a seguinte pergunta (que veja só você é *decorrente* da filosofia esposada por Virgílio) : « Qual foi o processo que gerou esses dados? » A amostra foi AAS? Como cada nível de cada fator foi determinado/calculado? Deixo-o com o seguinte

Re: [R-br] INDICADORES SOCIAIS OF TOPICS

2018-06-01 Por tôpico Cesar Rabak via R-br
Eucilene, Embora eu concorde *de plano* com sua colocação a respeito dos dados como primeiro passo para discutir um possível tratamento de dados, gostaria de enriquecer o debate, aproveitando a questão posta pelo OP, dizendo que em todas as ciências que se servem de Estatística há, em adição ao

Re: [R-br] INDICADORES SOCIAIS OF TOPICS

2018-05-28 Por tôpico Cesar Rabak via R-br
Edson, Malgrado a definição do quê seria uma técnica específica versus as mais comumente empregadas¹, vale a pena ver a *task view* Social Sciences no CRAN! HTH -- Cesar Rabak [1] Citando diretamente do link sugerido: « *Social scientists use a wide range of statistical methods, most of

Re: [R-br] Error: argument "ast" is missing, with no default (solve_model)

2018-05-28 Por tôpico Cesar Rabak via R-br
Em vez de usar o pipe você não devia chamar solve_model com a especificação do seu modelo, no seu caso *model*? On Mon, May 28, 2018 at 3:01 PM, Adriele Giaretta Biase via R-br < r-br@listas.c3sl.ufpr.br> wrote: > Boa tarde pessoal, > > estou trabalhando com modelos matemáticos de otimização

Re: [R-br] quebra de linha do título de tabelas geradas pela função xtable

2018-05-27 Por tôpico Cesar Rabak via R-br
Paulo, Você já recebeu sugestões, que pela natureza tendem a ser "quebra galhos" no grande esquema das coisas do LaTeX. As definições em relação à largura das *caption* tanto para tabelas como figuras, ou no idioma do LaTeX *floats* *deveria* ser dado pelo modelo de documento que você ou

Re: [R-br] Extraindo os raios e área de uma elipse ellipsoidhull() do pacote cluster

2018-05-24 Por tôpico Cesar Rabak via R-br
OK Alexandre! De fato, como eu não tenho conhecimento do domínio do problema que você estuda, me soa "não intuitivo" que se use uma ferramenta para medir uma área de influência sem levar em conta alguma 'resposta' dos tipos de tratamento... Boa sorte na busca de um pacote ou funções no R para

Re: [R-br] Extraindo os raios e área de uma elipse ellipsoidhull() do pacote cluster

2018-05-24 Por tôpico Cesar Rabak via R-br
Os valores retornados por EpointsA$d2 e EpointsB$d2 " average squared radius. " cfme. a documentação não lhe são suficientes? 2018-05-23 11:14 GMT-03:00 ASANTOS via R-br : > Caros Membros, > > Com fins ilustrativos e exploratório, tenho pontos correspondentes a

Re: [R-br] Modelo linear de segundo grau

2018-05-24 Por tôpico Cesar Rabak via R-br
André, O seu modelo não deveria ser especificado como: *> m3=rsm(y~ FO(x)+PQ(x))* > m3 Call: rsm(formula = y ~ FO(x) + PQ(x)) Coefficients: (Intercept)FO(x)PQ(x) -0.1333 6.7190 -0.8810 HTH -- Cesar Rabak 2018-05-21 14:45 GMT-03:00 Andre Oliveira via R-br

Re: [R-br] Emprego da PCA para seleção de variáveis mais correlacionadas no primeiro componente principal

2018-05-05 Por tôpico Cesar Rabak via R-br
Alexandre, Eu usei as funções do pacote "factoextra" . A partir do seu CMR você pode seguir com: > install.packages("factoextra") ## só vai ser necessário uma única vez >

Re: [R-br] Emprego da PCA para seleção de variáveis mais correlacionadas no primeiro componente principal

2018-05-04 Por tôpico Cesar Rabak via R-br
Sim Alexandre. Aliás, eu recomendo que você teste os resultados da pesquisa com os comandos (para os dados e resultados da PCA do seu CMR): > biplot(pca.object) e > biplot(pca.object, choice=2:3) Nessa representação você pode ver como as variáveis estão correlacionadas em relação às

Re: [R-br] Emprego da PCA para seleção de variáveis mais correlacionadas no primeiro componente principal

2018-05-02 Por tôpico Cesar Rabak via R-br
E uma observação que pode ajudá-lo a entender melhor o problema (ou a solução): Você escreveu: « Gostaria de selecionar as 5 variáveis resposta no primeiro componente principal mais correlacionadas . . . », mas o que o procedimento descrito no link faz é "Selecionar as cinco variáveis mais

Re: [R-br] Emprego da PCA para seleção de variáveis mais correlacionadas no primeiro componente principal

2018-05-02 Por tôpico Cesar Rabak via R-br
Você viu o resumo na primeira resposta no link do SO que você postou? On Wed, May 2, 2018 at 6:02 PM, ASANTOS via R-br wrote: > Prezados Membros, > > Gostaria de selecionar as 5 variáveis resposta no primeiro componente > principal mais correlacionadas usando a

Re: [R-br] Plot_ly background

2018-04-13 Por tôpico Cesar Rabak via R-br
Jobenil, Você já tentou fazer e por isso tem um CMR, especialmente para os dados do exemplo? 2018-04-10 16:55 GMT-03:00 Jobenil Luiz Magalhães Júnior via R-br < r-br@listas.c3sl.ufpr.br>: > Prezados(as), boa tarde > > Gostaria de construir um gráfico usando o pacote plotly de acordo com o >

Re: [R-br] Chi-quadrado para amostras grandes

2018-04-06 Por tôpico Cesar Rabak via R-br
Caro Elias, De fato o número grande de registros permite que eventuais correlações não importantes do ponto de vista do domínio do problema tornem-se *estatisticamente* "significantes". Não é um "problema" que possa ser resolvido por meio da matemática, mas sim do raciocínio no domínio do

Re: [R-br] Reproduçao dos resultados de Gabriel(1971)

2018-04-06 Por tôpico Cesar Rabak via R-br
Luiz, SE a pergunta for gerar exatamente a figura 2 que usa convenções e forma gráfica da época do artigo, provavelmente ninguém a reproduzirá com o R, agora se a pergunta é obter um biplot no R, há vários pacotes para essa finalidade! Qual a sua dificuldade para gerar o biplot? Aliás, você tem

Re: [R-br] Margem de gráfico

2018-04-05 Por tôpico Cesar Rabak via R-br
André: Considerando que o CMR é um exemplo do próprio pacote, *e* que o efeito parece ser decorrente da divisão por 500 da variável dependente: require(Publish) data(Diabetes) x=ci.mean(bp.2s~frame,data=Diabetes) plot(x) *Acho *que você deve reportar o problema para os mantenedores do pacote

Re: [R-br] Erro ao fazer boxplot - Arquivo SPSS

2018-03-26 Por tôpico Cesar Rabak via R-br
Edmar, Você posta os comandos mas não as respostas do R! Por favor envie-nos também a resposta aos comandos: > head(age) > summary(age) 2018-03-26 11:58 GMT-03:00 Edmar Caldas via R-br : > Eu estou usando RStudio Version 1.0.136. > > eu só fiz isso até agora. > > >

Re: [R-br] Dúvida TRUE/FALSE

2018-03-23 Por tôpico Cesar Rabak via R-br
Daniele, Sem o CMR é difícil dizer *exatamente* o quê pode estar acontecendo! Quando você compara um vetor, o resultado é um vetor: > set.seed(0) > > > vetor <- round(abs(rnorm(10) * 200)) > > > vetor <= 50 > > [1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE TRUE FALSE > > > > > SE você

Re: [R-br] comparação entre: ginv(x) com solve( t(x)%*% x ) %*% t(x)

2018-03-17 Por tôpico Cesar Rabak via R-br
A resposta a essa pergunta se centra na definição de "igualdade": 1) SE você deseja a igualdade com acurácia numérica exata, a resposta é "nunca"; 2) SE você aceita uma "faixa" de "precisão" dos resultados, a resposta é "na maioria dos casos". Para mais detalhes procure (tanto nos arquivos

Re: [R-br] Instalação do pacote brms

2018-03-14 Por tôpico Cesar Rabak via R-br
Pelas mensagens dá a impressão que o pacote rstan requer bibliotecas não FOSS (não livres no sentido Free Software) e seu OS não está habilitando a instalação delas... 2018-03-14 7:53 GMT-03:00 Andre Oliveira via R-br : > boa dia, > estou a procura de um pacote para

Re: [R-br] Interpretação dos dados

2018-03-09 Por tôpico Cesar Rabak via R-br
Fúlvia, Você realizou um teste que os gringos denominam *omnibus test* "não paramétrico", ou seja com relaxamento de certos pressupostos e obteve um resultado de valor-p menor que 0,05. Como você tem quatro tratamentos (df=3), então a segunda questão seria saber se *todos* os tratamento seriam

Re: [R-br] PCA

2018-03-09 Por tôpico Cesar Rabak via R-br
Veja estas obras: http://math.agrocampus-ouest.fr/infoglueDeliverLive/enseignement/support2cours/livres/ana.do.R https://livre.fnac.com/a9664406/Jerome-Pages-Analyses-factorielles-simples-et-multiples-5e-ed-Cours-et-etudes-de-cas

Re: [R-br] AJUDA: Teste de Homogeneidade

2018-03-08 Por tôpico Cesar Rabak via R-br
Lucas, Sem um CMR (v. o guia de postagem a respeito) este seu último e-mail vira o que os *nerds* denominaram jocosamente de « *confessional debugging* »! Por óbvio se você a variável contiver dados advindos dos resíduos de uma ANOVA o teste testará apenas eles. Por outro lado se a variável

Re: [R-br] Empilhar arquivos de um diretório em arquivo único

2018-02-26 Por tôpico Cesar Rabak via R-br
Yury, Como a resposta do Luis é uma solução diferente, apenas comento o seguinte: *caveat emptor* raciocínio em cima do código sem testar. A linha de código: * for(i in 1:length(dir(paste0(raiz, x{ * Ao ler o diretório n vezes, na segunda passada passará a ver um arquivo adicional

Re: [R-br] Experimento dados de contagem

2018-02-15 Por tôpico Cesar Rabak via R-br
Bom! Estamos convergindo! Os tratamentos são completamente categóricos ou podem se considerados alguma forma de "dosagem" de alguma substância ou agente? 2018-02-14 21:10 GMT-02:00 Maurício Lordêlo : > Isso Cesar. Considerar apenas uma variável dependente para este >

Re: [R-br] Experimento dados de contagem

2018-02-13 Por tôpico Cesar Rabak via R-br
OK, então a unidade experimental são larvas e o desfecho são as que morrem (ou o oposto as que morrem), portanto há no meu entender por essa nova informação apenas uma variável dependente, correto? 2018-02-11 16:07 GMT-02:00 Maurício Lordêlo via R-br < r-br@listas.c3sl.ufpr.br>: > > Fiz os

Re: [R-br] Experimento dados de contagem

2018-02-13 Por tôpico Cesar Rabak via R-br
Olá Maurício, A minha dúvida é que se descreve o experimento com palavras como "... mesma unidade experimental..." e pelo meu parco entendimento de entomologia pupa e larva são estágios diferentes da vida do inseto, por isso a minha questão era se poder-se-ia considerar como dois experimentos em

Re: [R-br] Experimento dados de contagem

2018-02-10 Por tôpico Cesar Rabak via R-br
Só para testar meu entendimento, as pupas e as larvas são objetos *diferentes*, então formalmente poder-se-ia dizer que foram dois experimentos feitos "em paralelo", correto? 2018-02-10 20:06 GMT-02:00 Maurício Lordêlo via R-br < r-br@listas.c3sl.ufpr.br>: > Caros, > Foi realizado um

Re: [R-br] Salvar saídas de um loop

2018-02-09 Por tôpico Cesar Rabak via R-br
Yuri, Seu *loop* na variável de indução "city", tem a chamada à write.table no escopo da variável "meses". Faça a seguinte correção e veja se os resultados fazem sentido para você: . #for(e in 1:length(estrategias)){ for(m in 1:length(meses)){ a = read.table(paste0(raiz, 'E1\\',

Re: [R-br] Assunto OFF

2018-02-05 Por tôpico Cesar Rabak via R-br
OK. Você quer criar uma curva de "calibração" do processo. Até aqui estou usando a palavra "processo" no sentido industrial do termo. O quê você precisa fazer Edson, é modelá-lo em um *processo estatístico.* A menos que o produto (no sentido da saída do processo "industrial") seja uma novação

Re: [R-br] Assunto OFF

2018-02-05 Por tôpico Cesar Rabak via R-br
Prez@do Edson, SE você se refere a um "gráfico de controle" no sentido usado em Qualidade, gráfico de média com uma medida de dispersão e no caso geral, assume-se que o processo siga a "Lei Normal" (gaussiana). O problema estaria no seu processo que "visualmente" estaria "fora de controle". A

Re: [R-br] Ordenar os fatores do eixo X no plot.group()

2018-02-02 Por tôpico Cesar Rabak via R-br
Fernando, É um bom começo, mas ainda não fica igual ao plot.group. . .樂 2018-02-02 16:49 GMT-02:00 Fernando Antonio de souza via R-br < r-br@listas.c3sl.ufpr.br>: > Para ficar igual a saída da função plot.group é necessário formatar alguns > parametros par() > >

Re: [R-br] Ordenar os fatores do eixo X no plot.group()

2018-02-02 Por tôpico Cesar Rabak via R-br
Rodrigo, O comando plot (groups) do agricolae é um barplot especializado que pega as médias do objeto retornado por LSD.test (no seu caso o 'comparison'). Você pode plotar as médias achando dentro dele onde elas estão armazenadas. No caso do objeto em questão elas estão em comparison$groups e

Re: [R-br] Como trabalhar com medida repetida para um mesmo indivíduo

2018-01-16 Por tôpico Cesar Rabak via R-br
Carlos, Para responder essa pergunta você é necessário precisar como você medirá a resposta: em princípio, ao você mencionar que a variável resposta seria a "altura de alimentação", pareceria que essa variável seria contínua. Outro ponto importante é que a análise estatística que você vai fazer

Re: [R-br] Problema ao gerar um mapa com dados estatísticos

2018-01-12 Por tôpico Cesar Rabak via R-br
Edson, O que a chamada: > summary(am$idh2000) retorna? 2018-01-12 13:25 GMT-02:00 Edson Lira via R-br : > Caros amigos, abaixo os arquivos (links) que estou usando para gerar um > mapa do estado do Amazonas com > o IDH Municipal de 2000 (está 2005, mas é 2000). > >

Re: [R-br] Distribuição de probabilidade

2018-01-08 Por tôpico Cesar Rabak via R-br
Caro Felipe, Algumas sugestões no âmbito da possível mecânica de como fazer isso já foram passadas por outros aqui na lista, por isso vou fazer um comentário sobre outra parte da sua questão: "meu objetivo é verificar qual melhor a distribuição para os meus dados e, a partir dessa "melhor",

Re: [R-br] Dúvida em componentes principais - gráfico

2018-01-04 Por tôpico Cesar Rabak via R-br
Ana, Nos casos como o seu o "costume" é se fazer o biplot, ou representação similar, para as duas componentes com maiores autovalores, que "explicam" mais a variabilidade dos seus dados. Caso a terceira componente acabe tendo importância na explicação do domínio do problema você pode

Re: [R-br] Panel xyplot Lattice

2017-12-05 Por tôpico Cesar Rabak via R-br
Fernando, Como já escrevi, não é o que eu vejo na figura que você mesmo enviou ou na gerada no meu ambiente usando seu script: O nome da variável no eixo Y está rente à borda da figura, e na direita os tick marks estão perto da borda mas não tanto que um cropping manual não fosse suficiente. A

Re: [R-br] Panel xyplot Lattice

2017-12-04 Por tôpico Cesar Rabak via R-br
Fernando, Se você constrói uma imagem que tem que ter três quadrados na horizontal, para manter o aspecto 1:1 eles ocuparão ~ 1/3 da largura da imagem... não vejo como isso pode ser um problema... Se você deseja que ocupem melhor sua página no seu relatório, copie e cole a imagem numa página

Re: [R-br] Sobreposição de BoxPlots

2017-12-03 Por tôpico Cesar Rabak via R-br
Leonard, A Jasmine está tentando mostrar a distribuição do consumo (na verdade "rendimento") dos veículos dos carros do dataframe mtcars com uma "dupla fatorização": nº de cilindros e nº de marchas desses veículos. Quando ela aplica o parâmetro "varwidth" aos boxplots a largura variável deles

Re: [R-br] usar separador de milhar

2017-11-15 Por tôpico Cesar Rabak via R-br
Adriele, Essa sua pergunta ou é muito elementar, e o caminho para a resposta seria você ler o documento (ou assemelhado) R-intro, e a documentação das funções: cat(), sprintf() e a melhor para começar format(). OU tem algo complicado aí que você não conseguiu "passar" na sua pergunta... HTH --

Re: [R-br] Corrigir ausência de dados

2017-11-14 Por tôpico Cesar Rabak via R-br
Edson, Considerando a última msg, sobre o problema em geral: Não existe uma *única* fórmula para resolver esse problema. No jargão atual a tentativa de resolvê-lo, além do óbvio de deletar os casos com as variáveis faltantes ou considerar uma a uma, há, a assim denominada, *imputação*. A

Re: [R-br] Gráfico com médias e erro agrupando var. preditora através de duas categorias

2017-11-13 Por tôpico Cesar Rabak via R-br
Veja os comandos stack() para colocar os seus dados no jeito necessário e a função plotCI() do pacote gplots. HTH -- Cesar Rabak 2017-11-13 12:51 GMT-02:00 Carlos Peixoto Neto via R-br < r-br@listas.c3sl.ufpr.br>: > Boa tarde a todos. Peço um auxílio para o caso : uma espécie de peixe >

Re: [R-br] Correlação entre variáveis categoricas

2017-11-10 Por tôpico Cesar Rabak via R-br
Elias, Umas observações: - Ao você « ... e eliminar os relacionamentos de variáveis com baixa corrrelação para evitar conexões causais espurias. » vai deixar as com alta correlação *apenas* e aí ter "colinearidade" para as variáveis? - A recomendação: « Para variáveis categoricas x

Re: [R-br] Função do ExpDes.pt

2017-11-10 Por tôpico Cesar Rabak via R-br
André, Em adição ao que se pede costumeiramente (e que a lista automaticamente adiciona sobre o CMR), neste caso ajudaria se você fornecesse informações sobre a plataforma, versão do R e a versão do pacote em questão... O pacote datado de 2013, versão 1.1.2 não tem essa função na minha

Re: [R-br] Digest R-br, volume 83, assunto 2

2017-11-08 Por tôpico Cesar Rabak via R-br
O próprio pacote não vem com nenhum exemplo? 樂 2017-11-07 1:09 GMT-02:00 Wilson Cohen via R-br : > Olá Pessoal > Alguém poderia compartilhar um exemplo simples de uso do modelo não-linear > usando o pacote nlme? > Obrigado. > > Em 6 de novembro de 2017 11:00,

Re: [R-br] Erro em script de apostila Embrapa Arroz e Feijão

2017-10-23 Por tôpico Cesar Rabak via R-br
Tá com cara que o script esperava alguma variável já "setada" de outro script (anterior) ou alguma coisinha com versão mais antiga do R. Se você puder seguir o guia de postagem e postar um CMR. . . 2017-10-20 10:57 GMT-02:00 Euriann Yamamoto via R-br < r-br@listas.c3sl.ufpr.br>: > Bom dia,

Re: [R-br] Preencher área entre duas linhas

2017-10-17 Por tôpico Cesar Rabak via R-br
Mesmo que você defina a linha de polygon() como de cor nula (transparente)? On Mon, Oct 16, 2017 at 5:26 PM, Delcio Rudinei Bortolanza via R-br < r-br@listas.c3sl.ufpr.br> wrote: > Olá! > Segue um exemplo hipotético: > > #=== > x<-1:5 > > set.seed(4) >

Re: [R-br] Merge tabelas não padronizados.

2017-10-17 Por tôpico Cesar Rabak via R-br
Ronaldo, Sua questão está demasiadamente difusa para ser entendida. A frase « . . .em algum local do campo tem o texto que necessito. . . » não permite que se construa nenhuma especificação de como achar o quê você necessita. Tente explicar melhor! HTH 2017-10-10 15:29 GMT-03:00 Ronaldo

Re: [R-br] Regressão Exponencial

2017-09-18 Por tôpico Cesar Rabak via R-br
Fabiana, Se você não se importar de fazer um passo da matemática dessa solução "na mão", a regressão exponencial é normalmente feita sem se necessitar de pacotes mais sofisticados (nada contra o uso deles) empregando a transformação da variável pra tornar regressão num caso de regressão linear:

Re: [R-br] ajustar escala e interpolação

2017-09-11 Por tôpico Cesar Rabak via R-br
+1 2017-09-11 10:19 GMT-03:00 Jônatan via R-br : > Só gostaria de acrescentar que vento é um grandeza vetorial. Se for > interpolar (independente do método) para os dados de uma maior para menor > resolução (upscaling) considere a velocidade e a direção do vento (ou as

Re: [R-br] Modelo Misto - plot (augPred()) mensagem de erro

2017-09-09 Por tôpico Cesar Rabak via R-br
Você não recebeu nenhum erro na chamada a nlsList() não? Você tentou ver o resultado, com summary ou print? 2017-09-08 19:19 GMT-03:00 Fernando Souza via R-br : > Caros colegas ao rodar a função augPred {nlme} recebo a seguinte mensagem: > > "Error in

Re: [R-br] Estimar último valor na sequência de 0s e 1s...

2017-09-09 Por tôpico Cesar Rabak via R-br
Leonard, Como você é um habitual participante da lista, vou, de plano, descartar a hipótese que essa questão seja um exercício ou "lição de casa", malgrado ela tem todos os aspectos de uma 樂! Essa questão diz respeito de detectar padrões em sequências e a resposta está nesse domínio do

Re: [R-br] ajustar escala e interpolação

2017-09-09 Por tôpico Cesar Rabak via R-br
Heloíse, Você se fez entender e de fato você tem duas particularidades: seu grid e coordenadas de vento além de terem granularidades diferentes têm um ligeiro *offset* que se mostra como você descreve. Se as coordenadas do vento (que têm a maior granularidade) fossem arrendondadas para décimo de

Re: [R-br] ajustar escala e interpolação

2017-09-06 Por tôpico Cesar Rabak via R-br
Heloise, Esta observação me deixa confuso: #Quando faço um merge() dos dois dataframes, obviamente meus valores de > #vento$media aparecem como NA pois não há nenhuma correspondência entre > as #coordenadas do grid e de vento. As coordenadas do grid e vento múltiplas de 0,5 grau deveriam "se

Re: [R-br] Como ver o que uma função faz

2017-09-01 Por tôpico Cesar Rabak via R-br
Pedro, Acho que nesse caso a forma mais prática seria você baixar o código do pacote e ver a dita cuja diretamente. . . 2017-09-01 9:38 GMT-03:00 Pedro Emmanuel Alvarenga Americano do Brasil via R-br : > library(rms) > rms:::Function.cph > rms:::Function.rms > >

Re: [R-br] Barplot com dois eixos

2017-09-01 Por tôpico Cesar Rabak via R-br
Olá Diogo. De fato vendo a sua explicação fica menos claro o efeito que você deseja. Para não ficar um interminável vai e vem de trechos de código e avaliações que vão começar a parecer o jogo do "tá quente, tá frio!", acho que poderíamos fazer o seguinte: Qual o gráfico que você quer gerar com

Re: [R-br] Barplot com dois eixos

2017-08-31 Por tôpico Cesar Rabak via R-br
Usando seu exemplo: > barras <- barplot(dadosest,ylim=c(0,150),beside=TRUE) > points(barras[,1],dados$y1) > points(barras[,2],dados$y2, pch=19) > lines(barras[,1],dados$y1) > lines(barras[,2],dados$y2, lty=2) HTH -- Cesar Rabak PS.: Se você quer ver as "bolinhas e as linhas" em outro gráfico do

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

2017-08-25 Por tôpico Cesar Rabak via R-br
t; RSQLServer não é boa alternativa, não sofre manutenção a tempos e depende >> de Java na máquina. Sugiro o pacote odbc. >> >> SQL Server é mestre em dar problema no R comigo, mas o problema é na >> autenticação. >> >> Em 1 de ago de 2017 9:14 AM, "Zhu Sha

Re: [R-br] não parametrica - kruskal e dunn

2017-08-25 Por tôpico Cesar Rabak via R-br
Christina, O pacote *agricolae* não tem o teste de Dunn como parte do repertório de funções dele e estranhamente reporta na sua documentação que os "testes post-hoc não paramétricos, citando inclusive o *kruskal*, estão usando o critério LSD de Fisher". Você terá que pegar a função do teste de

Re: [R-br] RES: Dúvidas

2017-08-24 Por tôpico Cesar Rabak via R-br
Rosana: Se « Meus dados apresentam colunas informando o dia, o ano e o valor do atributo » onde eles indicam a parte espacial? Você parece estar "perdida" e não leu a documentação do pacote spacetime. Veja o artigo no Journal of Statistical Software Vol 52, nº 7 (2012) doi

Re: [R-br] Análise de dados antes e depois

2017-08-08 Por tôpico Cesar Rabak via R-br
Cleber, Sua pergunta extrapola o objetivo deste grupo (não CMR, etc.). Ademais ela enseja muitas questões antes que possa ser iniciada uma discussão com assertividade. Minha recomendação é que você a refaça nos foruns mais voltados para questões de Estatística (p.ex.:

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

2017-07-31 Por tôpico Cesar Rabak via R-br
Correto Pedro! Contudo, não entendo porque. . . Minha forma de atacar esse problema seria mandar os dados para o BD por outra maneira (que ñ via R) e tentar conectar apenas para puxar os dados: aí a gente vê se é na criação do BD via R que ele deixa em aberto um tamanho default muito grande ou

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

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

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 >

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

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

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

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

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

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

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] 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] 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] 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

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] 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

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

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] 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:

Re: [R-br] SOMAR LINHAS

2017-05-04 Por tôpico Cesar Rabak via R-br
que possamos refazer > e lhe ajudar aqui. > > Att > > Em 3 de maio de 2017 22:36, Cesar Rabak via R-br <r-br@listas.c3sl.ufpr.br > > escreveu: > >> Veja a documentação para o comando agregate, especialmente a interface >> com uso de "fórmulas". >

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

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

2017-04-24 Por tôpico Cesar Rabak via R-br
s/bootstrapping.html> > > > 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 sim

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

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] 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

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 <-

  1   2   >