[R-br] Superfície resposta

[Adriano Carvalho Costa por (R-br)]

Boa noite pessoal, estou com um artigo de aquaponia em que estamos fazendo
a viabilidade econômica de acordo com o preço de venda dos produtos, sendo
as variáveis respostas TIR, VPL entre outros. Estou precisando de um
parceiro que interesse em participar como co-autor para fazer o gráfico.
O trabalho faz parte de uma dissertação de mestrado. Quem tiver interesse
meu whats 64 981186749.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Análise de dados medidas repetidas no tempo

[Adriano Carvalho Costa via R-br]

Bom dia, solicioto o scrip fazendo favor

Em 15 de junho de 2016 19:27, Luiz Roberto Martins Pinto via R-br <
r-br@listas.c3sl.ufpr.br> escreveu:

> Eu enviei para várias pessoas um script R para análise de dados de medidas
> repetidas no tempo, elaborado pelo prof. Ivan Bezerra (UESC).
>
> Solicito o favor de informar:
> 1) Se o script foi relevante para as análises.
> 2) Se podem informar o script que efetivamente foi utilizado.
> 3) Novas sugestões de análises.
>
> Aos que não solicitaram o script, mas que podem colaborar, agradeço se o
> fizerem.
>
> Abraços,
>
> Luiz Roberto.
>
>
>
>
> Luiz Roberto Martins Pinto
> Prof. Pleno/DCET/UESC
> Laboratório de Estatística Computacional
> Universidade Estadual de Santa Cruz
> Ilhéus-Bahia-Brasil
>
> luizroberto.u...@gmail.com
> skype: lrmpinto
> http://lattes.cnpq.br/2732314327604831
>
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Vídeo Aula R

[Adriano Carvalho Costa]

Parabéns, ajudará muitas pessoas.

Em 19 de fevereiro de 2015 16:02, Augusto Ribas 
escreveu:

> Muito legal. Parabéns pela iniciativa.
>
> Em 19 de fevereiro de 2015 13:55, Guilherme Moraes Ferraudo <
> guiferra...@gmail.com> escreveu:
>
> Muito obrigado Teo.
>>
>> Vou atrás disto...
>>
>> Parabens.
>>
>> Abs
>>
>> Em 18 de fevereiro de 2015 16:33, Téo Calvo 
>> escreveu:
>>
>>> Olá novamente, fico feliz com o retorno positivo dos companheiro em
>>> relação às vídeo aulas.
>>>
>>> O Sistema Operacional que uso, é o Manjaro (GNU/Linux), o programa de
>>> gravação de vídeo da desktop é o Kazam.
>>>
>>> Sempre tento dar preferência para softwares livres, acredito que para
>>> usuários Windows, há o Camtasia, mas não sei informar se é proprietários.
>>>
>>> E muito obrigado novamente, sem dúvida serve como motivação. :)
>>>
>>> Abraços.
>>>
>>> ___
>>> R-br mailing list
>>> R-br@listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>
>>
>>
>> --
>> Guilherme Moraes Ferraudo
>> São Paulo - SP
>> http://lattes.cnpq.br/2096118558794430
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
>
> --
> Grato
> Augusto C. A. Ribas
>
> Site Pessoal: http://recologia.com.br/ 
> Github: https://github.com/Squiercg
> Lattes: http://lattes.cnpq.br/7355685961127056
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Ajuda em análise de trilha (path analysis)

[Adriano Carvalho Costa]

Alisson,

Faça a média de cada tratamento pra cada varíavel e use ela pra fazer as
análises. Fiz um stepwise pra manter as variáveis independentes com uma
correlação aceitavél para evitar a multicolinearidade.











Em 17 de dezembro de 2013 14:18, Gustavo Carvalho
escreveu:

