- R-br - DINF UFPR

Operação com expand.grid
by ASANTOS 05 Nov '13

05 Nov '13

Boa tarde Pessoal, Estou com dificuldades de fazer um expand.grid() alterando-se as distância, como por exemplo se faço: # Ponto inicial p0<-c(1,1) c0 <- seq(p0[1],p0[1]+190,by=2) d0 <- seq(p0[2],p0[2]+47,by=3) d1 <- expand.grid(x=c0, y=d0) ## Coordenadas do grid espaçadas 3 x 2 plot(d1[,1],d1[,2]) #---------- Crio 16 pontos no eixo y e 32 pontos no eixo x, totalizando 1536 pontos, no entanto, eu gostaria de criar os mesmo 1536 pontos, mas gostaria de fazer da seguinte maneira: a cada 8 pontos na linha e 8 na coluna, distanciar em 4 unidades na linha e na coluna do próximo agrupamento de 8 pontos na linha e 8 na coluna para todos os pontos, formando ao final 24 agrupamentos de 8 pontos na linha por 8 na coluna separados de maneira equidistantes em 4 unidades. Alguém poderia me dar uma ajuda, pois fiz inúmeras tentativas sem sucesso? 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 ======================================================================

7 8

pseudo-inversa
by joão Rodrigo Castro 05 Nov '13

05 Nov '13

Pessoal, desejo fazer uma análise discriminante de Fisher e um dos passos é calcular a a inversa da matriz de covariância conjunta entre duas matrizes grandes. Essa matriz de covariância conjunta é singular e por isso não dá para calcular sua inversa, daí recebi uma orientação de usar a inversa de Moore-Penrose. Alguém sabe como efetuar esse cálculo no R? -- João Rodrigo de Castro Bolsista Programa de Educação Tutorial - PET Graduando em Meteorologia Universidade Federal de Pelotas

2 2

Dúvida Tutorial Geoestatística
by Helder Gramacho 04 Nov '13

04 Nov '13

Boa noite, sou iniciante na Geoestatística e no R, estava praticando numa espécie de tutorial situado no seguinte link: http://jurachumbo.blogspot.com.br/p/resultados_06.html Estava executando os comandos até que apareceu um erro que não estou conseguindo entender, seguem os comandos e o erro: data(jura) jura.grid=juragrid.dat coordinates(jura.grid)=~Xloc+Yloc gridded(jura.grid) = TRUE m<-vgm(0.12645599,"Exp", 0.1593586, 0.04550649) x<-krige(log(Pb)~1,jura, jura.grid, model=m) Error in krige(log(Pb) ~ 1, jura, jura.grid, model = m) : erro na avaliação do argumento 'locations' na seleção do método para a função 'krige': 'Error: object 'jura' not found ' Alguém saberia me explicar qual é o problema? Desde já agradeço, Hélder Gramacho dos Santos Engenheiro Agrônomo Mestrando Ciênc. Geod. Tec. Geoinformação-UFPE *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 2

Duvida - Ellipse
by Carlos Pereira 04 Nov '13

04 Nov '13

Preciso construir elipses de confiança para um conjunto de dados bidimensionais. Mas não quero assumir normalidade. Conheço as posteriores condicionais dos parâmetros, mas meu interesse e em uma variável defendida como a multiplicação da raiz quadrada dos valores singulares pelas coordenadas dos autovetores. Se eu fizer esta multiplicação nas cadeias consigo gerar uma amostra a posteriori, mas não conheço a distribuição conjunta destas variáveis. Desta forma tenho que construir elipses baseadas em agrupamento de pontos. E não sei como construí- las. Logo se alguém souber como posso fazer peço que me ajude. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MGEmail: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

3 2

Codificação de arquivo .nnw
by Andre Oliveira 04 Nov '13

04 Nov '13

Boa noite, ha 3 anos contrui uma apostila usando o Latex + Sweave .rnw e foi criado no Windows com TinnR. Resolvi atualisar a apostila e compilar o arquivo .rnw, porém, percebi que meu arquivo .rnw está com toda ascentuação comprometida. Onde tem ascentos, til tudo com ponto de interrogação. Recorro a vocês do grupo para tentar resolver este problema. Como mudei para Linux heguei abrir este arquivo no Linux com Rstudio. Alguma dica? obrigado André Oliveira Souza

4 4

retidada de tendencia sem UTM
by Hélio Gallo Rocha 04 Nov '13

04 Nov '13

Caros da lista, Estudando outros dados, não houve problema na validação cruzada para retirada de tendencia polinomial de primeiro grau, achei que fosse pela quantidade de pontos, mais de 300. Como minhas coordenadas estão em UTM, subtrai todas elas pelo menor valor. Rodei o CRM... sem problemas. Para quem está em fase de estudos da geoestatistica, ainda faltam pontos de dúvida, favor darem suas opiniões, listo aqui minhas dúvidas e acrescentem outras: 1. Sempre temos de retirar tendencia? 2. Poderia retirar tendencia com uma covariavel fenológica da planta? penso que sim, já que a retirada polinomial de primeiro grau utiliza as coordenadas. 3. Alguém teria algum artigo para embasar este tipo de situação? Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 3

