- R-br - DINF UFPR

Pintando os ladinhos de uma distribuição normal com o comando polygon {graphics}
by Augusto Ribas 01 Abr '11

01 Abr '11

Olá pessoal. Bem eu estava lendo um livro aqui, o The R Book do M. J Crawley(2007) E na pagina 212 tem uma figura que to tentando reproduzir mas num to conseguindo. A ideia é mais ou menos o seguinte. #Primeiro temos la um plot de uma curva de distribuição normal com 3 desvios: x<-seq(-3,3,0.01) y<-exp(-abs(x)^2) plot(x,y,type="l") #Dai se no exemplo dele ele fala que se trata de dados de alturas de pessoas, onde a media é 170cm e desvio padrão de 8cm #dai eu mudei o eixo X da seguinte forma altura<-170+(8*x) plot(altura,y,type="l") #dai o grafico fica com a a distribuição baseado no dados em questão altura e nao em unidades de desvio. #então chega o meu problema, ele nao deixo codigo de como faze o pintadinho de ver a chande de alguma coisa. #por exemplo, alturas menores de 162, qual a chance de ver alguem menor de 162. #eu vi naquelas listas de graficos que tem codigo que da pra marcar isso com o comando polygon #mas eu nao to conseguindo usalo direito pra isso, e no help dele os exemplos nao me ajudaram muito #algo assim plot(altura,y,type="l") polygon(c(146,162,162),pnorm(c(-3,-1,-3)),col="black") #mas eu nao quero fazer um triangulo, tinha que ser aquela curvinha seguindo os valores de y, mas ai como ter varios #valores de y ao longo dos de x pra fica bonitinho que eu nao to conseguindo, acho que olhando a figura da pra entender #minha duvida Obrigado a todos. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 1

Criar arquivo
by Victor Eduardo 01 Abr '11

01 Abr '11

Galera, além do Write, existe algum outro arquivo que faça eu criar um arquivo em txt? Quero armazenar as linhas de uns arquivos que selecionei, por exemplo, tem algumas das delas aqui: V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 9 8200 2010 5 22 16 0 -49.7 -29.4 4 -9999 1019 1019 106 1.6 294.6 15 8203 2010 5 22 16 0 -52.7 -31.4 464 -9999 967 1021 60 4.3 291.9 24 8238 2010 5 22 16 0 -54.5 -25.6 231 -9999 988 1014 68 2.2 296.8 26 8240 2010 5 22 16 0 -53.6 -26.3 810 -9999 924 1016 63 5.0 293.2 64 8209 2010 5 22 16 0 -49.3 -26.9 86 -9999 1009 1019 176 1.3 293.8 68 8204 2010 5 22 16 0 -53.5 -30.5 450 -9999 970 1022 67 7.5 293.4 Abraços!

3 5

string
by Daniel Marcelino 01 Abr '11

01 Abr '11

Caros, como eu poderia escrever um loop para retirar todos os caracateres que estão antes e depois da "/" na string abaixo e colocando-os em novas variáveis? Ex. data <-c("155 - VITAL DO REGO FILHO - PB - Senador", "1111 - - PP - - Deputado Federal / 25888 - ATAIDES MENDES PEDROSA - PB - Deputado Estadual", "1111 - - PP - - Deputado Federal / 22333 - EDNALDO PEREIRA DE SANTANA - PB - Deputado Estadual", "151 - JOSE WILSON SANTIAGO - PB - Senador", "45123 - ANTONIO HERVAZIO BEZERRA CAVALCANTI - PB - Deputado Estadual", "1212 - DAMIÃO FELICIANO DA SILVA - PB -Deputado Federal")

5 7

Res: Res: Res: [R_STAT] Dados não balanceados!
by Walmes Zeviani 01 Abr '11

01 Abr '11

