Dicas rápidas para R

POST EM CONSTANTE ATUALIZAÇÃO
PARA VER MAIS DICAS, ASSUNTOS E COMANDOS
BUSQUE NA BARRA ACIMA, À ESQUERDA

O ambiente de software estatístico R é muito difundido em várias áreas de ciêncie, indo de biologia até geofísica. A linguagem a princípio é bem dificil de se lidar por ser basicamente toda vetorial, mas eu gosto sempre de ter uma "colinha" para me ajudar.

Coloquei abaixo algumas dessas dicas, algumas bobas e outras que vivo esquecendo...

EXECUTAR UM SCRIPT NO PROMPT DO R
source(paste(dir,"/batchGRADE.R",sep=""))

QUEBRA UM VETOR CHARACTER EM VARIAS COMPONENTES
strsplit("ola tudo bem"," ")
resultado: um vetor de 3 posições: "ola", "tudo" e "bem"

PARA ACESSAR APENAS UMA PARTE DE UM MAPA E REESCREVER O SHP
shape0=c("0.shp")
library(maps)
library(maptools)
mapa0=readShapeLines(shape0)
writeLinesShape(mapa0[mapa0[[3]]=="RJ",],"testeRJ")

#MOSTRA APENAS AS RODOVIAS DO SHAPE
plot(mapa2[mapa2[[1]]=="Rodovia",])

#MOSTRA APENAS UM SUB-OBJETO DO SHAPE
plot((mapa2[1,]@lines[[1]])@Lines[[1]]@coords,type="l")

PARA RODAR UM SCRIPT NO R (tabajara version)
#!/bin/ksh
/usr/bin/time R --no-save --no-restore --no-site-file --no-init-file --no-environ --no-readline < script.R

# Usando o X Virtual Frame Buffer como servidor X para rodar o R em crontab
/usr/bin/Xvfb :7 -screen 0 1152x900x24&
export DISPLAY=":7.0"

GREP TABAJARA EM R
as.real(read.table("model.ctl",skip=5,nrows=1)[2])
(o exemplo pega o valor de tdef num ctl)

APAGA TODOS OS OBJETOS
rm(list = ls())

TAMANHO DO OBJETO NA MEMÓRIA
object.size(temp)

MEMÓRIA
memory.size()
gc()

EXEMPLO DE IF VETORIAL
sqrt(ifelse(x > 0, x, NA))

SEQUÊNCIAS E CLASSIFICAÇÃO
b1=rep(a1,length(a2))
b2=rep(a2,length(a1))
b2=sort(b2)
lon=seq(-45,,length.out=5)
llon=length(lon)
lon=rep(lon,llon)
dim(lon)=c(llon,llon)
lon=t(lon)

VIEWPORTS E ARROWS
plot.new()
plot(lon, lat, type="n",xlab="latitude", ylab="longitude",main="China
Sea Wind Speed/Direction and Temperature")
vps <- baseViewports()
par(new=TRUE)
pushViewport(vps$inner, vps$figure, vps$plot)
length <- 0.5 # "cm"
x1 <- unit(lon, "native") - unit(0.5*length*bu, "cm")
y1 <- unit(lat, "native") - unit(0.5*length*bv, "cm")
x2 <- unit(lon, "native") + unit(0.5*length*bu, "cm")
y2 <- unit(lat, "native") + unit(0.5*length*bv, "cm")
grid.arrows(grob=grid.segments(x1, y1, x2, y2),length=unit(1.5, "mm"),gp=gpar(col="red"))
popViewport(3)
plot(mapa_estado,lwd=2,xlim=c(a2[1],a2[5]),ylim=c(a1[1],a1[8]))
#plot(mapa_estado,lwd=2,xlim=c(-25,-18),ylim=c(-45,-41))



POST EM CONSTANTE ATUALIZAÇÃO
PARA VER MAIS DICAS, ASSUNTOS E COMANDOS
BUSQUE NA BARRA ACIMA, À ESQUERDA

Nenhum comentário: