- R-br - DINF UFPR

sobrescrito
by Peterson Gardin 20 Jun '13

20 Jun '13

Olá pessoal 1) Como colocar subscrito para fórmula de nitrato NO3, colocar o 3 em subescrito. 2) Como colocar sobrescrito no g-1 e h-1. Comandos gráficos abaixo. bar.group(comparison$groups,horiz=FALSE,col=rainbow(5),border="blue", xlim=c(0,6), ylim=c(0,2.8),las=1) title(cex.lab=1.2, xlab="Porcentagem de N", ylab="Micromol NO3 g-1 h-1") João Peterson Pereira Gardin Doutor em Fisiologia Vegetal, Eng. Agr., Pesquisador Fones: (49) 3533 5626 e 9933-4290 MSN: petersongardin(a)yahoo.com.br - SKYPE: petersongardin Endereço: Rua João Zardo, 1660, Bairro Campo Experimental, Videira - SC - CEP 89560-000

4 7

Comparar dois data frames
by Fátima Lima Paula 20 Jun '13

20 Jun '13

Olá pessoal. Existe alguma forma de comparar dois data frames para saber se há alguma diferença entre eles? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 6

Fábio, please add me to your LinkedIn network
by MÁRCIO SOUZA 20 Jun '13

20 Jun '13

LinkedIn ------------ MÁRCIO SOUZA requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from MÁRCIO SOUZA http://www.linkedin.com/e/-33g73l-hi6et0wh-1d/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of MÁRCIO SOUZA http://www.linkedin.com/e/-33g73l-hi6et0wh-1d/rso/265728365/Zfs1/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hi6et0wh-1d/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 20 Jun '13

20 Jun '13

Aos membros do R-br por favor, voltando ao meu propósito inicial eu preciso apenas de uma orientação em relação a uma função que retorne o índice (i,j) e o valor do pixel de todo raster no R, nada mais do que isto. O R possui bastante recursos para executar esta tarefa, mas infelizmente não encontrei solução. Depois disto exportar/converter para shapefile será mais simples. Grato, Alex Santos

3 3

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 20 Jun '13

20 Jun '13