Thiago, Se no seu modelo você tiver que estimar efeitos da interação clone:local, um dos efeitos será não estimável por não possuir informação amostral para ele, no coef() do modelo vai haver um NA para esse termo. Por outro lado, se essa interação clone:local for desprezada (assumida ser zero, não fizer sentido, etc), o desbalanceamento não vai impedir você de estimar os efeitos principais. No primeiro caso, caso exista interação, você terá dificuldade para testar efeitos principais, embora não exista o mínimo de sentido testar os efeitos principais na presença de interação, no entanto alguns softwares usam algumas restrições e coisa e tal e fazem testes para os efeitos principais na presença de interação em delineamentos com perda completa de casela, um esforço um tanto desnecessário. Se a interação é não significativa você pode ajustar um modelo com apenas os efeitos principais e estima-los e testar hipóteses a respeito tranquilamente. Outra saída é considerar o efeito de local e suas interações como aleatório. Daí a coisa saí fácil porque a estimação vai importar informação das demais observações. À 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: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== *De:* Thiago De paula protásio <depaulaprotasio(a)yahoo.com.br> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Quinta-feira, 31 de Março de 2011 19:20:24 *Assunto:* Re: Res: Res: [R_STAT] Dados não balanceados! Pessoal; Necessito da grande ajuda de vcs.... Tem um delineamento em esquema fatorial, mas os dados estão desbalaceados. Nesse caso, um dos fatores não está completo.... Eu tenho 4 clones, 3 locais e 2 espaçamentos, mas para um dos clones eu não tenho todos os locais. Como posso proceder com essa estatística no R? O Walmes já me deu alguns toques sobre algo parecido a algum tempo.... Conto com vcs! Muito obrigado! --- Em *sáb, 12/2/11, walmes zeviani <walmeszeviani(a)yahoo.com.br>* escreveu: De: walmes zeviani <walmeszeviani(a)yahoo.com.br> Assunto: Res: Res: [R_STAT] Dados não balanceados! Para: R_STAT(a)yahoogrupos.com.br Data: Sábado, 12 de Fevereiro de 2011, 13:21 Thiago, Como o Ivan disse, a coisa do teste dos termos do modelo ou anova, pode ser empregada a car::Anova(..., type="III"). Mas os cuidados com desbalanceamentos atravessam a fronteira da anova. Os testes de médias (ou demais inferências) não são mais baseados nas médias amostrais, uma vez que essas são viesadas devido a perda de ortogonalidade causada pelo desbalanceamento. O que deve ser obtido são as médias ajustadas que requerem manipulação matricial do vetor de solução do sistema linear, ou coef(modelo.lm). Antes de tudo isso, deve se saber se o seu desbalanceamento não comprometeu a estimação dos efeitos, pois existem duas "intensidades" de desbalanceamento: * apenas perda de repetições: aqui toda combinação entre níveis dos fatores tem pelo menos uma repetição com observação, e.g. fator A (A1,A2,A3) e B (B1,B2) todas as 6 combinações de A com B possuem pelo menos uma observação. Aqui todos os efeitos podem ser estimados embora com erros padrões diferentes em função do desbalanceamento. * perda de caselas: aqui pelo menos uma combinação dos níveis dos tratamentos não foi observada, e.g. fator A (A1,A2,A3) e B (B1,B2) alguma das 6 combinações foi completamente não observada (e.g. A1B1). Nesse caso não é possível estimar alguns efeitos justamente por não haver as observações experimentais correspondentes para estimá-los. No caso de dois fatores sem interação entre eles tem um exemplo na página de Rídiculas do LEG que fiz ontem para outro fim. Quem sabe possamos fazer uma Rídicula com o seu problema. http://www.leg.ufpr.br/doku.php/ridiculas À disposição. Walmes. ============================================================ Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br / @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes ============================================================ ------------------------------ *De:* Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Sábado, 12 de Fevereiro de 2011 10:38:53 *Assunto:* Res: [R_STAT] Dados não balanceados! Bom dia Thiago!! Nestes casos é aconselhável vc utilizar a soma de quadrados do tipo III. Para obter resultados semelhantes a alguns softwares consagrados vc pode utilizar a reparametrização do tipo soma zero. Para isso: options(contrasts = c("contr.sum", "contr.helmert")) library(car) modelo <- lm( y ~ A*B) Anova(modelo, type="III") Abraço. (S,f,P) Allaman *M.Sc Ivan Bezerra Allaman* Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA **email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9925-6428 ------------------------------ *De:* Thiago De paula protásio <depaulaprotasio(a)yahoo.com.br> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Sábado, 12 de Fevereiro de 2011 10:20:28 *Assunto:* [R_STAT] Dados não balanceados! Olá Pessoal; Estou precisando de algumas orientações para rodar um um fatorial com dados desbalaceados. Espero que alguém posso me ajudar com alguma função ou material. Muito obrigado! Bom final de semana a todos!!

