Fevereiro 2015 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

2 2

Estatísticas do R
by Andre Oliveira 26 Abr '16

26 Abr '16

Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

8 8

Fwd: Enc: model.matrix()
by Ze Henrique 23 Abr '15

23 Abr '15

> Assunto: model.matrix() > > Bom tarde pessoal, > > gostaria de saber se é possível transformar um nível de um fator em numérico no próprio vetor, de forma que ele não seja codificado na função model.matrix(). > > Segue abaixo um CMR: > > k = rep(0:3,each=3) > k = factor(k) > k = model.matrix(~-1+k) > > # Essa é a matriz gerada. > > k0 k1 k2 k3 > 1 1 0 0 0 > 2 1 0 0 0 > 3 1 0 0 0 > 4 0 1 0 0 > 5 0 1 0 0 > 6 0 1 0 0 > 7 0 0 1 0 > 8 0 0 1 0 > 9 0 0 1 0 > 10 0 0 0 1 > 11 0 0 0 1 > 12 0 0 0 1 > > # Entretanto gostaria de obter a seguinte matriz sem a necessidade de eliminar as colunas posteriormente. > k0 k3 > 1 1 0 > 2 1 0 > 3 1 0 > 4 0 0 > 5 0 0 > 6 0 0 > 7 0 0 > 8 0 0 > 9 0 0 > 10 0 1 > 11 0 1 > 12 0 1 > > Att. > > Tiago. > > ################################################################# > > Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES > > Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas > > #################################################################

7 22

Probit/CL50
by Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA 10 Abr '15

10 Abr '15

Olá pessoal, Preciso de orientação para fazer a análise de Probit, estimar a CL50 e fazer o gráfico. Meu experimento avaliou o efeito inseticida de doses de uma substância. Se puderem...será muito válido! Flávia Batista ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e 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.

5 11

Empilhar variáveis em colunas para linhas
by ASANTOS 22 Mar '15

22 Mar '15

Caros Listeiros, Tenho um conjunto de dados artificiais abaixo, constituído de dados de contagem de três espécies (Esp1, Esp2 e Esp3) e duas variáveis trat e tempo. As três espécies e as variáveis estão dispostas em colunas, sendo: #Dados artificiais ----------------------------------- da <- expand.grid(trat=gl(4,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) da$Esp1<-rpois(da[,1],5) da$Esp2<-rpois(da[,1],10) da$Esp3<-rpois(da[,1],15) head(da) trat tempo Esp1 Esp2 Esp3 1 1 1 4 13 8 2 2 1 7 7 17 3 3 1 3 12 16 4 4 1 8 8 13 5 1 2 4 13 22 6 2 2 4 8 10 Porém, eu gostaria de reorganizar o banco de dados e substituir as três espécies (Esp1, Esp2 e Esp3), por duas colunas, sendo uma coluna com o nome das espécies e uma nova coluna para os dados de contagem, representada por N, de forma que mantenha a correspondência entre as variáveis trat e tempo, ficando da seguinte maneira: head(new.da) trat tempo Esp N 1 1 1 Esp1 4 2 1 1 Esp2 13 3 1 1 Esp3 8 4 2 1 Esp1 7 5 2 1 Esp2 7 6 2 1 Esp3 17 Alguém poderia me indicar alguma função ou procedimento para fazer isso? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

3 5

krige.conv
by Felinto COSTA 07 Mar '15

07 Mar '15

Prezados. Tenho dois modelos teóricos (exponencial e esférico ) que foram ajustados aos dados por máxima verossimilhança. As coordenadas numéricas (UTM) dos dados foram reduzidas pela subtração da coordenada média segundo cada eixo. Em ambos considero a existência de uma tendência polinomial de 2ª ordem na média. Eles foram ajustados sem problemas, assim como a validação cruzada de ambos pode ser feita. O problema está ao querer estimar algum ponto(s) usando a função krige.conv. Apenas o modelo exponencial é aceito. Ao usar o esférico a função retorna a seguinte mensagem: krige.conv: model with mean given by a 2nd order polynomial on the coordinates krige.conv: anisotropy correction performed krige.conv: performing the Box-Cox data transformation Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 1.76896e-16 Alguém saberia como contornar esse problema? Antecipadamente agradeço a colaboração Felinto COSTA

