<div dir="ltr"><div dir="ltr"><div>Marcelo, acredito que isso pode te ajudar...</div><div><br></div><div><br></div><div>> dados <- textConnection("Bloco Clone Esp Arv CAP HT</div><div>+ 1 A 3x3 1 12,5 18</div><div>+ 1 A 3x3 2 10,1 11</div><div>+ 1 A 3x3 2 11 11,5</div><div>+ 1 A 3x3 3 19 21</div><div>+ 1 A 3x3 4 19 20</div><div>+ 1 A 3x3 5 9 8</div><div>+ 1 A 3x3 6 18 17</div><div>+ 1 A 3x3 6 17 18</div><div>+ 1 A 3x3 6 18,5 17,5</div><div>+ 1 A 3x3 7 8 10")</div><div>> x <- read.table(dados, header = T, dec = ",");x</div><div> Bloco Clone Esp Arv CAP HT</div><div>1 1 A 3x3 1 12.5 18.0</div><div>2 1 A 3x3 2 10.1 11.0</div><div>3 1 A 3x3 2 11.0 11.5</div><div>4 1 A 3x3 3 19.0 21.0</div><div>5 1 A 3x3 4 19.0 20.0</div><div>6 1 A 3x3 5 9.0 8.0</div><div>7 1 A 3x3 6 18.0 17.0</div><div>8 1 A 3x3 6 17.0 18.0</div><div>9 1 A 3x3 6 18.5 17.5</div><div>10 1 A 3x3 7 8.0 10.0</div><div>> CAP.sqrt <-function(x){ raiz = sqrt(sum(x^2))</div><div>+ return(raiz)}</div><div>> x$COVA <- paste(x$Bloco, x$Clone, x$Esp, x$Arv, sep = "")</div><div>> attach(x)</div><div>> CAPcalc <- tapply(CAP, COVA, CAP.sqrt) #note que aqui com o exemplo que você deu só teremos 7 resultados (possíveis combinações entre bloco, clone e espaçamento),</div><div>> # criando um novo data.frame com a variável CAPcalc</div><div>> detach(x)</div><div>> # TESTE</div><div>> newtable <- data.frame(BLOCO = rep(1, times = 7),</div><div>+ CLONE = rep("A", times = 7),</div><div>+ ESP = rep("3x3", times = 7),</div><div>+ ARV = 1:7,</div><div>+ CAPcalc)</div><div>> newtable</div><div> BLOCO CLONE ESP ARV CAPcalc</div><div>1A3x31 1 A 3x3 1 12.50000</div><div>1A3x32 1 A 3x3 2 14.93352</div><div>1A3x33 1 A 3x3 3 19.00000</div><div>1A3x34 1 A 3x3 4 19.00000</div><div>1A3x35 1 A 3x3 5 9.00000</div><div>1A3x36 1 A 3x3 6 30.90712</div><div>1A3x37 1 A 3x3 7 8.00000</div></div></div><br><div class="gmail_quote"><div dir="ltr">Em sáb, 3 de nov de 2018 às 12:00, <<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Enviar submissões para a lista de discussão R-br para <br>
<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para <br>
<a href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">r-br-request@listas.c3sl.ufpr.br</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
<a href="mailto:r-br-owner@listas.c3sl.ufpr.br" target="_blank">r-br-owner@listas.c3sl.ufpr.br</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
1. IC 95% (Edmar Caldas)<br>
2. Re: IC 95% (Marcus Nunes)<br>
3. Cálculo do CAP com base na repetição da árvore (Marcelo Laia)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sat, 3 Nov 2018 01:21:39 +0000 (UTC)<br>
From: Edmar Caldas <<a href="mailto:edmar_caldas@yahoo.com.br" target="_blank">edmar_caldas@yahoo.com.br</a>><br>
To: a lista Brasileira oficial de discussão do programa R.<br>
<<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] IC 95%<br>
Message-ID: <<a href="mailto:188698043.47848.1541208099198@mail.yahoo.com" target="_blank">188698043.47848.1541208099198@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Boa Noite!<br>
Pq meus intervalos de confiança deram negativos, se a média é positiva? Outra notação, valor do teste t também deu negativo?<br>
> t.test(Pesquisa$wwwhr~Pesquisa$sex, var.equal=FALSE)<br>
Welch Two Sample t-test<br>
data: Pesquisa$wwwhr by Pesquisa$sext = -4.7853, df = 1314.3, p-value = 1.9e-06alternative hypothesis: true difference in means is not equal to 095 percent confidence interval: -3.083359 -1.290332sample estimates: mean in group FEMININO mean in group MASCULINO 4.920046 7.106892 <br>
<br>
<br>
Edmar<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4d191f27/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4d191f27/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sat, 3 Nov 2018 07:43:11 -0300<br>
From: Marcus Nunes <<a href="mailto:marcus.nunes@gmail.com" target="_blank">marcus.nunes@gmail.com</a>><br>
To: a lista Brasileira oficial de discussão do programa R.<br>
<<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] IC 95%<br>
Message-ID:<br>
<CA+QGQvuJ2brzVLqirK1-thcuEMRf+=Vfda4Nh3-w0ZAfoY=<a href="mailto:C8w@mail.gmail.com" target="_blank">C8w@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Este é um teste t de comparação da diferença entre as médias de dois<br>
grupos. Por padrão, o R assume que os níveis dos grupos são dados em ordem<br>
alfabética. Portanto, no teu caso, Grupo 1 é FEMININO e Grupo 2 é<br>
MASCULINO. Desta forma, tua hipótese nula é<br>
<br>
H_0: mu_1 - mu_2 = 0<br>
<br>
Ou seja,<br>
<br>
H_0: mu_FEMININO - mu_MASCULINO = 0<br>
<br>
O estimador será X_barra_FEMININO - X_barra_MASCULINO. Portanto,<br>
4.920046 - 7.106892<br>
= -2.186846. A partir disso, é fácil ver porque o intervalo de confiança é<br>
todo negativo (dica: há diferença entre as médias, com FEMININO <<br>
MASCULINO) e porque a estatística do teste deu negativa também.<br>
--<br>
Marcus Nunes<br>
Professor Adjunto<br>
Universidade Federal do Rio Grande do Norte<br>
Centro de Ciências Exatas e da Terra<br>
Departamento de Estatística<br>
Laboratório de Estatística Aplicada<br>
<a href="mailto:marcus.nunes@ccet.ufrn.br" target="_blank">marcus.nunes@ccet.ufrn.br</a><br>
<a href="https://marcusnunes.me/" rel="noreferrer" target="_blank">https://marcusnunes.me/</a><br>
<br>
<br>
<br>
On Fri, Nov 2, 2018 at 10:22 PM Edmar Caldas por (R-br) <<br>
<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>> wrote:<br>
<br>
> Boa Noite!<br>
><br>
> Pq meus intervalos de confiança deram negativos, se a média é positiva?<br>
> Outra notação, valor do teste t também deu negativo?<br>
><br>
> > t.test(Pesquisa$wwwhr~Pesquisa$sex, var.equal=FALSE)<br>
><br>
> Welch Two Sample t-test<br>
><br>
> data: Pesquisa$wwwhr by Pesquisa$sex<br>
> t = -4.7853, df = 1314.3, p-value = 1.9e-06<br>
> alternative hypothesis: true difference in means is not equal to 0<br>
> 95 percent confidence interval:<br>
> -3.083359 -1.290332<br>
> sample estimates:<br>
> mean in group FEMININO mean in group MASCULINO<br>
> 4.920046 7.106892<br>
><br>
><br>
><br>
> Edmar<br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4bc0b32e/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4bc0b32e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Sat, 3 Nov 2018 10:40:56 -0300<br>
From: Marcelo Laia <<a href="mailto:marcelolaia@gmail.com" target="_blank">marcelolaia@gmail.com</a>><br>
To: a lista Brasileira oficial de discussão do programa R.<br>
<<a href="mailto:r-br@listas.c3sl.ufpr.br" target="_blank">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] Cálculo do CAP com base na repetição da árvore<br>
Message-ID: <20181103134056.GC2176@localhost><br>
Content-Type: text/plain; charset=iso-8859-1<br>
<br>
Bom dia!<br>
<br>
Em Engenharia Florestal é comum termos numa mesma cova duas ou mais hastes<br>
(árvores), a qual chamamos de bifurcadas (para duas), trifurcadas (para três),<br>
etc..<br>
<br>
Assim, podemos ter o seguinte:<br>
<br>
Bloco Clone Esp Arv CAP HT<br>
1 A 3x3 1 12,5 18<br>
1 A 3x3 2 10,1 11<br>
1 A 3x3 2 11 11,5<br>
1 A 3x3 3 19 21<br>
1 A 3x3 4 19 20<br>
1 A 3x3 5 9 8<br>
1 A 3x3 6 18 17<br>
1 A 3x3 6 17 18<br>
1 A 3x3 6 18,5 17,5<br>
1 A 3x3 7 8 10<br>
(...)<br>
5 C 6x1,5 66 20 21<br>
(...)<br>
6 E 6x1,5 66 18 19<br>
<br>
<br>
Para prosseguir com os cálculos, é comum reduzir as plantas bifurcadas,<br>
trifurcadas e quadrifurcadas a uma única planta, da seguinte maneira:<br>
<br>
CAP = RAIZ(CAP1^2 + CAP2^2 + CAPn^2)<br>
<br>
Assim, o exemplo acima ficaria assim:<br>
<br>
CAPcalc = RAIZ(10,1^2+11^2) = 14,93<br>
<br>
CAPcalc = RAIZ(18^2 + 17^2 + 18,5^2)<br>
<br>
Bloco Clone Esp Arv CAP HT CAPcalc<br>
1 A 3x3 1 12,5 18<br>
1 A 3x3 2 10,1 11 14,93<br>
1 A 3x3 2 11 11,5<br>
1 A 3x3 3 19 21<br>
1 A 3x3 4 19 20<br>
1 A 3x3 5 9 8<br>
1 A 3x3 6 18 17 30,91<br>
1 A 3x3 6 17 18<br>
1 A 3x3 6 18,5 17,5<br>
1 A 3x3 7 8 10<br>
(...)<br>
5 C 6x1,5 66 20 21<br>
(...)<br>
6 E 6x1,5 66 18 19<br>
<br>
São seis blocos, cinco clones diferentes e seis espaçamentos.<br>
<br>
Pergunto: alguém já otimizou isso no R? Tenho 7.473 linhas de dados para<br>
analisar e já fiz isso na unha (Calc Libreoffice) para os 3 anos anteriores.<br>
Como ainda teremos mais 3 anos de medição desse experimento, estou buscando<br>
ajuda. Não tenho nenhum código em mente para enviar. O pedido de ajuda para<br>
este caso e para o outro que irei enviar em seguida parte do zero.<br>
<br>
Obrigado<br>
<br>
-- <br>
Marcelo<br>
<br>
<br>
------------------------------<br>
<br>
Subject: Legenda do Digest<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
------------------------------<br>
<br>
Fim da Digest R-br, volume 95, assunto 3<br>
****************************************<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="text-align:left"><a href="http://lattes.cnpq.br/3271223159819454" style="font-family:georgia,serif" target="_blank"><b>Caio Cézar Guedes Corrêa</b></a><br></div><div style="text-align:left"><font face="georgia, serif"><span style="color:rgb(68,68,68)">Universidade Estadual do Norte Fluminense Darcy Ribeiro - UENF</span><br></font></div><div style="text-align:left"><span style="color:rgb(68,68,68)"><font face="georgia, serif">Genética e Melhoramento de Plantas (Plant proteomics)</font></span></div></div></div></div></div></div></div>