Enviar submissões para a lista de discussão R-br para
r-br@listas.c3sl.ufpr.br
Para se cadastrar ou descadastrar via WWW, visite o endereço
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou
corpo da mensagem para
r-br-request@listas.c3sl.ufpr.br
Você poderá entrar em contato com a pessoa que gerencia a lista pelo
endereço
r-br-owner@listas.c3sl.ufpr.br
Quando responder, por favor edite sua linha Assunto assim ela será
mais específica que "Re: Contents of R-br digest..."
Tópicos de Hoje:
1. IC 95% (Edmar Caldas)
2. Re: IC 95% (Marcus Nunes)
3. Cálculo do CAP com base na repetição da árvore (Marcelo Laia)
----------------------------------------------------------------------
Message: 1
Date: Sat, 3 Nov 2018 01:21:39 +0000 (UTC)
From: Edmar Caldas <edmar_caldas@yahoo.com.br>
To: a lista Brasileira oficial de discussão do programa R.
<r-br@listas.c3sl.ufpr.br>
Subject: [R-br] IC 95%
Message-ID: <188698043.47848.1541208099198@mail.yahoo.com>
Content-Type: text/plain; charset="utf-8"
Boa Noite!
Pq meus intervalos de confiança deram negativos, se a média é positiva? Outra notação, valor do teste t também deu negativo?
> t.test(Pesquisa$wwwhr~Pesquisa$sex, var.equal=FALSE)
Welch Two Sample t-test
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
Edmar
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4d191f27/attachment-0001.html>
------------------------------
Message: 2
Date: Sat, 3 Nov 2018 07:43:11 -0300
From: Marcus Nunes <marcus.nunes@gmail.com>
To: a lista Brasileira oficial de discussão do programa R.
<r-br@listas.c3sl.ufpr.br>
Subject: Re: [R-br] IC 95%
Message-ID:
<CA+QGQvuJ2brzVLqirK1-thcuEMRf+=Vfda4Nh3-w0ZAfoY=C8w@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Este é um teste t de comparação da diferença entre as médias de dois
grupos. Por padrão, o R assume que os níveis dos grupos são dados em ordem
alfabética. Portanto, no teu caso, Grupo 1 é FEMININO e Grupo 2 é
MASCULINO. Desta forma, tua hipótese nula é
H_0: mu_1 - mu_2 = 0
Ou seja,
H_0: mu_FEMININO - mu_MASCULINO = 0
O estimador será X_barra_FEMININO - X_barra_MASCULINO. Portanto,
4.920046 - 7.106892
= -2.186846. A partir disso, é fácil ver porque o intervalo de confiança é
todo negativo (dica: há diferença entre as médias, com FEMININO <
MASCULINO) e porque a estatística do teste deu negativa também.
--
Marcus Nunes
Professor Adjunto
Universidade Federal do Rio Grande do Norte
Centro de Ciências Exatas e da Terra
Departamento de Estatística
Laboratório de Estatística Aplicada
marcus.nunes@ccet.ufrn.br
https://marcusnunes.me/
On Fri, Nov 2, 2018 at 10:22 PM Edmar Caldas por (R-br) <
r-br@listas.c3sl.ufpr.br> wrote:
> Boa Noite!
>
> Pq meus intervalos de confiança deram negativos, se a média é positiva?
> Outra notação, valor do teste t também deu negativo?
>
> > t.test(Pesquisa$wwwhr~Pesquisa$sex, var.equal=FALSE)
>
> Welch Two Sample t-test
>
> data: Pesquisa$wwwhr by Pesquisa$sex
> t = -4.7853, df = 1314.3, p-value = 1.9e-06
> alternative hypothesis: true difference in means is not equal to 0
> 95 percent confidence interval:
> -3.083359 -1.290332
> sample estimates:
> mean in group FEMININO mean in group MASCULINO
> 4.920046 7.106892
>
>
>
> Edmar
> _______________________________________________
> 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181103/4bc0b32e/attachment-0001.html>
------------------------------
Message: 3
Date: Sat, 3 Nov 2018 10:40:56 -0300
From: Marcelo Laia <marcelolaia@gmail.com>
To: a lista Brasileira oficial de discussão do programa R.
<r-br@listas.c3sl.ufpr.br>
Subject: [R-br] Cálculo do CAP com base na repetição da árvore
Message-ID: <20181103134056.GC2176@localhost>
Content-Type: text/plain; charset=iso-8859-1
Bom dia!
Em Engenharia Florestal é comum termos numa mesma cova duas ou mais hastes
(árvores), a qual chamamos de bifurcadas (para duas), trifurcadas (para três),
etc..
Assim, podemos ter o seguinte:
Bloco Clone Esp Arv CAP HT
1 A 3x3 1 12,5 18
1 A 3x3 2 10,1 11
1 A 3x3 2 11 11,5
1 A 3x3 3 19 21
1 A 3x3 4 19 20
1 A 3x3 5 9 8
1 A 3x3 6 18 17
1 A 3x3 6 17 18
1 A 3x3 6 18,5 17,5
1 A 3x3 7 8 10
(...)
5 C 6x1,5 66 20 21
(...)
6 E 6x1,5 66 18 19
Para prosseguir com os cálculos, é comum reduzir as plantas bifurcadas,
trifurcadas e quadrifurcadas a uma única planta, da seguinte maneira:
CAP = RAIZ(CAP1^2 + CAP2^2 + CAPn^2)
Assim, o exemplo acima ficaria assim:
CAPcalc = RAIZ(10,1^2+11^2) = 14,93
CAPcalc = RAIZ(18^2 + 17^2 + 18,5^2)
Bloco Clone Esp Arv CAP HT CAPcalc
1 A 3x3 1 12,5 18
1 A 3x3 2 10,1 11 14,93
1 A 3x3 2 11 11,5
1 A 3x3 3 19 21
1 A 3x3 4 19 20
1 A 3x3 5 9 8
1 A 3x3 6 18 17 30,91
1 A 3x3 6 17 18
1 A 3x3 6 18,5 17,5
1 A 3x3 7 8 10
(...)
5 C 6x1,5 66 20 21
(...)
6 E 6x1,5 66 18 19
São seis blocos, cinco clones diferentes e seis espaçamentos.
Pergunto: alguém já otimizou isso no R? Tenho 7.473 linhas de dados para
analisar e já fiz isso na unha (Calc Libreoffice) para os 3 anos anteriores.
Como ainda teremos mais 3 anos de medição desse experimento, estou buscando
ajuda. Não tenho nenhum código em mente para enviar. O pedido de ajuda para
este caso e para o outro que irei enviar em seguida parte do zero.
Obrigado
--
Marcelo
------------------------------
Subject: Legenda do Digest
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
------------------------------
Fim da Digest R-br, volume 95, assunto 3
****************************************