3 3

Extrair coordenadas de um raster dada uma condição
by ASANTOS 05 Mar '15

05 Mar '15

Caros Listeiros, Fiz a seleção da coordenada dos pixels (p.cen) em um raster que continham pontos de interesse (pontos) no CRM abaixo, gostaria de saber como impor uma condição para selecionar as coordenadas dos pixels que não contém nenhum ponto. Tentei criar um novo objeto com todas as coordenadas no raster e remover as coordenadas p.cen sem sucesso, alguém poderia me ajudar? Obrigado, ### <code r> require(raster); require(sp) ## Criando um raster r <- raster(nc=10, nr=10) r <- setValues(r, round(runif(ncell(r))* 255)) ##Pontos selecionados x <- c(44.25,33.06,63.22,71.35) y <- c(22.00,81.90,15.18,71.03) pontos <- SpatialPoints(cbind(x,y)) cells <- cellFromXY(r, pontos); cells coordinates(r)[cells,] p.cen <- xyFromCell(r, cells); p.cen plot(r) points(pontos, col="red") points(p.cen) values <- extract(r,p.cen) ##Extrai as coordendas dados<-cbind(p.cen,cells,values) dados ### </code> -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

3 4

como salvar gráfico após o iso do comando identify
by Antonio Silva 03 Mar '15

03 Mar '15

Olá a todos Quero salvar um gráfico após o uso do comando identify mas não consegui. Utilizo o comando identify para adicionar etiquetas aos pontos de forma manual, para que os números não fiquem embolados ou sobre as linhas. Exemplo x<-c(1,2,3,4,5,6) y<-c(20,30,15,7,25,40) plot(x,y,type="b",ylim=c(0,45)) identify(x,y,labels=y) No entanto, quando escrevo png(filename="teste.png",width=20,height=20,units="cm",res=300) x<-c(1,2,3,4,5,6) y<-c(20,30,15,7,25,40) plot(x,y,type="b",ylim=c(0,45)) identify(x,y,labels=y) dev.off() o gráfico é gravado sem aparecer na tela e consequentemente sem a identificação dos pontos. Utilizo o R integrado ao gedit no UBUNTU Trusty Tahr 64 bits. Gostaria de fazer o procedimento pela linha de comando, sem a utilização de GUIs. Há alguma forma de gravar o gráfico com os pontos identificados? Obrigado! Antônio Olinto

3 3

Mover os vértices de um polígono
by ASANTOS 03 Mar '15

03 Mar '15

Caros Listeiros, Criei um polígono (c1) abaixo com os dados dos vértices da área em UTM e gostaria de saber se alguém conhece alguma função para mover todos os vértices 30m para dentro e reduzir a área do polígono? #Contorno da área -------------------------------------------------------------- x.coords1<-c(371299.9, 371266.4, 371205.6, 371111.8, 371047.6, 371018.2, 371014.0, 371009.3, 370983.1, 370919.7, 370853.6, 370785.6, 370748.8, 370711.8, 370687.8, 370696.4, 370785.9, 370885.5, 371035.8, 371148.1, 371205.2, 371231.7, 371236.5, 371240.3, 371285.8, 371326.5, 371397.2, 371417.1, 371432.9, 371445.0, 371455.7, 371466.4, 371476.6, 371502.6, 371536.0, 371550.0, 371546.8, 371528.3, 371470.0, 371393.3, 371299.9, 371299.9) y.coords1<-c(8246589, 8246560, 8246508, 8246428, 8246373, 8246349, 8246348, 8246352, 8246385, 8246465, 8246551, 8246638, 8246685, 8246732, 8246764, 8246771, 8246846, 8246932, 8247062, 8247160, 8247209, 8247230, 8247224, 8247221, 8247160, 8247107, 8247016, 8246991, 8246967, 8246939, 8246914, 8246892, 8246875, 8246846, 8246821, 8246809, 8246802, 8246785, 8246735, 8246669, 8246589, 8246589) # Cria um poligono t05<-unique(data.frame(x=rev(x.coords1),y=rev(y.coords1))) c1 = cbind (t05$x, t05$y) r1 = rbind (c1, c1[1, ]) # junta plot(c1) polygon(c1) # -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 1