Janeiro 2012 - R-br

Criando um processo de pontos sobre uma área regular
by ASANTOS 10 Jan '12

10 Jan '12

Boa tarde pessoal, Estou tentando criar um processo de pontos aleatório sobre uma área regular representada por árvores regularmente espaçadas, sendo que o processo de pontos só pode ocorrer onde não existam árvores e estou quebrando a cabeça com isso, bom primeiro criei a área com as árvores: require(plotrix) ##Cria o limite da área limx<-c(2,2,20.5,20.5) limy<-c(2,20.5,20.5,2) plot(limx,limy) lim=cbind(limx,limy) lim<-as.matrix(lim) polygon(lim) ##Aqui criei a posição das árvores e com área igual ao seu diâmetro sq=seq(2.5,20,2.5) for (k in sq) { for(j in sq){ draw.circle(k,j,radius=0.3,col="red") }} Bom agora eu queria jogar sobre essa área o processo de pontos: ##Criando o processo de pontos x <- runif(n=500,min=0, max=20) y <- runif(n=500,min=0, max=20) xy<-cbind(x,y) points(xy) ## Porém, os pontos também caem nas mesmas posições ocupadas pelas árvores (obviamente), onde gostaria de saber se alguém conheceria alguma solução para eu fazer com que não haja sobreposição entre os raios que representam os diâmetros das árvores e os pontos criados? Na verdade, preciso definir uma maneira de fazer com que as coordenadas contidas nos diâmetros das árvores não sejam validas no momento da geração do processo de pontos, atualmente estou tentando soluções com o pacote spatstat mais sem sucesso, pois não consigo fazer com que a área criada seja um objeto ppp valido. Obrigado, -- Alexandre DOS SANTOS Engenheiro Florestal, Msc. Laboratório de Entomologia Florestal Departamento de Entomologia Universidade Federal de Lavras Caixa Postal 3037 37200-000 - Lavras - Minas Gerais - Brasil Tel: +55 35 92230304

6 9

Fatorial com parcela subdivida
by Wenceslau Teixeira 10 Jan '12

10 Jan '12

Caros colegas da lista do R, um colega me pediu ajuda para analisar um experimento com um delineamento (complexo para meus conhecimentos) E um fatorial 2 x 3 x 3 com parcela subvididida (micronutrientes), feito em blocos com 3 repeticoes. ** Sera que alguem poderia nos indicar algum material do R sobre como fazer o modelo para este tipo de delineamento. Agradeco antecipadamente a ajuda Wenceslau Teixeira Embrapa Solos ** * * *FV GL .F* Blocos 2**** N 1**** P 2**** K 2**** NxP 2**** NxK 2* *** PxK 4**** NxPxK 4**** ** Erro a ?* *** ** ** Micro (?) 1 **** * * *Erro b (?) ? . * Total 53**** ** **

2 1

Como faço para atualizar minha versão do R
by andre neves 10 Jan '12

10 Jan '12

Prezados, Como faço para atualizar minha versão do R. Um abraço a todos. André =========================== Aviso de confidencialidade: =========================== - Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal n. 5.851, de 7 de dezembro de 1972, eh enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada eh ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. ===================== Confidentiality note: ===================== - "This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure".

3 4

Linux: Duas versões do R
by Jose Claudio Faria 09 Jan '12

09 Jan '12

Pessoal, É possível ter 2 versões do R (no caso 2.13.1 e 2.14.0-dev) instaladas e funcionais no notebook? A 2.14.0 eu compilo dos fontes, a outra baixo já compilada. Sempre usei as dev no Linux, mais ela está apresentando conflito com alguns pacotes que necessito para rodar uma análise. Ando procurando essa informação mas ainda não encontrei: Qualquer dica é bem vinda. Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 3

post: invocando R através do Stata
by Daniel Marcelino 09 Jan '12

09 Jan '12

A quem possa interessar, acabo de publicar um post sobre como rodar o R através do Stata em Mac OS. Assim, talvez o post sirva para quem use os dois pacotes statisticos em Mac OS. http://danielmarcelino.com/

1 0

post: invocando o R através do Stata
by Daniel Marcelino 09 Jan '12

09 Jan '12

A quem possa interessar, acabo de publicar um post sobre como rodar o R através do Stata em Mac OS. Assim, talvez o post sirva para quem use os dois pacotes statisticos em Mac OS. http://danielmarcelino.com/

1 0

Erro em shapiro.test
by Flávia Machado 06 Jan '12

06 Jan '12

Pessoal, fui fazer um novo teste de shapiro e encontrei outro erro. Fiz: > dados<-read.table(file.choose(),header=T) > attach(dados) > dados apareceu: Erro em shapiro.test(FrMad, FrImat) : argumento(s) não utilizado(s) (FrImat) Coloquei em anexo a planilha que estou usando. Dessa vez não esqueci de "attachar" o arquivo. Alguém saberia me dizer o motivo do erro? Att. -- Flávia Machado

3 3

Software publico brasileiro
by Pedro Rafael 06 Jan '12

06 Jan '12

http://linuxmagazine.com.br/lm/noticia/software_publico_brasileiro_na_linux… -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

ajuda boxplo
by Ana Paula Sousa 06 Jan '12

06 Jan '12

Colegas, Fiz um gráfico boxplot de séries temporais ajustadas por mês. Na escala horizontal, ele entende os valores com 1,2...,12. Preciso que 1=jan;2=fev e assim por diante. Como faço? Atenciosamente,

3 2

Ajuda com pacote "seas"
by Inajara Welter Macedo 06 Jan '12

06 Jan '12

Ola, Estou trabalhando com dados de precipitacao, de 10 anos, e quero obter resultados de periodos de seca. Preciso contar quantos dias consecutivos sem chuva houve, ao todo, e por ano e depois por mes. Meus dados eram uma matriz de 4 colunas, chamado "mod1", cada coluna corresponde ao ano, mes, dia, e precipitacao, mas transformei para data.frame, para poder utilizar o pacote "seas": as.data.frame(mod1) Tentei usar a funcao "interarrival": CCS<-interarrival (mod1, var = "Precip", p.cut = 0, inv = TRUE) # CCS= casos consecutivos secos Mas, ocorre um erro: Erro em seas.df.check(x, orig, var) : a ‘date’ column must exist in ‘mod1’ Nao sei como colocar essa coluna de data que pede, porque existem as colunas de ano , mes, e dia. Como consigo converter essas colunas em uma unica coluna de data?? Obrigada, Inajara - UFPEL/RS

5 4