> Veja o lavaan também. Eu gosto do livro "Cause and correlation in biology".
>
> 2013/12/17 Tiago Souza Marçal :
> > O pacote agricolae realiza a análise de trilha.
> >
> > Uma boa abordagem teórica você vai encontrar no livro "Modelos
> Biométricos
> > Aplicados ao melhoramento"
> >
> > Att.
> >
> > Tiago.
> >
> >
> > #
> >
> > Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES
> >
> > Bolsista de Iniciação Científica da área de Genética e Melhoramento de
> > Plantas
> >
> > #
> >
> >
> > 
> > Date: Mon, 16 Dec 2013 13:55:14 -0800
> > From: alisson...@yahoo.com.br
> > To: r-br@listas.c3sl.ufpr.br
> > Subject: [R-br] Ajuda em análise de trilha (path analysis)
> >
> >
> > Caro Colegas,
> > Boa noite.
> >
> > Eu estou fazendo uma análise de trilha (path analysis) para um colega,
> > entretanto nunca realizei uma análise desse tipo. Assim gostaria de uma
> > ajuda de vocês com sugestões e materiais para consulta. Esse colega
> conduziu
> > um experimento em fatorial com 4 níveis de solo e 6 níveis de amenizantes
> > para metais pesados com cultivo de planta indicadoras. Ele determinou os
> > metais nos solos e mediu os parâmetros de crescimento das plantas que
> foram
> > cultivadas nesses solos com os amenizantes. Ele estudou relações causa e
> > efeito com as variáveis mediadas.
> >
> > Obrigado.
> >
> > Alisson Lucrécio da Costa
> >
> > ___ R-br mailing list
> > R-br@listas.c3sl.ufpr.br
> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de
> > postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo
> > reproduz�vel.
> >
> > ___
> > R-br mailing list
> > R-br@listas.c3sl.ufpr.br
> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código
> > mínimo reproduzível.
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] ANOVA

[Adriano Carvalho Costa]

anova(fit)


Em 4 de novembro de 2013 15:13, geovane barbosa
escreveu:

>
> Olá pessoal tudo bem, boa tarde
>
>
>
> Estou querendo realizar uma ANOVA em um conjunto de dados, o código segue
> abaixo
>
>
> fit <- aov(Pressão ~  Fumante, data=mydataframe)
>
>
>
> Nesse caso ele faria anova entre o grupo de fumantes (3 categorias) do
> grupo Raça(3 categorias).
>
>
>
> como posso fazer esse procedimento?
>
>
>
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Correlações

[Adriano Carvalho Costa]

Boa noite,

Estou estudando sobre correlações, alguém ai tem material sobre correlação:
 - Variáveis qualitativo e variáveis qualitativa
 - Variáveis discretas e variáveis discretas
- Variáveis discretas e variáveis qualitativas
- Variáveis qualitativas e variáveis continuas
- Variáveis discretas e variáveis contínuas

Obrigado!
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Escalas

[Adriano Carvalho Costa]

Eric

plot(c(5,20), c(8,30), type='n')

 demo(graphics)


*Adriano Carvalho Costa*
Doutorando em Produçao Animal / Melhoramento - UFLA
Coordenador do Grupo de Melhoramento Animal e Biotecnologia






Em 10 de maio de 2012 09:53, Paulo Justiniano escreveu:

> ylim no plot()
>
>
>
> On Thu, 10 May 2012, Eric Keven wrote:
>
>  Boa noite,
>> Como eu faço para alterar a escala dos gráficos e fazer com que o programa
>> comece num valor que eu fixar para a escala?
>> Obrigado!
>>
>> Att,
>>
>>
>>
>>
>>
>> Eric Keven Silva
>>
>> Coordenador do Núcleo de Estudos de Processos da Indústria de Alimentos
>> (PROALI)
>>
>> Graduando em Engenharia de Alimentos
>>
>> Universidade Federal de Lavras
>>
>>
>>
>>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Contraste

[Adriano Carvalho Costa]

Boa noite,

Gostaria de saber se alguém tem um scrip para contraste.

Obrigado!
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Obtenção do ponto de inflexão em modelo logístico

[Adriano Carvalho Costa]

Boa noite,

Venho solicitar a ajuda de voces, pois estou tentando calcular o ponto de
inflexão de um modelo logístico e não estou conseguindo, tenho realizado os
seguintes comandos.

> ajnll.gg1<-nls(PESO ~
A/(1+B*exp(-k*IDADE)),data=gg1.dat,start=c(A=900,B=1,k=0.01))
> s.ajnll.gg1<-summary(ajnll.gg1)
> s.ajnll.gg1
Formula: PESO ~ A/(1 + B * exp(-k * IDADE))
Parameters:
   Estimate Std. Error t value Pr(>|t|)
A 9.284e+02  1.322e+01  70.202  < 2e-16 ***
B 1.801e+01  2.210e+00   8.148 2.68e-14 ***
k 2.826e-02  1.337e-03  21.140  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 81.91 on 222 degrees of freedom
Number of iterations to convergence: 8
Achieved convergence tolerance: 2.787e-06