Elias, informo que minhas amostas não são pixels. Minha metodologia está consolidada, preciso apenas de uma diretriz para encontrar está função no R (pois acredito que o mesmo seja capaz para realizar o que desejo. Abaixo, está em linha gerais minha metodologia, mas por favor, não se atente a mesma pois meu foco é encontra/contrair a função que extraia o índice (i,j) de cada pixel. 1 - Eu tenho uma articulação do mapeamento sistemático que recobre minha área de interesse (no formato shapefile - poligono). Sobre estas áreas que houve a amostragem sistemática. (obs: Numa solução paralela utilizando o FME (software comercial) já agreguei as informações de (i,j) em cada quadrícula, podendo utilizar esta informação como X, Y e representa minhas quadrículas de forma contínua. Mas quero substituir o FME (comercial) pelo R (livre), pois acredito no potencial do R e acho ele capaz de executar esta tarefa. 2 - converti a articulação para raster, já com a proposta de encontrar as informações do índice (i,j) de cada pixel. Sendo que cada pixel possui o mesmo tamanho de minhas quadrículas. Mas espero que compreenda, pois o que desejo realmente é uma função que leia um arquivo raster - extraia o índice (i,j) de cada pixel (compondo uma matriz com os índices mais o valor do atributo) e exporte esta matriz para um shapefile. obs.: Adianto que o mesmo não será georreferenciado e nem há necessidade neste momento, pois haverá o georreferenciamento após a krigagem através dos extremos da articulação original. Att, Alex Santos

2 1

sortear 100x e correlacionar
by Hélio Gallo Rocha 20 Jun '13

20 Jun '13

Bom dia a todos, Peço ajuda na seguinte questão destes dados que comparam dois modelos de previsão de safra em café, onde : colheita = dados da parcela colhida modelo 1 de previsão modelo 2 de previsão CRM coleta=c(rep(1,25),rep(2,25)) colheita=c(10.2,9.22,8.6,15.57,11.93,12.28,12.43,16.47,18.72,12.77,16.03,13.4,16.95, 15.4,18.64,14.38,15.02,16.98,14.38,16.35,18.28,16.06,17.92,12.85,13.83,4.48,3.35,7.24, 7.73,8.45,4.83,3.8,2.2,5.99,4.24,4.67,8.31,4.25,7.32,5.45,4.24,5.72,5.09,4.84,5.62,5.5, 4.54,4.72,7.69,3.45) modelo1=c(17.08,12.24,13.44,18.52,11.82,11.37,16.17,16.92,16.66,13.32,18.95,14.78, 15.54,16.12,25.62,18.2,15.86,20.5,20.47,20.45,17.15,22.15,20.29,13.82,12.86,3.13, 1.84,6.9,6.35,8.9,5.64,4.79,2,4.64,5.06,4.85,5.61,4.82,6.59,4.15,4.46,5.1,4.78, 4.05,5.02,5.65,3.04,3.61,9.48,4.13) modelo2=c(12.43,10.99,10.32,12.46,9.85,10.15,12.68,12.46,12.16,9.92,12.58,10.63, 11.18,10.6,15.71,13.22,10.82,12.85,13.48,14.13,11.41,14.09,12.23,10.28,8.75, 5.2,3.56,9.89,9.03,12,8.16,7.15,3.2,6.58,6.84,7.37,8.5,7.87,8.38,5.58,6.62, 6.36,6.84,6.48,7.58,7.35,4.99,5.3,11.99,6.21) dados=cbind(coleta,colheita,modelo1,modelo2);dados x=cbind(colheita,modelo1,modelo2) cor(x) #sorteio de 20 medidas, r=do.call(rbind, by(dados, dados[,1], function(x, ...) x[sample(1:nrow(x), ...), ], size=10)) y=cbind(r$colheita,r$modelo1,r$modelo2);y cor(y) Gostaria de sortear 100 vezes com 20 amostras, sendo o resultado a média da correlação e a média dos dados Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 4

Funções beta-generalizada, beta-generalizada mononuclear, Richards ou Beta-Richards no R
by Ederson Civardi 20 Jun '13

20 Jun '13

Bom dia a todos. Gostaria de obter ajuda de algum de vocês que já trabalharam ou utilizaram alguma dessas funções (beta-generalizada, beta-generalizada mononuclear, Richards ou Beta-Richards) ou outras similares/adequadas para ajustar curvas de regressão não linear para o R. Estou avaliando tamanho de lesão de uma doença em função da temperatura (10 - 15 -20 -25 e 30 graus C) e molhamento foliar de 4 h até 32h, com intervalo de 4h. Preciso gerar superfície de resposta para estes fatores e não sei como fazer isso no R. Além já trabalhou com isso? E como posso fazer no R? Desde já agradeço a qualquer resposta obtida. Ederson --

1 0

teste wilcoxon
by alanarocha＠sapo.pt 20 Jun '13

20 Jun '13

----- Mensagem encaminhada de alanarocha(a)sapo.pt ----- Data: Tue, 18 Jun 2013 17:29:28 +0100 De: alanarocha(a)sapo.pt Assunto: teste wilcoxon Para: r-br(a)listas.c3sl.ufpr.br Olá a todos, executei a base de dados original: DoseT-9 DoseT-6 DoseT-3 DoseT+3 DoseT+6 DoseT+9 1 54020104 930,0 900,0 900,0 390,0 1.170,0 2 54020178 1.860,0 930,0 900,0 420,0 840,0 3 54020303 1.860,0 2.790,0 2.700,0 1.620,0 1.440,0 4 54020492 930,0 930,0 900,0 420,0 840,0 5 54020637 1.080,0 480,0 870,0 2.970,0 2.700,0 6 54020698 2.790,0 1.590,0 930,0 390,0 720,0 NA 1 2.700,0 2 630,0 3 1.530,0 4 720,0 5 2.520,0 6 2.340,0 depois tentei executar só duas colunas da base de dados onde quero aplicar o teste de wilcoxon para elas o resultado não entendo: > setwd("f:/Ana_trabalho/Cinacalset") > dados<-read.table("Cinacalset.csv",header=TRUE,sep=";",dec=",") > colnames(dados)<-c("DoseT-9", "DoseT-6", "DoseT-3", "DoseT+3", > "DoseT+6", "DoseT+9") > dados[,-c(1,2,5,6,7)] DoseT-3 DoseT+3 1 900,0 900,0 2 930,0 900,0 3 2.790,0 2.700,0 4 930,0 900,0 5 480,0 870,0 6 1.590,0 930,0 7 930,0 900,0 8 750,0 780,0 9 3.720,0 3.600,0 10 420,0 60,0 11 2.790,0 2.790,0 12 930,0 900,0 13 3.720,0 3.720,0 14 930,0 900,0 15 930,0 1.800,0 16 900,0 930,0 17 1.860,0 1.860,0 18 1.800,0 1.800,0 19 390,0 270,0 20 750,0 780,0 21 930,0 900,0 22 930,0 930,0 23 540,0 450,0 24 390,0 390,0 25 1.860,0 1.860,0 26 2.340,0 2.340,0 27 210,0 390,0 28 390,0 930,0 29 1.440,0 1.350,0 30 930,0 900,0 31 1.410,0 1.560,0 32 840,0 930,0 33 900,0 930,0 34 1.860,0 1.710,0 35 900,0 600,0 36 900,0 1.860,0 37 900,0 930,0 38 900,0 930,0 39 900,0 930,0 40 900,0 930,0 41 390,0 390,0 42 900,0 930,0 43 930,0 900,0 44 2.700,0 2.790,0 45 900,0 930,0 46 900,0 930,0 47 900,0 930,0 48 900,0 930,0 49 900,0 1.860,0 50 1.800,0 930,0 51 1.800,0 930,0 52 3.600,0 3.720,0 53 1.800,0 1.860,0 54 1.860,0 1.860,0 55 1.860,0 1.860,0 56 930,0 900,0 57 360,0 930,0 58 390,0 390,0 59 1.800,0 1.860,0 60 900,0 1.860,0 61 900,0 930,0 62 930,0 930,0 63 930,0 930,0 64 930,0 900,0 65 930,0 900,0 66 900,0 930,0 67 930,0 900,0 68 2.790,0 2.790,0 69 1.860,0 1.800,0 70 390,0 900,0 71 930,0 900,0 72 360,0 390,0 73 900,0 930,0 74 930,0 930,0 75 2.520,0 2.040,0 76 900,0 840,0 77 1.680,0 2.820,0 78 1.680,0 1.860,0 79 840,0 930,0 80 930,0 930,0 81 1.680,0 1.860,0 82 1.860,0 1.800,0 83 900,0 840,0 84 930,0 1.860,0 85 930,0 900,0 86 930,0 390,0 87 360,0 360,0 88 930,0 930,0 89 930,0 900,0 90 930,0 900,0 91 840,0 1.860,0 92 360,0 360,0 93 1.800,0 1.680,0 wilcox.test(dados[,-c(1,2,5,6,7)]$DoseT-3, + dados[,-c(1,2,5,6,7)]$DoseT+3, paired=TRUE) Error in wilcox.test.default(dados[, -c(1, 2, 5, 6, 7)]$DoseT - 3, dados[, : not enough (finite) 'x' observations Cumprimentos Ana Quoting r-br-request(a)listas.c3sl.ufpr.br: > Enviar submissões para a lista de discussão R-br para > r-br(a)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(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)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. trabalho em R test wilcox ou regressão (alanarocha(a)sapo.pt) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 18 Jun 2013 11:06:30 +0100 > From: alanarocha(a)sapo.pt > To: r-br <r-br(a)listas.c3sl.ufpr.br>, fjsequeira(a)fc.ul.pt, > sandro.matos(a)vodafone.com, rcoster(a)gmail.com, > cezar-fonseca(a)hotmail.com > Subject: [R-br] trabalho em R test wilcox ou regressão > Message-ID: > <20130618110630.Horde.PzJaJaeWCxtRwDEmmikGXKA(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > > Bom dia tenho a seguinte questão a resolver, alguém me sabe dizer se > devo utilizar o teste de wilcoxon emparelhado ou um modelo de > regressão linear > 1. Avaliar as médias para cada caso, considerando 2 períodos > a. Período antes é igual aos Doset-9 ao Doset-3, que passará a ser > designado como H > b. Período depois é igual aos Doset+3 ao Doset+9, que passará a ser > designado como C > DoseT-9 DoseT-6 DoseT-3 DoseT+3 DoseT+6 DoseT+9 > 930,0 900,0 900,0 390,0 1.170,0 2.700,0 > 1.860,0 930,0 900,0 420,0 840,0 630,0 > 1.860,0 2.790,0 2.700,0 1.620,0 1.440,0 1.530,0 > 930,0 930,0 900,0 420,0 840,0 720,0 > 1.080,0 480,0 870,0 2.970,0 2.700,0 2.520,0 > 2.790,0 1.590,0 930,0 390,0 720,0 2.340,0 > 930,0 930,0 900,0 390,0 390,0 390,0 > 1.860,0 750,0 780,0 720,0 780,0 720,0 > 3.720,0 3.720,0 3.600,0 840,0 840,0 840,0 > 780,0 420,0 60,0 390,0 390,0 420,0 > 2.700,0 2.790,0 2.790,0 720,0 780,0 360,0 > 930,0 930,0 900,0 420,0 420,0 360,0 > 2.700,0 3.720,0 3.720,0 1.170,0 1.170,0 540,0 > 1.860,0 930,0 900,0 420,0 420,0 720,0 > ..... > no total são 93 dados para cada coluna. > > a minha ideia é fazer o calculo das médias entre doset-9 e doset+9 > entre -6 e +6 e entre -3 e +3. > Ana > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 30, assunto 24 > ***************************************** ----- Fim de mensagem reenviada -----

6 8

Re: [R-br] Como extrair índices linha e coluna de um raster
by Alex 20 Jun '13

20 Jun '13

Elias, parece simples, mas não posso utilizar o centróide do pixel, pois haverá descontinuidade nas minhas amostras, conceitualmente está errado. As quadrículas (pixels) adjacentes são contínuas, mas entre um centróide e outro há descontinuidade, conceitualmente está inconsistente. Como mencionei, existe uma solução num software comercial FME (única que encontrei até então, nem mesmo no ArcGIS encontrei tal função) e o resultado da krigagem foi coerente, realmente estimei quadrículas (áreas) não amostradas. Att, Alex Santos

2 1

teste de wilcoxon
by alanarocha＠sapo.pt 20 Jun '13

20 Jun '13

Bom dia, alguém me pode dizer como resolvo este erro sem mudar os nomes das variáveis? setwd("f:/Ana_trabalho/Cinacalset") dados<-read.table("Cinacalset_DoseT.csv",header=TRUE,sep=";",dec=",",stringsAsFactors=FALSE) colnames(dados)<-c("DoseT-9", "DoseT-6", "DoseT-3", "DoseT+3", "DoseT+6", "DoseT+9") dados[,-c(1,2,5,6,7)] DoseT-3 DoseT+3 1 900,0 900,0 2 930,0 900,0 3 2.790,0 2.700,0 4 930,0 900,0 5 480,0 870,0 6 1.590,0 930,0 7 930,0 900,0 8 750,0 780,0 9 3.720,0 3.600,0 10 420,0 60,0 11 2.790,0 2.790,0 12 930,0 900,0 13 3.720,0 3.720,0 14 930,0 900,0 15 930,0 1.800,0 16 900,0 930,0 17 1.860,0 1.860,0 18 1.800,0 1.800,0 19 390,0 270,0 20 750,0 780,0 21 930,0 900,0 22 930,0 930,0 23 540,0 450,0 24 390,0 390,0 25 1.860,0 1.860,0 26 2.340,0 2.340,0 27 210,0 390,0 28 390,0 930,0 29 1.440,0 1.350,0 30 930,0 900,0 31 1.410,0 1.560,0 32 840,0 930,0 33 900,0 930,0 34 1.860,0 1.710,0 35 900,0 600,0 36 900,0 1.860,0 37 900,0 930,0 38 900,0 930,0 39 900,0 930,0 40 900,0 930,0 41 390,0 390,0 42 900,0 930,0 43 930,0 900,0 44 2.700,0 2.790,0 45 900,0 930,0 46 900,0 930,0 47 900,0 930,0 48 900,0 930,0 49 900,0 1.860,0 50 1.800,0 930,0 51 1.800,0 930,0 52 3.600,0 3.720,0 53 1.800,0 1.860,0 54 1.860,0 1.860,0 55 1.860,0 1.860,0 56 930,0 900,0 57 360,0 930,0 58 390,0 390,0 59 1.800,0 1.860,0 60 900,0 1.860,0 61 900,0 930,0 62 930,0 930,0 63 930,0 930,0 64 930,0 900,0 65 930,0 900,0 66 900,0 930,0 67 930,0 900,0 68 2.790,0 2.790,0 69 1.860,0 1.800,0 70 390,0 900,0 71 930,0 900,0 72 360,0 390,0 73 900,0 930,0 74 930,0 930,0 75 2.520,0 2.040,0 76 900,0 840,0 77 1.680,0 2.820,0 78 1.680,0 1.860,0 79 840,0 930,0 80 930,0 930,0 81 1.680,0 1.860,0 82 1.860,0 1.800,0 83 900,0 840,0 84 930,0 1.860,0 85 930,0 900,0 86 930,0 390,0 87 360,0 360,0 88 930,0 930,0 89 930,0 900,0 90 930,0 900,0 91 840,0 1.860,0 92 360,0 360,0 93 1.800,0 1.680,0 wilcox.test(dados[,-c(1,2,5,6,7)]$DoseT-3, + dados[,-c(1,2,5,6,7)]$DoseT+3, paired=TRUE) Error in wilcox.test.default(dados[, -c(1, 2, 5, 6, 7)]$DoseT - 3, dados[, : not enough (finite) 'x' observations Ana

2 1