<div><div>Oi Walmes, </div><div><br></div><div>São 10 blocos, 5 deles em remansos e outros cinco em corredeiras.</div><div><br></div><div>Em cada bloco existem dois plots colocados lado a lado, sendo um com exclusão elétrica e outro sem exclusão (ou seja, controle). Dentro de cada um desses dois plots existem substratos de acrílico lisos e rugosos.</div>
<div><br></div><div>Preciso avaliar o efeito do habitat (remansos x corredeiras), o efeito da exclusão (exclusão x controle) e o efeito do tipo de substrato (liso x rugoso) levando em conta o design em spli-plot.</div><div>
<br></div><div>Ficou mais claro? </div><div><br></div><div>Obrigada,</div><div>Fabiana</div></div><div><br></div><div><br></div><div><br></div><div><br></div><br><br><div class="gmail_quote">Em 13 de dezembro de 2011 17:38, <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a>></span> escreveu:<br>
<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">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" 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">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">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. split-plot (Fabiana Schneck)<br>
2. Re: Transformar Variáveis discretas em Numeros Binários -><br>
Pre-processamento para Redes Neurais (Elias T. Krainski)<br>
3. Rotina order (Edson Lira)<br>
4. Re: Rotina order (Elias T. Krainski)<br>
5. Re: MQM (Elias T. Krainski)<br>
6. Re: Modificar valores de dataframe com base em outro<br>
dataframe (Elias T. Krainski)<br>
7. Categorização de variáveis (Edson Lira)<br>
8. Re: Categorização de variáveis (Luis Iván Ortiz Valencia)<br>
9. Re: Categorização de variáveis (Cesar Rabak)<br>
10. Re: Categorização de variáveis (Elias T. Krainski)<br>
11. Re: Categorização de variáveis (Edson Lira)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 13 Dec 2011 12:25:13 -0200<br>
From: Fabiana Schneck <<a href="mailto:fabiana.schneck@gmail.com">fabiana.schneck@gmail.com</a>><br>
To: r-br <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] split-plot<br>
Message-ID:<br>
<<a href="mailto:CALjK-AcaG3AELQAOopMDZ1X7oFYOzNg4957pn3-3BG9J%2Bt07oA@mail.gmail.com">CALjK-AcaG3AELQAOopMDZ1X7oFYOzNg4957pn3-3BG9J+t07oA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="windows-1252"<br>
<br>
Prezados colegas,<br>
<br>
tenho um experimento com o seguinte design split-plot:<br>
<br>
fator 'ambiente': 5 corredeiras e 5 remansos<br>
em cada plot de 'ambiente' (num total de 10) há outro fator: 'peixes' com<br>
dois tratamentos: exclusão e controle<br>
ainda, dentro de cada sub-plot do fator 'peixes' existe o fator 'substrato'<br>
com dois tratamentos: liso e rugoso.<br>
<br>
A variável resposta é quantidade de clorofila em cada unidade experimental<br>
(cada substrato).<br>
Um exemplo dos dados:<br>
<br>
<br>
<br>
clorofila<br>
<br>
ambiente<br>
<br>
peixes<br>
<br>
substrato<br>
<br>
bloco<br>
<br>
0.114<br>
<br>
1<br>
<br>
1<br>
<br>
1<br>
<br>
1<br>
<br>
0.160<br>
<br>
1<br>
<br>
1<br>
<br>
2<br>
<br>
1<br>
<br>
0.126<br>
<br>
1<br>
<br>
2<br>
<br>
1<br>
<br>
1<br>
<br>
0.149<br>
<br>
1<br>
<br>
2<br>
<br>
2<br>
<br>
1<br>
<br>
0.068<br>
<br>
1<br>
<br>
1<br>
<br>
1<br>
<br>
2<br>
<br>
0.137<br>
<br>
1<br>
<br>
1<br>
<br>
2<br>
<br>
2<br>
<br>
0.080<br>
<br>
1<br>
<br>
2<br>
<br>
1<br>
<br>
2<br>
<br>
0.160<br>
<br>
1<br>
<br>
2<br>
<br>
2<br>
<br>
2<br>
....<br>
....<br>
<br>
0.080<br>
<br>
2<br>
<br>
1<br>
<br>
1<br>
<br>
10<br>
<br>
0.068<br>
<br>
2<br>
<br>
1<br>
<br>
2<br>
<br>
10<br>
<br>
0.114<br>
<br>
2<br>
<br>
2<br>
<br>
1<br>
<br>
10<br>
<br>
0.124<br>
<br>
2<br>
<br>
2<br>
<br>
2<br>
<br>
10<br>
<br>
Considerei que os fatores 'peixes' e 'substrato' estão em blocos (10),<br>
enquanto o fator 'ambiente' está acima dos blocos (5 blocos são no ambiente<br>
corredeiras e 5 no ambiente remansos).<br>
<br>
Analise os dados utilizando dois modelos distintos:<br>
<br>
*aov (log(cloro)~(amb+peixes+substrato)^2+Error(bloco/peixes))*<br>
*<br>
*<br>
<br>
Error: bloco<br>
<br>
Df Sum Sq Mean Sq F value Pr(>F)<br>
<br>
amb 1 0.1162 0.11618 0.128 0.7297<br>
<br>
Residuals 8 7.2599 0.90749<br>
<br>
<br>
<br>
Error: bloco:peixes<br>
<br>
Df Sum Sq Mean Sq F value Pr(>F)<br>
<br>
peixes 1 2.96190 2.96190 14.6147 0.005069 **<br>
<br>
amb:peixes 1 0.30121 0.30121 1.4862 0.257528<br>
<br>
Residuals 8 1.62133 0.20267<br>
<br>
---<br>
<br>
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1<br>
<br>
<br>
<br>
Error: Within<br>
<br>
Df Sum Sq Mean Sq F value Pr(>F)<br>
<br>
substrato 1 1.9924 1.99244 10.2620 0.005209 **<br>
<br>
amb:substrato 1 0.0013 0.00134 0.0069 0.934840<br>
<br>
peixes:substrato 1 0.5934 0.59335 3.0560 0.098473 .<br>
<br>
Residuals 17 3.3007 0.19416<br>
<br>
<br>
*<br>
*<br>
*lme(log(cloro)~(amb+peixes+**substrato)^2, random=~1|bloco/peixes))*<br>
<br>
Linear mixed-effects model fit by REML<br>
Data: NULL<br>
AIC BIC logLik<br>
85.58484 100.5499 -32.79242<br>
<br>
Random effects:<br>
Formula: ~1 | bloco<br>
(Intercept)<br>
StdDev: 0.4197679<br>
<br>
Formula: ~1 | peixes %in% bloco<br>
(Intercept) Residual<br>
StdDev: 0.06522495 0.4406333<br>
<br>
Fixed effects: log(cloro) ~ (amb + peixes + substrato)^2<br>
Value Std.Error DF t-value p-value<br>
(Intercept) -2.4342471 0.2647064 17 -9.196027 0.0000<br>
amb2 -0.2929029 0.3611524 8 -0.811023 0.4408<br>
peixes2 0.6142687 0.2448449 8 2.508808 0.0364<br>
substrato2 0.6783940 0.2413448 17 2.810891 0.0120<br>
amb2:peixes2 0.3471056 0.2847218 8 1.219104 0.2575<br>
amb2:substrato2 0.0231237 0.2786810 17 0.082976 0.9348<br>
peixes2:substrato2 -0.4871762 0.2786810 17 -1.748150 0.0985<br>
Correlation:<br>
(Intr) amb2 peixs2 sbstr2 amb2:p2 amb2:s2<br>
amb2 -0.682<br>
peixes2 -0.462 0.229<br>
substrato2 -0.456 0.223 0.329<br>
amb2:peixes2 0.269 -0.394 -0.581 0.000<br>
amb2:substrato2 0.263 -0.386 0.000 -0.577 0.000<br>
peixes2:substrato2 0.263 0.000 -0.569 -0.577 0.000 0.000<br>
<br>
Standardized Within-Group Residuals:<br>
Min Q1 Med Q3 Max<br>
-2.42000225 -0.41789060 -0.02451742 0.44159262 1.87090777<br>
<br>
Number of Observations: 40<br>
Number of Groups:<br>
bloco peixes %in% bloco<br>
10 20<br>
<br>
<br>
<br>
Tenho dúvidas se estou definindo de forma correta os modelos e também se<br>
devo usar aov ou lme, uma vez que os resultados dos termos isolados diferem<br>
entre as duas funções (porém os resultados são indênticos para as<br>
interações).<br>
<br>
Desde já obrigada pela atenção,<br>
<br>
Fabiana<br>
<br>
--<br>
-----<br>
Fabiana Schneck<br>
Programa de Pós-Graduação em Ecologia<br>
Universidade Federal do Rio Grande do Sul<br>
CP 15007 - CEP 91501-970<br>
Porto Alegre, RS, Brasil<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/6ca4ee6e/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/6ca4ee6e/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 13 Dec 2011 06:30:03 -0800 (PST)<br>
From: "Elias T. Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Transformar Variáveis discretas em Numeros<br>
Binários -> Pre-processamento para Redes Neurais<br>
Message-ID:<br>
<<a href="mailto:1323786603.22667.YahooMailNeo@web114716.mail.gq1.yahoo.com">1323786603.22667.YahooMailNeo@web114716.mail.gq1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Fernando, talvez eu não tenha entendido exatamente o que você quer. Veja dois exemplos, considerando as colunas 2 e 3 de 'warpbreaks' (disponivel no 'datasets', como o 'iris').<br>
<br>
<br>
### Considerando que um 'factor' é um inteiro, apenas retiramos os 'levels', e vc pode manipular<br>
### Se ha 3 ou mais classes, o resultado nao é binário, é inteiro de 1 ao número de classes<br>
sapply(warpbreaks[,2:3], unclass)<br>
sapply(warpbreaks[,2:3], unclass)-1<br>
<br>
<br>
### converte em colunas binárias, como dummies em qualquer modelo regressão<br>
<br>
### se há mais de duas classes, é criado mais de uma coluna<br>
model.matrix(~wool + tension, warpbreaks)[,-1]<br>
<br>
<br>
Att.<br>
<br>
Elias T. Krainski<br>
<br>
<br>
>________________________________<br>
> De: Fernando Neto <<a href="mailto:fernandoneto7@gmail.com">fernandoneto7@gmail.com</a>><br>
>Para: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
>Enviadas: Segunda-feira, 12 de Dezembro de 2011 21:32<br>
>Assunto: [R-br] Transformar Variáveis discretas em Numeros Binários -> Pre-processamento para Redes Neurais<br>
><br>
><br>
>Boa noite, R-anos!<br>
>Por favor, estou precisando ajustar os arquivos de uma base de dados para processamento de uma RNA.<br>
>Estou usando a função decodeClassLabels(...) mas a operacao tá ficando muito custosa.<br>
>Por exemplo<br>
>__<br>
>Masculino Programador<br>
>Masculino Médico<br>
>Feminino Programador<br>
>__<br>
>tem que ser<br>
>__<br>
>0 1 0 1<br>
>1 0 1 0<br>
>__<br>
><br>
><br>
>Só que eu trato cada variável de cada vez. E acho que a operacao de concatenacao que eu uso:<br>
>cbind(...) tá tornando a operacao custosa (além de que a memória passa de 2GB).<br>
><br>
><br>
>PS.: eu faço o seguinte:<br>
><br>
><br>
>for (i in (VETOR COM AS COLUNAS A SEREM TRANSFORMADAS)){<br>
> colAux = decodeClassLabels(vector[,i])<br>
> vector = vector [,-i]<br>
> vector = cbind(vector, colAux) <br>
>} <br>
><br>
><br>
>Alguém sugere algo?<br>
><br>
><br>
>Muito grato!!!<br>
>Fernando<br>
><br>
><br>
>--<br>
>----------------------------------------------------------------<br>
><br>
>Fernando Neto<br>
>Twitter: @fernandompneto<br>
>Facebook: <a href="http://facebook.com/fernandompneto" target="_blank">facebook.com/fernandompneto</a><br>
><br>
>Tecnologia de Ponte<br>
><a href="http://tecnologiadeponte.blogspot.com" target="_blank">http://tecnologiadeponte.blogspot.com</a><br>
>----------------------------------------------------------------<br>
>fmpn2 @ CIn - UFPE<br>
><a href="http://cin.ufpe.br/~fmpn2" target="_blank">http://cin.ufpe.br/~fmpn2</a><br>
><br>
>- Engenharia da Computação - Turma 2009.2 - CIn, UFPE.<br>
>- Monitor de Estatistica e Probabilidade Para Engenharia da Computacao<br>
><br>
><br>
>----------------------------------------------------------------<br>
>Confidencialidade*: *A informação contida nesta mensagem de e-mail,<br>
>incluindo quaisquer anexos, é confidencial e está reservada apenas à pessoa ou entidade para a qual foi endereçada. Se você não é o destinatário ou a pessoa responsável por encaminhar esta mensagem ao destinatário, você está, por meio desta, notificado que não deverá rever, retransmitir, imprimir, copiar, usar ou distribuir esta mensagem de e-mail ou quaisquer anexos. Caso você tenha recebido esta mensagem por engano, por favor, contate o remetente imediatamente e apague esta mensagem de seu computador ou de qualquer outro banco de dados. Muito obrigado.<br>
>Confidentiality Notice: The information contained in this email message, including any attachment, is confidential and is intended only for the person or entity to which it is addressed. If you are neither the intended recipient nor the employee or agent responsible for delivering this message to the intended recipient, you are hereby notified that you may not review, retransmit, convert to hard copy, copy, use or distribute this email message or any attachments to it. If you have received this email in error, please contact the sender immediately and delete this message from any computer or other data bank. Thank you.<br>
><br>
>--------<br>
><br>
>_______________________________________________<br>
>R-br mailing list<br>
><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
><br>
><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/19fb4b33/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/19fb4b33/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 13 Dec 2011 07:22:54 -0800 (PST)<br>
From: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
To: R-br Lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] Rotina order<br>
Message-ID:<br>
<<a href="mailto:1323789774.75434.YahooMailNeo@web160704.mail.bf1.yahoo.com">1323789774.75434.YahooMailNeo@web160704.mail.bf1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Bom dia a todos!<br>
<br>
Com a rotina abaixo eu consigo ordenar a tabela.<br>
<br>
inap<-table(tip_doa$INAPTIDAO,tip_doa$mês)<br>
inap2 = inap[order(rowSums(inap), decreasing=TRUE),]<br>
ftable(addmargins(inap2, c(1,2),<br>
FUN=list(Total = sum)))<br>
<br>
Que me dá a saida abaixo(as primeiras linhas)<br>
<br>
jan fev mar abr mai jun jul ago set out nov Total<br>
<br>
ESTEVE EM ZONA DE MALARIA 465 443 387 362 343 347 343 385 550 397 409 4431<br>
HEMATOCRITO IRREGULAR 299 387 478 473 576 393 266 307 362 381 333 4255<br>
VARIAÇÃO DE PARCEIRO SEM PRESERVATIVO (MAIS DE 4 EM 12 MESES) 249 329 240 152 100 183 150 86 110 124 106 1829<br>
VARIAÇÃO DE PARCEIROS SEM PRESERVATIVO (MAIS DE 1 EM 6 MESES) 190 206 147 83 73 104 87 58 103 56 102 1209<br>
VARIAÇÃO DE PARCEIROS SEXUAIS SEM PRESERVATIVO 0 0 20 106 238 135 109 138 134 142 160 1182<br>
GRIPE 119 112 125 98 105 112 75 96 90 84 102 1118<br>
<br>
Gostaria de acrescentar o ano e ordenar. Quando acrescento o ano me dá o seguinte problema.<br>
<br>
> inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês)<br>
> inap2 = inap[order(rowSums(inap),decreasing=TRUE),]<br>
Erro em inap[order(rowSums(inap), decreasing = TRUE), ] :<br>
número incorreto de dimensões<br>
> inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês)<br>
> inap2 = inap[order(rowSums(inap),decreasing=TRUE),]<br>
Erro em inap[order(rowSums(inap), decreasing = TRUE), ] :<br>
número incorreto de dimensões<br>
> ftable(addmargins(inap2, c(1,2,3),<br>
+ FUN=list(Total = sum)))<br>
Erro em newdimnames[[margin]] : índice fora de limites<br>
<br>
Qual o problema? Alguém pode me ajudar?<br>
<br>
[]'s.<br>
<br>
<br>
Edson Lira<br>
Estatístico<br>
Manaus-Amazonas<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/d6d0f94c/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/d6d0f94c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Tue, 13 Dec 2011 07:43:37 -0800 (PST)<br>
From: "Elias T. Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>>, Edson Lira<br>
<<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
Subject: Re: [R-br] Rotina order<br>
Message-ID:<br>
<<a href="mailto:1323791017.65689.YahooMailNeo@web114705.mail.gq1.yahoo.com">1323791017.65689.YahooMailNeo@web114705.mail.gq1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Com<br>
<br>
inap <- table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês)<br>
'inap' deixa de ter duas dimensões. Portanto, rowSums() é inadequado.<br>
<br>
<br>
Elias T. Krainski<br>
<br>
<br>
>________________________________<br>
> De: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
>Para: R-br Lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
>Enviadas: Terça-feira, 13 de Dezembro de 2011 13:22<br>
>Assunto: [R-br] Rotina order<br>
><br>
><br>
>Bom dia a todos!<br>
><br>
>Com a rotina abaixo eu consigo ordenar a tabela.<br>
><br>
>inap<-table(tip_doa$INAPTIDAO,tip_doa$mês)<br>
>inap2 = inap[order(rowSums(inap), decreasing=TRUE),]<br>
>ftable(addmargins(inap2, c(1,2),<br>
> FUN=list(Total = sum)))<br>
><br>
>Que me dá a saida abaixo(as primeiras linhas)<br>
><br>
> jan fev mar abr mai jun jul ago set out nov Total<br>
> <br>
>ESTEVE EM ZONA DE MALARIA 465 443 387 362 343 347 343 385 550 397 409 4431<br>
>HEMATOCRITO IRREGULAR 299 387 478 473 576 393 266 307 362 381 333 4255<br>
>VARIAÇÃO DE PARCEIRO SEM PRESERVATIVO (MAIS DE 4 EM 12 MESES) 249 329 240 152 100 183 150 86 110 124 106 1829<br>
>VARIAÇÃO DE PARCEIROS SEM PRESERVATIVO (MAIS DE 1 EM 6 MESES) 190 206 147 83 73 104 87 58 103 56 102 1209<br>
>VARIAÇÃO DE PARCEIROS SEXUAIS SEM PRESERVATIVO 0 0 20 106 238 135 109 138 134 142 160 1182<br>
>GRIPE 119 112 125 98 105 112 75 96 90 84 102 1118<br>
><br>
>Gostaria de acrescentar o ano e ordenar. Quando acrescento o ano me dá o seguinte problema.<br>
><br>
>><br>
inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês)<br>
>> inap2 = inap[order(rowSums(inap),decreasing=TRUE),]<br>
>Erro em inap[order(rowSums(inap), decreasing = TRUE), ] :<br>
> número incorreto de dimensões<br>
>> inap<-table(tip_doa$INAPTIDAO,tip_doa$ano,tip_doa$mês)<br>
>> inap2 = inap[order(rowSums(inap),decreasing=TRUE),]<br>
>Erro em inap[order(rowSums(inap), decreasing = TRUE), ] :<br>
> número incorreto de dimensões<br>
>> ftable(addmargins(inap2, c(1,2,3),<br>
>+ FUN=list(Total = sum)))<br>
>Erro em newdimnames[[margin]] : índice fora de limites<br>
><br>
>Qual o problema? Alguém pode me ajudar?<br>
><br>
> []'s.<br>
><br>
><br>
>Edson Lira<br>
>Estatístico<br>
>Manaus-Amazonas<br>
>_______________________________________________<br>
>R-br mailing list<br>
><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
><br>
><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/53dc517b/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/53dc517b/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Tue, 13 Dec 2011 08:00:29 -0800 (PST)<br>
From: "Elias T. Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] MQM<br>
Message-ID:<br>
<<a href="mailto:1323792029.83765.YahooMailNeo@web114716.mail.gq1.yahoo.com">1323792029.83765.YahooMailNeo@web114716.mail.gq1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Se em<br>
<br>
A b = x<br>
A tem dimensão 3x3, então vc quer encontrar 6 valore. Para isso, vc tem 3 valores em b e 3 valores em x. Isso é possível? No usual, quando vc quer saber os 3 valores de b e conhece A e x, vc tem bem mais valores conhecidos que desconhecidos.<br>
<br>
<br>
### Exemplo<br>
<br>
### considerando o contexto de LM<br>
n <- 5<br>
X <- cbind(1, runif(5))<br>
b <- c(10, 5)<br>
y <- rnorm(n, X%*%b, 0.5)<br>
<br>
<br>
### estima b<br>
(XX <- crossprod(X))<br>
(Xy <- crossprod(X, y))<br>
(bhat <- solve(XX, Xy))<br>
<br>
### com bhat e Xy podemos obter XX?<br>
### considere que<br>
### XX %*% bhat = Xy<br>
### para compatibilidade de operações seja<br>
### XX %*% (bhat %*% bhat') = Xy %*% bhat'<br>
<br>
require(Matrix)<br>
bb <- forceSymmetric(tcrossprod(bhat))<br>
solve(bb, tcrossprod(Xy, bhat)) ## nao é igual a XX<br>
<br>
pode ser que este 'sistema' possua mais de uma solução...<br>
<br>
<br>
Att.<br>
<br>
Elias T. Krainski<br>
<br>
<br>
>________________________________<br>
> De: FHRB Toledo <<a href="mailto:fernandohtoledo@gmail.com">fernandohtoledo@gmail.com</a>><br>
>Para: R-Br <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
>Enviadas: Domingo, 11 de Dezembro de 2011 12:07<br>
>Assunto: [R-br] MQM<br>
><br>
><br>
>Senhores,<br>
><br>
>Tenho uma dúvida que não necessáriamente se relaciona com o programa, mas que em muito tem a ver com teoria estatística e algebra de matrizes, espero que aqui ache pelo menos uma sugestão para a resposta ou uma recomendação de literatura.<br>
><br>
>O problema é dado o sistema de equações na forma matricial: Xb = y, pelo MQM, dada a matriz X e o vetor y, obtenho o vetor de soluções por X^-1y, respeitando todas as pressuposições.<br>
><br>
>Se caso conheço ao invés de X e y, conheço os vetores b e y, existe alguma forma de conhecer a matriz X, podendo apenas assegurar que ela é uma matriz quadrada simétrica e positiva definida?<br>
><br>
>att,<br>
>FH<br>
><br>
>_______________________________________________<br>
>R-br mailing list<br>
><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
><br>
><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/018859e7/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/018859e7/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 6<br>
Date: Tue, 13 Dec 2011 08:06:38 -0800 (PST)<br>
From: "Elias T. Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Modificar valores de dataframe com base em outro<br>
dataframe<br>
Message-ID:<br>
<<a href="mailto:1323792398.50539.YahooMailNeo@web114717.mail.gq1.yahoo.com">1323792398.50539.YahooMailNeo@web114717.mail.gq1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
d1 <- data.frame(snp1=c("A","A","A","A","T"), snp2=c("A","G","G","A","G"),<br>
<br>
snp3=c("C","T","C","C","T"), snp4=c("T","G","G","T","T"))<br>
<br>
<br>
d2 <- data.frame(snp=paste("snp",1:4,sep=""), anc=c("A","G","C","T"))<br>
<br>
sapply(1:ncol(d1), function(i) 1-(as.character(d1[,i])==as.character(d2[i,2])))<br>
<br>
<br>
Elias T. Krainski<br>
<br>
<br>
>________________________________<br>
> De: João José de Simoni Gouveia <<a href="mailto:joao.gouveia@univasf.edu.br">joao.gouveia@univasf.edu.br</a>><br>
>Para: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
>Enviadas: Sábado, 10 de Dezembro de 2011 22:21<br>
>Assunto: [R-br] Modificar valores de dataframe com base em outro dataframe<br>
><br>
>Prezados,<br>
><br>
>Possuo dois dataframes (nomeados aqui de 1 e 2). Preciso modificar os valores do dataframe1 com base no dataframe2.<br>
>Mais precisamente:<br>
>Se os valores da coluna snp1 (dataframe1) forem iguais ao valor da linha snp1 (dataframe2), preciso converter estes valores para 0 (zero), caso contrário o valor deve ser convertido para 1 (um).<br>
>Alguém tem alguma dica?<br>
>Att<br>
>JJ<br>
>Abaixo uma ilustração do que eu tenho<br>
>DataFrame 1 --><br>
> snp1 snp2 snp3 snp4<br>
>1 A A C T<br>
>2 A G T G<br>
>3 A G C G<br>
>4 A A C T<br>
>5 T G T T<br>
><br>
>DataFrame 2 --><br>
> snp anc<br>
>1 snp1 A<br>
>2 snp2 G<br>
>3 snp3 C<br>
>4 snp4 T<br>
><br>
>E o que eu preciso ter:<br>
>Dataframe final --><br>
> snp1 snp2 snp3 snp4<br>
>1 0 1 0 0<br>
>2 0 0 1 1<br>
>3 0 0 0 1<br>
>4 0 1 0 0<br>
>5 1 0 1 0<br>
><br>
>--<br>
>============================================<br>
>João José de Simoni Gouveia <br>
>Médico Veterinário, Msc.<br>
>Doutorando em Zootecnia - PDIZ (UFC/UFPB/UFRPE) <br>
>Professor Assistente I - Colegiado Acadêmico de Zootecnia <br>
>Universidade Federal do Vale do São Francisco (UNIVASF) <br>
>----------------------------------------------------------------------------------------<br>
>Endereço: UNIVASF (Campus Ciências Agrárias) <br>
>Rodovia BR 407, Km 12, Lote 543 <br>
>Projeto de Irrigação Senador Nilo Coelho, s/nº-C1 <br>
>Petrolina, PE - CEP 56.300-000 <br>
>Telefone: (87)3986-3800/3804-3801<br>
>Skype: joao_jose_de_simoni_gouveia<br>
>E-mail alternativo: <a href="mailto:jjsgouveia@gmail.com">jjsgouveia@gmail.com</a> <br>
>============================================<br>
>_______________________________________________<br>
>R-br mailing list<br>
><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
><br>
><br>
><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/7ca47fe5/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/7ca47fe5/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 7<br>
Date: Tue, 13 Dec 2011 09:54:57 -0800 (PST)<br>
From: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
To: R-br Lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] Categorização de variáveis<br>
Message-ID:<br>
<<a href="mailto:1323798897.47220.YahooMailNeo@web160703.mail.bf1.yahoo.com">1323798897.47220.YahooMailNeo@web160703.mail.bf1.yahoo.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Boa tarde, caros r-istas, estou com alguns probleminhas na categorização de variáveis.<br>
<br>
Tenho este banco.<br>
pad pas sexo pg<br>
155 134 83 feminino 15<br>
156 72 99 feminino 30<br>
157 82 151 feminino 10<br>
158 23 70 feminino 20<br>
159 75 128 masculino 13<br>
160 73 125 feminino 05<br>
161 125 80 masculino 08<br>
162 65 110 feminino 10<br>
163 132 92 masculino 13<br>
164 114 68 masculino 11<br>
165 83 153 feminino 26<br>
166 166 107 masculino 22<br>
167 114 76 feminino 03<br>
<br>
<br>
<br>
Preciso classificar da seguinte forma:<br>
se pad<85 e pas<130 é normal<br>
se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
e assim por diante.<br>
<br>
<br>
<br>
a outra classificação é a seguinte:<br>
<br>
Para homens o percentual de gordura<br>
<br>
entre 15 e 20 é normal<br>
<br>
entre 20 e 25 levemente obeso<br>
entre 25 e 30 excessivamente obeso<br>
entre 10 e 14 gordura ideal<br>
3 gordura essencial<br>
<br>
<br>
<br>
<br>
<br>
Para homens o percentual de gordura<br>
<br>
entre 25 e 30 é normal<br>
<br>
entre 30 e 35 levemente obeso<br>
entre 35 e 40 excessivamente obeso<br>
entre 14 e 18 gordura ideal<br>
12 gordura essencial<br>
<br>
<br>
<br>
<br>
Como fazer estaS categorizações?<br>
<br>
[]'s.<br>
<br>
Edson Lira<br>
Estatístico<br>
Manaus-Amazonas<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/ddf3f837/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/ddf3f837/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 8<br>
Date: Tue, 13 Dec 2011 15:58:43 -0200<br>
From: Luis Iván Ortiz Valencia <<a href="mailto:liov2067@gmail.com">liov2067@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>, Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
Subject: Re: [R-br] Categorização de variáveis<br>
Message-ID:<br>
<<a href="mailto:CAJ%2BS-MHuH4RMQ35u092HTUwcCfMP3nqkR%2BWdzyouauFD9Yhhng@mail.gmail.com">CAJ+S-MHuH4RMQ35u092HTUwcCfMP3nqkR+WdzyouauFD9Yhhng@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
cut()<br>
<br>
Em 13 de dezembro de 2011 15:54, Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>>escreveu:<br>
<br>
> Boa tarde, caros r-istas, estou com alguns probleminhas na categorização<br>
> de variáveis.<br>
><br>
> Tenho este banco.<br>
> pad pas sexo pg<br>
> 155 134 83 feminino 15<br>
> 156 72 99 feminino 30<br>
> 157 82 151 feminino 10<br>
> 158 23 70 feminino 20<br>
> 159 75 128 masculino 13<br>
> 160 73 125 feminino 05<br>
> 161 125 80 masculino 08<br>
> 162 65 110 feminino 10<br>
> 163 132 92 masculino 13<br>
> 164 114 68 masculino 11<br>
> 165 83 153 feminino 26<br>
> 166 166 107 masculino 22<br>
> 167 114 76 feminino 03<br>
><br>
><br>
> Preciso classificar da seguinte forma:<br>
> se pad<85 e pas<130 é normal<br>
> se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
> e assim por diante.<br>
><br>
><br>
> a outra classificação é a seguinte:<br>
><br>
> Para homens o percentual de gordura<br>
><br>
> entre 15 e 20 é normal<br>
> entre 20 e 25 levemente obeso<br>
> entre 25 e 30 excessivamente obeso<br>
> entre 10 e 14 gordura ideal<br>
> 3 gordura essencial<br>
><br>
><br>
><br>
> Para homens o percentual de gordura<br>
><br>
> entre 25 e 30 é normal<br>
> entre 30 e 35 levemente obeso<br>
> entre 35 e 40 excessivamente obeso<br>
> entre 14 e 18 gordura ideal<br>
> 12 gordura essencial<br>
><br>
><br>
> Como fazer estaS categorizações?<br>
><br>
> []'s.<br>
><br>
> Edson Lira<br>
> Estatístico<br>
> Manaus-Amazonas<br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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>
><br>
<br>
<br>
<br>
--<br>
Luis Iván Ortiz Valencia<br>
Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ<br>
Estatístico Msc.<br>
Spatial Analyst Msc.<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/2702a502/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20111213/2702a502/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 9<br>
Date: Tue, 13 Dec 2011 15:14:02 -0300<br>
From: Cesar Rabak <<a href="mailto:cesar.rabak@gmail.com">cesar.rabak@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Categorização de variáveis<br>
Message-ID: <<a href="mailto:4EE795EA.6090407@gmail.com">4EE795EA.6090407@gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Edson,<br>
<br>
Use o pacote car e a função recode.<br>
<br>
Em 13/12/2011 14:54, Edson Lira escreveu:<br>
> Boa tarde, caros r-istas, estou com alguns probleminhas na categorização<br>
> de variáveis.<br>
><br>
> Tenho este banco.<br>
> pad pas sexo pg<br>
> 155 134 83 feminino 15<br>
> 156 72 99 feminino 30<br>
> 157 82 151 feminino 10<br>
> 158 23 70 feminino 20<br>
> 159 75 128 masculino 13<br>
> 160 73 125 feminino 05<br>
> 161 125 80 masculino 08<br>
> 162 65 110 feminino 10<br>
> 163 132 92 masculino 13<br>
> 164 114 68 masculino 11<br>
> 165 83 153 feminino 26<br>
> 166 166 107 masculino 22<br>
> 167 114 76 feminino 03<br>
><br>
><br>
> Preciso classificar da seguinte forma:<br>
> se pad<85 e pas<130 é normal<br>
> se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
> e assim por diante.<br>
><br>
><br>
> a outra classificação é a seguinte:<br>
><br>
> Para homens o percentual de gordura<br>
><br>
> entre 15 e 20 é normal<br>
> entre 20 e 25 levemente obeso<br>
> entre 25 e 30 excessivamente obeso<br>
> entre 10 e 14 gordura ideal<br>
> 3 gordura essencial<br>
><br>
><br>
><br>
> Para homens o percentual de gordura<br>
><br>
> entre 25 e 30 é normal<br>
> entre 30 e 35 levemente obeso<br>
> entre 35 e 40 excessivamente obeso<br>
> entre 14 e 18 gordura ideal<br>
> 12 gordura essencial<br>
><br>
><br>
> Como fazer estaS categorizações?<br>
><br>
> []'s.<br>
><br>
> Edson Lira<br>
> Estatístico<br>
> Manaus-Amazonas<br>
><br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
<br>
<br>
--<br>
Cesar Rabak<br>
GNU/Linux User 52247.<br>
Get counted: <a href="http://counter.li.org/" target="_blank">http://counter.li.org/</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 10<br>
Date: Tue, 13 Dec 2011 10:22:32 -0800 (PST)<br>
From: "Elias T. Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Categorização de variáveis<br>
Message-ID:<br>
<<a href="mailto:1323800552.38193.YahooMailNeo@web114707.mail.gq1.yahoo.com">1323800552.38193.YahooMailNeo@web114707.mail.gq1.yahoo.com</a>><br>
Content-Type: text/plain; charset=iso-8859-1<br>
<br>
<br>
<br>
se vc enviasse um CMR seria mais fácil... veja este exemplo<br>
<br>
n <- 10<br>
x1 <- rnorm(n)<br>
<br>
### classificação considerando apenas uma variavel<br>
f1 <- factor(findInterval(x1, c(-Inf, -0.5, 0.5, Inf), labels=c("Ruim", "Médio", "Bom"))<br>
<br>
quando vc tem:<br>
<br>
se pad<85 e pas<130 é normal<br>
>se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
o que deve ser se pad<85 mas pas>130 por exemplo?<br>
<br>
Att.<br>
Elias T. Krainski<br>
<br>
<br>
<br>
>________________________________<br>
> De: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
>Para: R-br Lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
>Enviadas: Terça-feira, 13 de Dezembro de 2011 15:54<br>
>Assunto: [R-br] Categorização de variáveis<br>
><br>
><br>
>Boa tarde, caros r-istas, estou com alguns probleminhas na categorização de variáveis.<br>
><br>
><br>
>Tenho este banco.<br>
> pad pas sexo pg<br>
>155 134 83 feminino 15<br>
>156 72 99 feminino 30<br>
>157 82 151 feminino 10<br>
>158 23 70 feminino 20<br>
>159 75 128 masculino 13<br>
>160 73 125 feminino 05<br>
>161 125 80 masculino 08<br>
>162 65 110 feminino 10<br>
>163 132 92 masculino 13<br>
>164 114 68 masculino 11<br>
>165 83 153 feminino 26<br>
>166 166 107 masculino 22<br>
>167 114 76 feminino 03<br>
><br>
><br>
><br>
><br>
>Preciso classificar da seguinte forma:<br>
>se pad<85 e pas<130 é normal<br>
>se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
>e assim por diante.<br>
><br>
><br>
><br>
><br>
>a outra classificação é a seguinte:<br>
><br>
><br>
>Para homens o percentual de gordura<br>
><br>
><br>
>entre 15 e 20 é normal<br>
><br>
>entre 20 e 25 levemente obeso<br>
>entre 25 e 30 excessivamente obeso<br>
>entre 10 e 14 gordura ideal<br>
> 3 gordura essencial<br>
><br>
> <br>
><br>
><br>
><br>
><br>
><br>
>Para homens o percentual de gordura<br>
><br>
><br>
>entre 25 e 30 é normal<br>
><br>
>entre 30 e 35 levemente obeso<br>
>entre 35 e 40 excessivamente obeso<br>
>entre 14 e 18 gordura ideal<br>
> 12 gordura essencial<br>
><br>
> <br>
><br>
><br>
><br>
>Como fazer estaS categorizações?<br>
><br>
><br>
>[]'s.<br>
><br>
><br>
>Edson Lira<br>
>Estatístico<br>
>Manaus-Amazonas<br>
>_______________________________________________<br>
>R-br mailing list<br>
><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
><br>
> <br>
<br>
<br>
------------------------------<br>
<br>
Message: 11<br>
Date: Tue, 13 Dec 2011 15:37:02 -0400<br>
From: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>>, "Elias T.<br>
Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>><br>
Subject: Re: [R-br] Categorização de variáveis<br>
Message-ID: <<a href="mailto:63652F35-F09C-4C24-8D38-0F462C8DBBF2@yahoo.com.br">63652F35-F09C-4C24-8D38-0F462C8DBBF2@yahoo.com.br</a>><br>
Content-Type: text/plain; charset=utf-8<br>
<br>
O problema é que você tem duas variáveis, pas é pressão sistólica<br>
e pad é pressão diastólica. Se pas é menor que 130 e pad é menor que 180 o paciente é considerado normal e assim por diante.<br>
[].s<br>
Edson Lira<br>
Estatístico<br>
Ma-Am<br>
<br>
Em 13/12/2011, às 14:22, "Elias T. Krainski" <<a href="mailto:eliaskrainski@yahoo.com.br">eliaskrainski@yahoo.com.br</a>> escreveu:<br>
<br>
><br>
><br>
> se vc enviasse um CMR seria mais fácil... veja este exemplo<br>
><br>
> n <- 10<br>
> x1 <- rnorm(n)<br>
><br>
> ### classificação considerando apenas uma variavel<br>
> f1 <- factor(findInterval(x1, c(-Inf, -0.5, 0.5, Inf), labels=c("Ruim", "Médio", "Bom"))<br>
><br>
> quando vc tem:<br>
><br>
> se pad<85 e pas<130 é normal<br>
> >se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
> o que deve ser se pad<85 mas pas>130 por exemplo?<br>
><br>
> Att.<br>
> Elias T. Krainski<br>
><br>
><br>
><br>
>> ________________________________<br>
>> De: Edson Lira <<a href="mailto:edinhoestat@yahoo.com.br">edinhoestat@yahoo.com.br</a>><br>
>> Para: R-br Lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
>> Enviadas: Terça-feira, 13 de Dezembro de 2011 15:54<br>
>> Assunto: [R-br] Categorização de variáveis<br>
>><br>
>><br>
>> Boa tarde, caros r-istas, estou com alguns probleminhas na categorização de variáveis.<br>
>><br>
>><br>
>> Tenho este banco.<br>
>> pad pas sexo pg<br>
>> 155 134 83 feminino 15<br>
>> 156 72 99 feminino 30<br>
>> 157 82 151 feminino 10<br>
>> 158 23 70 feminino 20<br>
>> 159 75 128 masculino 13<br>
>> 160 73 125 feminino 05<br>
>> 161 125 80 masculino 08<br>
>> 162 65 110 feminino 10<br>
>> 163 132 92 masculino 13<br>
>> 164 114 68 masculino 11<br>
>> 165 83 153 feminino 26<br>
>> 166 166 107 masculino 22<br>
>> 167 114 76 feminino 03<br>
>><br>
>><br>
>><br>
>><br>
>> Preciso classificar da seguinte forma:<br>
>> se pad<85 e pas<130 é normal<br>
>> se 85<pad<89 e 130<pas<139 Normal limítrofe<br>
>> e assim por diante.<br>
>><br>
>><br>
>><br>
>><br>
>> a outra classificação é a seguinte:<br>
>><br>
>><br>
>> Para homens o percentual de gordura<br>
>><br>
>><br>
>> entre 15 e 20 é normal<br>
>><br>
>> entre 20 e 25 levemente obeso<br>
>> entre 25 e 30 excessivamente obeso<br>
>> entre 10 e 14 gordura ideal<br>
>> 3 gordura essencial<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> Para homens o percentual de gordura<br>
>><br>
>><br>
>> entre 25 e 30 é normal<br>
>><br>
>> entre 30 e 35 levemente obeso<br>
>> entre 35 e 40 excessivamente obeso<br>
>> entre 14 e 18 gordura ideal<br>
>> 12 gordura essencial<br>
>><br>
>><br>
>><br>
>><br>
>><br>
>> Como fazer estaS categorizações?<br>
>><br>
>><br>
>> []'s.<br>
>><br>
>><br>
>> Edson Lira<br>
>> Estatístico<br>
>> Manaus-Amazonas<br>
>> _______________________________________________<br>
>> R-br mailing list<br>
>> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
>><br>
>><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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.<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><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>
Fim da Digest R-br, volume 10, assunto 19<br>
*****************************************<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><font style="color:rgb(102,102,102);font-family:arial,helvetica,sans-serif" color="#999999" size="2"><div>-----</div><div>Fabiana Schneck</div></font><font style="color:rgb(102,102,102);font-family:arial,helvetica,sans-serif" color="#999999" size="2">Programa de Pós</font><font style="color:rgb(102,102,102);font-family:arial,helvetica,sans-serif" color="#999999" size="2">-Graduação em Ecologia</font><br style="color:rgb(102,102,102)">
<font style="color:rgb(102,102,102);font-family:arial,helvetica,sans-serif" color="#999999" size="2">Universidade Federal do Rio Grande do Sul</font> <span style="color:rgb(102,102,102)"><br>CP 15007 - CEP 91501-970</span> <br style="color:rgb(102,102,102)">
<span style="color:rgb(102,102,102)">Porto Alegre, RS, Brasil</span><br style="color:rgb(102,102,102)"><span style="color:rgb(102,102,102)"></span><span style="color:rgb(102,102,102)"></span><br><font style="color:rgb(102,102,102);font-family:arial,helvetica,sans-serif" color="#999999" size="2"></font><br style="color:rgb(102,102,102)">
<font color="#999999"><br></font><br>