> logi <- expression(A/(1+B*exp(-K*x)))
>  D(logi,'x')
A * (B * (exp(-K * x) * K))/(1 + B * exp(-K * x))^2

Substituindo na equação
 > GGtc= 928.4*(18.01*(exp(-0.028*x)*0.028))/(1+18.01*exp(-0.028*x))^2
> GGtc
[1] 1.328331 1.361869 1.396139 1.431151 1.466914
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Teste de Dunnett

[Adriano Carvalho Costa]

Boa noite,

Venho pedir ajuda de voces para perguntar se alguém conhece alguma rotina
para realizar um teste de Dunnet. Pois estou com um experimento fatorial
com um tratamento adicional e devo comparar o adicional versus os demais
(um a um).

Muito obrigado!
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Fwd: Duvida para an[alise

[Adriano Carvalho Costa]

Walmes, completando o ultimo e-mail...

Para Prop. Ca e CaCl eu tenho um nível com zero.

Em 19 de dezembro de 2011 19:48, Adriano Carvalho Costa <
acarvalhoco...@gmail.com> escreveu:

> Walmes,
>
> Mas não é recomendado regressão acima de quatro pontos (doses ou níveis)?
>
> Em 19 de dezembro de 2011 18:54, Walmes Zeviani 
> escreveu:
>
>>  Se a análise dos resíduos não recusar o modelo, pode usar para predição
>> sim, claro, para valores dentro do intervalo investigado.
>>
>> À disposição.
>> Walmes.
>>
>>
>> ==
>> Walmes Marques Zeviani
>> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
>> Departamento de Estatística - Universidade Federal do Paraná
>> fone: (+55) 41 3361 3573
>> VoIP: (3361 3600) 1053 1173
>> e-mail: wal...@ufpr.br
>> twitter: @walmeszeviani
>> homepage: http://www.leg.ufpr.br/~walmes
>> linux user number: 531218
>> ==
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Fwd: Duvida para an[alise

[Adriano Carvalho Costa]

Walmes,

Mas não é recomendado regressão acima de quatro pontos (doses ou níveis)?

Em 19 de dezembro de 2011 18:54, Walmes Zeviani
escreveu:

> Se a análise dos resíduos não recusar o modelo, pode usar para predição
> sim, claro, para valores dentro do intervalo investigado.
>
> À disposição.
> Walmes.
>
>
> ==
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: wal...@ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Fwd: Duvida para an[alise

[Adriano Carvalho Costa]

Walmes,

Neste caso eu vou fazer a regressão com objetivo comparar os betas apenas.
O modelo gerado não vou utilizar pois tenho apenas três pontos. Certo?






















































































































































Em 19 de dezembro de 2011 13:58, Walmes Zeviani
escreveu:

> Walmes,
>
> Acho esse conjunto de contrastes uma maneira de complicar o simples. Eu
> ajustataria modelos de regressão da contínua separado pela categórica e
> compararia os interceptos e inclinações, ou apenas a colocaria em gráfico
> com bandas de confiança para o valor predito. Com as devidas considerações,
> um output mais ou menos assim
> http://ridiculas.wordpress.com/2011/08/19/uma-alternativa-para-apresentar-bandas-de-confianca-em-regressao/
> .
>
> À disposição.
> Walmes.
>
> ==
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: wal...@ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Fwd: Duvida para an[alise

[Adriano Carvalho Costa]

-- Mensagem encaminhada --
De: Adriano Carvalho Costa 
Data: 19 de dezembro de 2011 13:03
Assunto: Duvida para an[alise
Para: r-br-requ...@listas.c3sl.ufpr.br


Boa tarde,
Estou com uma duvida para analisar este experimento. Alguem pode me ajudar.
Penso da seguinte forma.

Contrastes
1. Pectina e Amido vs Prop.Ca e CaCl
2. Pectina vs Amido
3. Pectina 1,5 vs Pectina 7,5 e Pectina 15
4. Pectina 7,5 vs Pectina 15
5. Amido 15 vs amido 45 e amido 75
5. Amido 45 vs amido 75
6. Prop. Ca vs CaCl
7. Prop. Ca 0 vc Prop. Ca 15  e Prop. Ca 30
8. Prop. Ca 15 vs Prop. Ca 30
9. CaCl 0 vc CaCl 15 e CaCl30
10.CaCl 15 vs CaCl Ca 30

  tratamento Concentrações  Pectina 1,5  Pectina 7,5  Pectina 15  amido 15
amido 45  amido 75  Prop. Ca 0  Prop. Ca 15  Prop. Ca 30  CaCl 0  CaCl 15
CaCl 30
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dados sem distribuição

[Adriano Carvalho Costa]

Cezar, utilizei box-cox, pois ela já resolveu algumas vezes.

Obrigado!!

Em 15 de agosto de 2011 19:06, Cesar Rabak  escreveu:

> Os seus dados têm uma distribuição assimétrica se descontarmos a grande
> frequência na classe (0;10], ou melhor dizendo ela me parece uma mistura de
> duas distribuições.
>
> Qual é o pressuposto científico (não estatístico) que lhe motivou a tenta a
> transf. Box-Cox?
>
>
> Em 15/8/2011 18:25, Adriano Carvalho Costa escreveu:
>
>>  Boa noite,
>> Estou com a base de dados acima e estou fazendo transformação box-cox e
>> os dados não atenderam a distribuição normal, voces tem alguma rotina de
>> análise não paramétrica?
>> Script
>> dados<-read.table("final.txt", h=T)
>> dados
>> names(dados)
>> attach(dados)
>> head(dados)
>> dim(dados)
>> Transformação Box-Cox
>> require(MASS)
>> library(MASS)
>> motil
>> Trat
>> boxcox(motil+0.01~Trat, data=dados)
>> lambida=0.35
>> motili=(motil^(0.39) - 1)/0.39
>> motili
>> plot(motil~Trat)
>> plot(motili~Trat)
>> anava <- aov(motili ~ Trat, data = dados)
>> plot(anava)
>> shapiro.test(anava$res)
>>
>>
>> __**_
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-br<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br>
>> Leia o guia de postagem 
>> (http://www.leg.ufpr.br/r-br-**guia<http://www.leg.ufpr.br/r-br-guia>)
>> e forneça código mínimo reproduzível.
>>
>
>
> --
> Cesar Rabak
> GNU/Linux User 52247.
> Get counted: http://counter.li.org/
> __**_
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-br<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br>
> Leia o guia de postagem 
> (http://www.leg.ufpr.br/r-br-**guia<http://www.leg.ufpr.br/r-br-guia>)
> e forneça código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dados sem distribuição

[Adriano Carvalho Costa]

Walmes, muito obrigado pela orientação,

Dados http://www.datafilehost.com/download-d2b09ba3.html.

Nunca trabalhei com glm, como posso continuar a análise.

Obrigado!


Em 15 de agosto de 2011 19:18, Adriano Carvalho Costa <
acarvalhoco...@gmail.com> escreveu:

> Walmes, o Trat são todos os tratamentos que eu tenho, no caso todas as
> combinações possiveis mais o adicional.
>
> A variavel resposta motilidade realmente esta em porcentagem com escalas
> que mudam de 5 em 5 , e vai de 0 a 100 e a varial quali tem 6 niveis de 0 a
> 6 mesmo.
>
>
>
> Em 15 de agosto de 2011 19:02, Walmes Zeviani 
> escreveu:
>
>>  Adriano,
>>
>> No seu *.txt o fator Trat é númerico. Você converteu para fator para rodar
>> as análises? No seu script não tem conversão. Outro ponto é que sua variável
>> motil tem cara de ser uma porcentagem, números multiplos de 5 entre 0 e 100,
>> o que sugere um número de 20 ensaios bernoulli. Se for o caso, poderia-se
>> usar glm binomial. A sua variável quali parece ser nota numa escala de 6
>> níveis (0 à 6) e que por essa característica pode não atender os
>> pressupostos da anova. Para tornar seu realmente código reproduzível (sem
>> etapa de download), envie o dput() dos dados (poucos dados), ou hospede os
>> dados passando o link (muitos dados), e não dados em anexo. Veja o guia de
>> postagem (rodapé da mensagem) para instruções de como enviar dados.
>>
>> À disposição.
>> Walmes.
>>
>> ==
>> Walmes Marques Zeviani
>> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
>> Departamento de Estatística - Universidade Federal do Paraná
>> fone: (+55) 41 3361 3573
>> VoIP: (3361 3600) 1053 1173
>> e-mail: wal...@ufpr.br
>> twitter: @walmeszeviani
>> homepage: http://www.leg.ufpr.br/~walmes
>> linux user number: 531218
>> ==
>>
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dados sem distribuição

[Adriano Carvalho Costa]

Walmes, o Trat são todos os tratamentos que eu tenho, no caso todas as
combinações possiveis mais o adicional.

A variavel resposta motilidade realmente esta em porcentagem com escalas que
mudam de 5 em 5 , e vai de 0 a 100 e a varial quali tem 6 niveis de 0 a 6
mesmo.



Em 15 de agosto de 2011 19:02, Walmes Zeviani escreveu:

> Adriano,
>
> No seu *.txt o fator Trat é númerico. Você converteu para fator para rodar
> as análises? No seu script não tem conversão. Outro ponto é que sua variável
> motil tem cara de ser uma porcentagem, números multiplos de 5 entre 0 e 100,
> o que sugere um número de 20 ensaios bernoulli. Se for o caso, poderia-se
> usar glm binomial. A sua variável quali parece ser nota numa escala de 6
> níveis (0 à 6) e que por essa característica pode não atender os
> pressupostos da anova. Para tornar seu realmente código reproduzível (sem
> etapa de download), envie o dput() dos dados (poucos dados), ou hospede os
> dados passando o link (muitos dados), e não dados em anexo. Veja o guia de
> postagem (rodapé da mensagem) para instruções de como enviar dados.
>
> À disposição.
> Walmes.
>
> ==
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: wal...@ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Dados sem distribuição

[Adriano Carvalho Costa]

Esqueci de falar, é um experimento com dois fatores (dilui e osm) e com um
tratamento adicional

Em 15 de agosto de 2011 18:25, Adriano Carvalho Costa <
acarvalhoco...@gmail.com> escreveu:

> Boa noite,
>
> Estou com a base de dados acima e estou fazendo transformação box-cox e os
> dados não atenderam a distribuição normal, voces tem alguma rotina de
> análise não paramétrica?
>
> Script
>
> dados<-read.table("final.txt", h=T)
> dados
> names(dados)
> attach(dados)
> head(dados)
> dim(dados)
> Transformação Box-Cox
> require(MASS)
> library(MASS)
> motil
> Trat
> boxcox(motil+0.01~Trat, data=dados)
> lambida=0.35
> motili=(motil^(0.39) - 1)/0.39
> motili
> plot(motil~Trat)
> plot(motili~Trat)
> anava <- aov(motili ~ Trat, data = dados)
> plot(anava)
> shapiro.test(anava$res)
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Dados sem distribuição

[Adriano Carvalho Costa]

Boa noite,

Estou com a base de dados acima e estou fazendo transformação box-cox e os
dados não atenderam a distribuição normal, voces tem alguma rotina de
análise não paramétrica?

Script

dados<-read.table("final.txt", h=T)
dados
names(dados)
attach(dados)
head(dados)
dim(dados)
Transformação Box-Cox
require(MASS)
library(MASS)
motil
Trat
boxcox(motil+0.01~Trat, data=dados)
lambida=0.35
motili=(motil^(0.39) - 1)/0.39
motili
plot(motil~Trat)
plot(motili~Trat)
anava <- aov(motili ~ Trat, data = dados)
plot(anava)
shapiro.test(anava$res)
Tratdilui   rep osm motil   quali
1   nacl1   45  95  5
2   nacl1   90  95  4
3   nacl1   135 90  4
4   nacl1   180 90  4
5   nacl1   225 95  4
6   nacl1   270 85  4
7   nacl1   315 70  3
8   nacl1   360 1   1
9   nacl1   405 0   0
10  nacl1   450 0   0
1   nacl2   45  90  4
2   nacl2   90  95  5
3   nacl2   135 95  5
4   nacl2   180 95  5
5   nacl2   225 90  4
6   nacl2   270 85  4
7   nacl2   315 65  4
8   nacl2   360 0   0
9   nacl2   405 0   0
10  nacl2   450 0   0
1   nacl3   45  85  4
2   nacl3   90  95  4
3   nacl3   135 90  5
4   nacl3   180 95  5
5   nacl3   225 85  4
6   nacl3   270 75  4
7   nacl3   315 70  4
8   nacl3   360 10  2
9   nacl3   405 0   0
10  nacl3   450 0   0
1   nacl4   45  75  4
2   nacl4   90  80  4
3   nacl4   135 85  4
4   nacl4   180 90  5
5   nacl4   225 95  5
6   nacl4   270 95  5
7   nacl4   315 25  3
8   nacl4   360 0   0
9   nacl4   405 0   0
10  nacl4   450 0   0
1   nacl5   45  60  4
2   nacl5   90  75  4
3   nacl5   135 85  4
4   nacl5   180 95  5
5   nacl5   225 95  5
6   nacl5   270 75  5
7   nacl5   315 50  4
8   nacl5   360 5   2
9   nacl5   405 0   0
10  nacl5   450 0   0
1   nacl6   45  80  4
2   nacl6   90  85  4
3   nacl6   135 90  4
4   nacl6   180 95  5
5   nacl6   225 85  5
6   nacl6   270 65  4
7   nacl6   315 60  4
8   nacl6   360 0   0
9   nacl6   405 0   0
10  nacl6   450 0   0
1   nacl7   45  65  4
2   nacl7   90  75  4
3   nacl7   135 80  4
4   nacl7   180 85  4
5   nacl7   225 80  4
6   nacl7   270 70  4
7   nacl7   315 50  3
8   nacl7   360 1   1
9   nacl7   405 0   0
10  nacl7   450 0   0
11  glicose 1   45  100 5
12  glicose 1   90  100 5
13  glicose 1   135 95  4
14  glicose 1   180 95  4
15  glicose 1   225 90  4
16  glicose 1   270 90  4
17  glicose 1   315 70  3
18  glicose 1   360 0   0
19  glicose 1   405 0   0
20  glicose 1   450 0   0
11  glicose 2   45  90  4
12  glicose 2   90  95  5
13  glicose 2   135 95  5
14  glicose 2   180 95  5
15  glicose 2   225 95  5
16  glicose 2   270 90  5
17  glicose 2   315 60  4
18  glicose 2   360 0   0
19  glicose 2   405 0   0
20  glicose 2   450 0   0
11  glicose 3   45  80  4
12  glicose 3   90  90  4
13  glicose 3   135 95  4
14  glicose 3   180 80  4
15  glicose 3   225 80  4
16  glicose 3   270 50  3
17  glicose 3   315 45  3
18  glicose 3   360 5   2
19  glicose 3   405 0   0
20  glicose 3   450 0   0
11  glicose 4   45  90  4
12  glicose 4   90  95  4
13  glicose 4   135 95  5
14  glicose 4   180 95  5
15  glicose 4   225 95  5

Re: [R-br] Experimento fatorial com tratamento adicional sem normalidade

[Adriano Carvalho Costa]

Muito obrigado Ivan!

Em 27 de julho de 2011 13:26, Ivan Bezerra Allaman
escreveu:

>  Adriano,
>
> Esse é o caminho das pedras,
>
> http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/
>
> Allaman
> (S,f,P)
>
> *Ivan Bezerra Allaman*
> Doutor em Zootecnia/UFLA
> email e msn - ivanala...@yahoo.com.br
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Experimento fatorial com tratamento adicional sem normalidade

[Adriano Carvalho Costa]

Ivan,

fa <- read.table("http://www.leg.ufpr.br/~walmes/cursoR/fat-adi.txt";,
header=TRUE)
attach(fa)
fa
names(fa)
concentração=concentraÃ.Ã.o
fa <- transform(fa, concentração=factor(concentração), bloc=factor(bloc))
m0 <- aov(media~bloc+trat, fa)
anova(m0)
boxplot (media~bloc+trat, data=fa)
library(MASS)
boxcox(media+0.01~bloc+trat,data=fa)
locator()
#Agora basta elevar cada valor da unidade resposta a 1,75?



Em 27 de julho de 2011 11:38, Adriano Carvalho Costa <
acarvalhoco...@gmail.com> escreveu:

> Ivan, não tentei a metodologia boxcox, voce pode me passar ela?
>
> Em 27 de julho de 2011 07:46, Jose Claudio Faria <
> joseclaudio.fa...@gmail.com> escreveu:
>
> Concordo com a sugestão do Ivan!
>>
>> Abs,
>> --
>> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>> Jose Claudio Faria
>> Estatistica - Prof. Pleno
>> UESC/DCET/Brasil
>> joseclaudio.faria at gmail.com
>> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>>
>> Em 27 de julho de 2011 07:36, Ivan Bezerra Allaman
>>  escreveu:
>>  > Adriano,
>> > Se a variável não tem distribuição normal e você não quer utilizar
>> nenhuma
>> > transformação, terá que recorrer a metodologia dos Modelos Lineares
>> > Generalizados (GLM) e encontrar a família de distribuição que se encaixe
>> aos
>> > seus dados. Eu o aconselho tentar uma transformação primeiro ser você
>> não
>> > estiver familiarizado com a metodologia GLM. Já tentou boxcox?
>> > Allaman
>> > (S,f,P)
>> >
>> > Ivan Bezerra Allaman
>> > Doutor em Zootecnia/UFLA
>> > email e msn - ivanala...@yahoo.com.br
>> >
>> > ___
>> > R-br mailing list
>> > R-br@listas.c3sl.ufpr.br
>> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código
>> > mínimo reproduzível.
>> ___
>> R-br mailing list
>> R-br@listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] Experimento fatorial com tratamento adicional sem normalidade

[Adriano Carvalho Costa]

Ivan, não tentei a metodologia boxcox, voce pode me passar ela?

Em 27 de julho de 2011 07:46, Jose Claudio Faria <
joseclaudio.fa...@gmail.com> escreveu:

> Concordo com a sugestão do Ivan!
>
> Abs,
> --
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
> Jose Claudio Faria
> Estatistica - Prof. Pleno
> UESC/DCET/Brasil
> joseclaudio.faria at gmail.com
> ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
>
> Em 27 de julho de 2011 07:36, Ivan Bezerra Allaman
>  escreveu:
>  > Adriano,
> > Se a variável não tem distribuição normal e você não quer utilizar
> nenhuma
> > transformação, terá que recorrer a metodologia dos Modelos Lineares
> > Generalizados (GLM) e encontrar a família de distribuição que se encaixe
> aos
> > seus dados. Eu o aconselho tentar uma transformação primeiro ser você não
> > estiver familiarizado com a metodologia GLM. Já tentou boxcox?
> > Allaman
> > (S,f,P)
> >
> > Ivan Bezerra Allaman
> > Doutor em Zootecnia/UFLA
> > email e msn - ivanala...@yahoo.com.br
> >
> > ___
> > R-br mailing list
> > R-br@listas.c3sl.ufpr.br
> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código
> > mínimo reproduzível.
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

[R-br] Experimento fatorial com tratamento adicional sem normalidade

[Adriano Carvalho Costa]

Boa noite, solicitar a ajuda de voces para realizar uma análise de um
experimento fatorial com um tratamento adicional sem normalidade.

Utilizando a base de dados abaixo (Fonte: Walmes) e considerando que ela não
tem distribuição normal como posso fazer a análise sem utilizar
transformação?


fa <- read.table("http://www.leg.ufpr.br/~walmes/cursoR/fat-adi.txt";,
header=TRUE)
attach(fa)
fa
names(fa)
concentração=concentraÃ.Ã.o
fa <- transform(fa, concentração=factor(concentração), bloc=factor(bloc))
m0 <- aov(media~bloc+trat, fa)
anova(m0)

m1 <- aov(media~bloc+origem*concentração, data=subset(fa, trat!="TEST"))
summary(m1)

Obrigado!
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Re: [R-br] número mínimo de casos em glm e betafit

[Adriano Carvalho Costa]

Edson,

A interpretação do Daniel esta correta?


Em 10 de junho de 2011 10:30, Daniel Marcelino escreveu:

> Obrigado pela dica de livro.
> Se eu entendi o seu exemplo, ele está indicando se eu tiver usando 4
> parâmetros e meu n for 5, minha precisão será de apenas 35%. Esse
> valor sobre para 90% se, por outro lado meu n for 15, correto?
>
> Daniel
>
> 2011/6/10 Edson Lira :
>  > Daaniel, dê uma olhada nestes exemplos do R.
> >
> >
> > power.anova.test(groups=4, n=5, between.var=1, within.var=3)
> > # Power = 0.3535594
> >
> > power.anova.test(groups=4, between.var=1, within.var=3,
> >  power=.80)
> > # n = 11.92613
> >
> > ## Assume we have prior knowledge of the group means:
> > groupmeans <- c(120, 130, 140, 150)
> > power.anova.test(groups = length(groupmeans),
> >  between.var=var(groupmeans),
> >  within.var=500, power=.90) # n = 15.18834
> >
> > No crcpress.com, tem um livro que talvez possa te ajudar muito.
> Probability
> > and Statistics with R.
> >
> >
> >
> >
> > Edson Lira
> > Estatístico
> > Manaus-Amazonas
> > 
> > De: Daniel Marcelino 
> > Para: r-br@listas.c3sl.ufpr.br
> > Enviadas: Quinta-feira, 9 de Junho de 2011 12:52
> > Assunto: Re: [R-br] número mínimo de casos em glm e betafit
> >
> > Acho que não é o caso, o modelo tem 6 parâmetros e 42 e 35 obs. Minha
> > dúvida surgiu porque um colega usou esses modelos, mas sem a certeza
> > do tamanho da amostra.
> >
> > Obrigado,
> > Daniel
> >
> > 2011/6/9 Paulo Justiniano :
> >> Bem Daniel, is vai depender de quantoas parametros(p)  voce estiver
> >> tentando
> >> estimar (quantas covariáveis tem etc).
> >> Voce sempre precisa de mais observacoes (n) do que parametros no modelo
> >> para
> >> poder estimar e podem obter os erros padroes das estimativas.
> >> Afora isto, os algoritmos devem rodar mas evidentemente a incerteza e
> >> mesmo
> >> estabilidade numerica vao ser comprometidas se n ~ p
> >>
> >>
> >> On Thu, 9 Jun 2011, Daniel Marcelino wrote:
> >>
> >>> Caros, alguém saberia me informar se existe algum limite mínimo de
> >>> casos para uma análise de glm, gls e betafit? Não achei nenhuma
> >>> informação que me levasse diretamente a este ponto.
> >>> ___
> >>> R-br mailing list
> >>> R-br@listas.c3sl.ufpr.br
> >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >>
> >> ___
> >> R-br mailing list
> >> R-br@listas.c3sl.ufpr.br
> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >>
> >>
> >
> >
> >
> > --
> > Daniel Marcelino
> > http://danielmarcelino.zip.net
> > Skype: d_marcelino
> > ___
> > R-br mailing list
> > R-br@listas.c3sl.ufpr.br
> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >
> >
> >
> > ___
> > R-br mailing list
> > R-br@listas.c3sl.ufpr.br
> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> >
> >
>
>
>
> --
> Daniel Marcelino
> http://danielmarcelino.zip.net
> Skype: d_marcelino
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Re: [R-br] teste de Homogeneidade e Normalidade

[Adriano Carvalho Costa]

boção, ele respondeu um negocio pra mim me chamando de Cesar e eu nem tinha
perguntado nada!!
rsrsrsrsrs

Em 2 de junho de 2011 10:56, Ivan Bezerra Allaman
escreveu:

>  Adriano,
>
> Para dizer que alguém está doido você precisa ter argumentos mais forte do
> que o dele. E pelo o que eu conheço dos dois, você precisará estudar
> muitoo para ter um argumento melhor do que a do Prof.
> Walmes.
>
> Allaman
> (S,f,P)
>
> *M.Sc Ivan Bezerra Allaman*
> Zootecnista
> Doutorando em Produção Animal/Aquicultura - UFLA
> **email e msn - ivanala...@yahoo.com.br
> Tel: (35)3826-6608/9900-2924
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Re: [R-br] teste de Homogeneidade e Normalidade

[Adriano Carvalho Costa]

Tá doido Walmes!!

Em 2 de junho de 2011 10:31, Walmes Zeviani escreveu:

> Cesar,
>
> As coisas estão meio confusas (pelo menos para mim). Uma coisa é
> distribuição marginal dos dados f(Y), outra coisa é a distribuição
> condicional dos dados como função das regressoras f(Y|x). Quando assumimos
> normalidade, estamos nos referindo a f(Y|x), e pelo fato do erro (E) ser
> aditivo nesses modelos, então f(E) é normal. A f(Y) não vem ao caso. Usando
> o seu exemplo, é necessário que f(consumo|velocidade) tenha distribuição
> normal, mas sobre a f(consumo) não há nenhuma suposição. Muitas pessoas
> pensam que a normalidade deve ser nos dados, e elas aplicam testes de
> normalidade aos dados (Y), mas a normalidade deve estar presente em Y|x.
>
> Walmes.
>
> ==
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: wal...@ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==
>
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br