Inserir desvio padrão em gráfico de barras
by Rodrigo Muniz 04 Nov '13

04 Nov '13

Prezados membros da lista gostaria de pedir o auxílio de vocês. A partir dos seguintes comando: dados=cbind(c("Mata Nativa", "Citrus 37", "Citrus 37", "Citrus 24", "Citrus 18", "Citrus 31","Pasto 37"), tapply(arenosa[1:21, "ECA"], arenosa[1:21, "Trat"], mean), tapply(arenosa[1:21, "ECB"], arenosa[1:21, "Trat"], mean), id=1:nrow(tapply(arenosa[1:21, "ECA"], arenosa[1:21, "Trat"], mean))) colnames(dados)<-c("Classes", "0-10", "10-20", "id") dados= as.data.frame(dados) barchart(value~Classes, groups=variable, data=melt(dados, measure.var=2:3), auto.key=list(columns=2), horizon=FALSE, ylab="Estoque Carbono", main="Área 3") Eu gero a seguinte tabela de dados com o respectivo gráfico de barra. Classes 0-10 10-20 id 301 Mata Nativa 10.04 7.6 1 312 Citrus 37 12.36 23.12 2 322 Citrus 37 15.32 27.58 3 332 Citrus 24 4.38 8.85 4 342 Citrus 18 10.43 21.57 5 352 Citrus 31 9.19 18.03 6 364 Pasto 37 6.69 12.53 7 Eu gostaria de inserir neste gráfico de barras o desvio padrão de cada barra. Como eu poderia fazer isso? Segue a tabela com os respectivos desvios padrão. 0-10 10-20 Mata Nativa 3.4250012 0.7715569 Citrus 37 2.0617064 4.1486182 Citrus 37 2.8501287 4.7112737 Citrus 24 0.1931321 0.3666515 Citrus 18 1.7760349 4.9343726 Citrus 31 0.6730775 1.6246333 Pasto 37 2.4274060 2.8266116 Desde já agradeço. -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

2 1

Dúvida no Mongo DB
by Savano Pereira 04 Nov '13

04 Nov '13

Prezados, Estou acessando um MongoDB remotamente utilizando o R. Gero uma consulta, porém, quando crio um data frame com a consulta ele está vazio, alguém sabe me dizer onde está o problema?. Estou começando a me inteirar sobre o MongoDB. Envio o log do programa... Desde já agradeço. > library(rmongodb) > > # parâmetros > host <- "xxxx" > username <- "xxx" > password <- "xxxx" > db <- "censo2010" > > #conectando a base > teste <- mongo.create(db="censo2010", host="xxxx", username="xxxx",password="xxxx") > > #Obtendo a lista de bancos no MongoDBserver > print(mongo.get.databases(teste)) NULL > > #obtendo a lista de collection > print(mongo.get.database.collections(teste, db)) [1] "censo2010.censo2010" > > #Returns a fresh mongo.bson.buffer object ready to have data appended onto it > buf <- mongo.bson.buffer.create() > > #Start a subobject within a mongo.bson.buffer > mongo.bson.buffer.start.object(buf, 'V0001') [1] TRUE > #Append a name/value pair into a mongo.bson.buffer > mongo.bson.buffer.append(buf, '$eq', '11') [1] TRUE > #Finish a subobject or array within a mongo.bson.buffer > mongo.bson.buffer.finish.object(buf) [1] TRUE > > #Convert a mongo.bson.buffer object to a mongo.bson object > query <- mongo.bson.from.buffer(buf) > > #Find records in a collection > cur <- mongo.find(teste, 'censo2010.censo2010', query = query) > > V0001 <- NULL > while (mongo.cursor.next(cur)) { + value <- mongo.cursor.value(cur) + V0001 <- rbind(V0001, mongo.bson.value(value, 'V0001')) + } > data2 <- data.frame(V0001=V0001) # AQUI O DATAFRAME ESTÁ VAZIO -- Savano Pereira TEL: 21 8771-9063

1 0

R Markdown através do Rcmdr
by pjobenil＠gmail.com 04 Nov '13

04 Nov '13

Boa tarde pessoal Passei um trabalho simples de análise de dados para uns alunos de engenharia e espero que eles usem o R Markdown através do Rcmdr para gerar um relatório com as análise para me entregar. Porém o ambiente não aceita caracteres latinos (ã, õ, í, etc) e o texto em html sai com uma série de cacos. Alguém conhece alguma forma de solucionar isso? Agradeço antecipadamente, Atenciosamente Jobenil Júnior Enviado do Email do Windows

4 6

Interface em R
by Tiago Souza Marçal 04 Nov '13

04 Nov '13

Boa tarde pessoal, Andei dando uma olhada no pacote tclk e suas funcionalidades na criação de interfaces. No entanto, gostaria de saber se existe algum plugin que possa ser instalado no computador para que a interface pudesse ser construida via netbeans. 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 #################################################################

3 3