<div> </div>
<div> </div>
<div> </div>
<div> </div>
<div><span style="color: #0000ff; font-size: medium;">Por que não estou conseguindo rodar isto?</span></div>
<div> </div>
<div>> dados <- textConnection("Bloco Clone Esp Arv CAP HT</div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 1 12,5 18</span></div>
<div>+ + 1 A 3x3 2 10,1 11</div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 2 11 11,5</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 3 19 21</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 4 19 20</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 5 9 8</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 6 18 17</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 6 17 18</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 6 18,5 17,5</span></div>
<div><span style="font-size: 12.8px;">+ + 1 A 3x3 7 8 10")</span></div>
<div><span style="font-size: 12.8px;">> x <- read.table(dados, header = T, dec = ",")</span></div>
<div><span style="font-size: 12.8px;">Error in read.table(dados, header = T, dec = ",") : </span></div>
<div><span style="font-size: 12.8px;">  duplicate 'row.names' are not allowed</span></div>
<div> </div>
<div><br /><br /></div>
<div><br />> dados <- textConnection("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 />> x <- read.table(dados, header = T, dec = ",");x<br /> Bloco Clone Esp Arv CAP HT<br />1 1 A 3x3 1 12.5 18.0<br />2 1 A 3x3 2 10.1 11.0<br />3 1 A 3x3 2 11.0 11.5<br />4 1 A 3x3 3 19.0 21.0<br />5 1 A 3x3 4 19.0 20.0<br />6 1 A 3x3 5 9.0 8.0<br />7 1 A 3x3 6 18.0 17.0<br />8 1 A 3x3 6 17.0 18.0<br />9 1 A 3x3 6 18.5 17.5<br />10 1 A 3x3 7 8.0 10.0<br />> CAP.sqrt <-function(x){ raiz = sqrt(sum(x^2))<br />+ return(raiz)}<br />> x$COVA <- paste(x$Bloco, x$Clone, x$Esp, x$Arv, sep = "")<br />> attach(x)<br />> CAPcalc <- tapply(CAP, COVA, CAP.sqrt) #note que aqui com o exemplo que<br />você deu só teremos 7 resultados (possíveis combinações entre bloco, clone<br />e espaçamento),<br />> # criando um novo data.frame com a variável CAPcalc<br />> detach(x)<br />> # TESTE<br />> newtable <- data.frame(BLOCO = rep(1, times = 7),<br />+ CLONE = rep("A", times = 7),<br />+ ESP = rep("3x3", times = 7),<br />+ ARV = 1:7,<br />+ CAPcalc)<br />> newtable<br /> BLOCO CLONE ESP ARV CAPcalc<br />1A3x31 1 A 3x3 1 12.50000<br />1A3x32 1 A 3x3 2 14.93352<br />1A3x33 1 A 3x3 3 19.00000<br />1A3x34 1 A 3x3 4 19.00000<br />1A3x35 1 A 3x3 5 9.00000<br />1A3x36 1 A 3x3 6 30.90712<br />1A3x37 1 A 3x3 7 8.00000<br /><br />Em sáb, 3 de nov de 2018 às 12:00, <r-br-request@listas.c3sl.ufpr.br><br />escreveu:<br /><br />> Enviar submissões para a lista de discussão R-br para<br />> r-br@listas.c3sl.ufpr.br<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" 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 />> r-br-request@listas.c3sl.ufpr.br<br />><br />> Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br />> endereço<br />> r-br-owner@listas.c3sl.ufpr.br<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 <edmar_caldas@yahoo.com.br><br />> To: a lista Brasileira oficial de discussão do programa R.<br />> <r-br@listas.c3sl.ufpr.br><br />> Subject: [R-br] IC 95%<br />> Message-ID: <188698043.47848.1541208099198@mail.yahoo.com><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 é<br />> 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 =<br />> 1.9e-06alternative hypothesis: true difference in means is not equal to 095<br />> percent confidence interval: -3.083359 -1.290332sample estimates: mean in<br />> group FEMININO mean in group MASCULINO 4.920046<br />> 7.106892<br />><br />><br />> Edmar<br />> -------------- Próxima Parte ----------<br />> Um anexo em HTML foi limpo...<br />> URL: <<br />> <a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4d191f27/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4d191f27/attachment-0001.html</a><br />> ><br />><br />> ------------------------------<br />><br />> Message: 2<br />> Date: Sat, 3 Nov 2018 07:43:11 -0300<br />> From: Marcus Nunes <marcus.nunes@gmail.com><br />> To: a lista Brasileira oficial de discussão do programa R.<br />> <r-br@listas.c3sl.ufpr.br><br />> Subject: Re: [R-br] IC 95%<br />> Message-ID:<br />> <CA+QGQvuJ2brzVLqirK1-thcuEMRf+=Vfda4Nh3-w0ZAfoY=<br />> C8w@mail.gmail.com><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 />> marcus.nunes@ccet.ufrn.br<br />> <a href="https://marcusnunes.me/" 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 />> r-br@listas.c3sl.ufpr.br> 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 />> > R-br@listas.c3sl.ufpr.br<br />> > <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" 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: <<br />> <a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4bc0b32e/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4bc0b32e/attachment-0001.html</a><br />> ><br />><br />> ------------------------------<br />><br />> Message: 3<br />> Date: Sat, 3 Nov 2018 10:40:56 -0300<br />> From: Marcelo Laia <marcelolaia@gmail.com><br />> To: a lista Brasileira oficial de discussão do programa R.<br />> <r-br@listas.c3sl.ufpr.br><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<br />> 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<br />> 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 />> R-br@listas.c3sl.ufpr.br<br />> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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 />><br /><br /><br />-- <br />*Caio Cézar Guedes Corrêa* <<a href="http://lattes.cnpq.br/3271223159819454" target="_blank">http://lattes.cnpq.br/3271223159819454</a>><br />Universidade Estadual do Norte Fluminense Darcy Ribeiro - UENF<br />Genética e Melhoramento de Plantas (Plant proteomics)<br />_______________________________________________<br />R-br mailing list<br />R-br@listas.c3sl.ufpr.br<br /><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forne�a c�digo m�nimo reproduz�vel.</div>