1 0

Implementação em C
by Fernando Henrique Toledo 01 Abr '11

01 Abr '11

Lista, boa noite. Gostaria de pedir a quem possa algumas dicas e tutoriais que me posibilitem um maior entendimento do processo de implementação de funções em C e uso dessas mesmas no ambiente R! Na verdade não programo em C, mas estou enfrentando certos problemas com algumas funções que crio no R, pois o tempo de processamento delasestá um pouco impeditivo. Li alguns tutoriais do R e do GCC, mas ainda não tenho conhecimento suficiente para enveredar por esse lado. Assim gostaria de saber da experiencia de vocês quanto a isso. fico no aguardo, att, Fernando H

7 31

Gráfico de surface
by Ana Lúcia Souza Silva Mateus 01 Abr '11

01 Abr '11

Como fazer um gráfico de surface no R? Grata Ana

2 2

Testar ângulos de componentes principais
by Ivan Bezerra Allaman 31 Mar '11

31 Mar '11

Boa tarde senhores! Mais uma vez venho recorrer ao auxílio dos senhores. Estou querendo testar se o ângulo entre componentes principais de 3 populações são iguais a zero. Suponhamos 3 populações: POP1, POP2 e POP3. Para cada população eu calculei um componente principal (PC). O interesse é se avaliar se um determinado estudo pode ser feito utilizando um único componente principal, considerando os dados agrupados, ou se é necessário fazer o estudo em separado para cada população. Segundo o material que estou seguindo, o autor calculou o ângulo entre cada componente principal de cada população e testou se estes são iguais a zero por meio do procedimento bootstrap. Chamando de @ o ângulo entre dois componentes principais, a hipótese de nulidade é: H0: @(POP1,POP2) = @(POP1,POP3) = @(POP2,POP3) = 0 Segundo o autor, o uso do procedimento bootstrap é direto. Nas suas palavras: "The bootstrap test procedure is straightforward: use the data to generate a modified data set that conforms to the null hypothesis, then repeatedly draw bootstrap samples and compute the test statistic, and finally compare the empirical distribution of the test statistic to the value calculated for the original data (Efron and Tibshirani, 1993)." É a primeira vez que utilizo o procedimento bootstrap, e estou tendo dificuldades de digerir o procedimento descrito pelo autor acima. Já li o livro do Efron (o tópico teste de hipóteses), mas não consigo entender como o procedimento consegue testar a hípotese de nulidade. Se for possível me darem uma luz eu agradeço. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Blocos incompletos no R
by Julia Maria 31 Mar '11

31 Mar '11

O experimento é em parcela subdividida com tratamentos principais dispostos em blocos incompletos balanceados. (Tese Iemma) Ele se dispõe da seguinte forma: V tratamentos principais dispostos em A blocos incompletos balanceados, com R repetições para cada tratamento principal, K tratamentos principais por bloco e U tratamentos secundários. Assim no experimento eu teria: V=4 tratamentos principais no experimento todo; K=2 tratamentos principais por bloco; R=3 repetições para cada tratamento principal; A=6 blocos; U=5 tratamentos secundários no ensaio. Assim, na primeira repetição do experimento eu teria somente 2 tratamentos primarios ao invés dos 4 (o que esta caracterizando o BIB) e dentro de cada tratamento primario, tenho os 5 secundários. Na segunda repetição, eu teria outro par de tratamento primario e dentro deles os 5 secndarios. Estou enviando uma figura com um tipo de croqui, do jeito que se encontra na tese.

3 3

Margem de saída de gráfico
by Ana Lúcia Souza Silva Mateus 31 Mar '11

31 Mar '11

Bom dia .... Tenho necessidade de colocar 22 gráficos em uma página de saída do R, do jeito que estou fazendo so cabe 9 (para uma boa visualização gráfica). Qdo peço um par(mfrow=c(5,5)) ele da a seguinte mensagem: margens da figura mto grande. Existe uma forma de aumentá-la. Grata Ana

2 2

{off topic) ubuntu usando wubi?
by Luis Iván Ortiz Valencia 31 Mar '11

31 Mar '11

Boa noite... Estou tentando instalar o Ubuntu no meu micro com windows. Voces recomendam usar o WUBI para instalar ou e melhor usar o caminho usual de particionar o hd ( este ultimo sendo muito confuso para mim ainda)? ivan

5 5