OperacionesSuma +, Resta -, Multiplicación * y División /Potenciación ^ y raíz cuadrada sqrtDivisión entera %/% y modulo %%Logaritmos log,log2,log10 y exponencial exp Trigonométricas sin, cos, tan, ....Operadores lógicos <,>,<=,>=,==,!=Estadísticas sum,mean, median,sd,...
Operaciones con matricesA %*% B : producto de matricest(A) : transpuesta de la matriz Asolve(A,b) : solución del sistema de ecuaciones Ax=b.solve(A) : inversa de la matriz Adiag(A) : matriz diagonal (A es una matriz)det(A): Determinante de la matriz A
Asignaciónvalores<-c(1,2,3,4)suma<-sum(valores)promedio<-suma/4
MatricesA<-c(1,2,3,4)dim(A)<-c(2,2)B<-c(5,6,7,8)dim(B)<-c(2,2)
Tipos de Archivo “nativos” de R
Datos (.RData)
Scripts (.r)
Importar datosLibreria (foreign)
DBFread.dbf(file)
Stataread.dta(file)
Epi inforead.epiinfo(file)
Minitabread.mtp(file)
SPSSread.spss(file)
SASread.ssd(libname, sectionnames, tmpXport=tempfile(), tmpProgLoc=tempfile(), sascmd="sas")
Libreria(xlsx)
read.xlsx(file)
read.csv(file)
Distribuciones de ProbabilidadDistribucion Parametros en RNormal m, sd normt-Student n tChi-Cuadrado n chisqF n,m fExponencial l expUniforme min,max
unifBeta p,q betaCauchy t,s cauchyLogística t,s logisLognormal m,sd lnormGamma p,s gammaWeibull p,s weibullGumbel p,s gumbel
Distribucion Parametros en RBinomial n,p binomBinomial neg. n,p nbinomGeometrica p geomHipergeometrica m, n, k hyperPoisson l pois
CONTINUAS
DISCRETAS
Función de Densidad dFunción de Distribución pQuantil qGenerar muestra aleatoria r
Nemotécnicos para funciones de distribución
Sintaxis en RNemo+Distribucion(Valor,parametros)
D1. Distribución Normal. Cambio de parámetros
d1_DISTRIBUCION_NORMAL.r
Suponiendo que Z ~ N(0,1), calcule la probabilidad de que: (haga gráficas de lo que se pide)a) Z sea menor que 1.48b) Z sea mayor que 1.90c) Z tenga un valor entre 1.48 y 1.90d) Z tenga un valor entre -1.48 y 1.90¿Cuál es el valor de Z sie) La probabilidad de observar valores mayores que él es de 0.60?f) La probabilidad de observar valores menores que él es de 0.10?g) ¿Cuáles son los valores de los cuantiles de tal manera que la probabilidad de que Z tome valores entre estos dos cuantiles es de 0.90? (distribuidos simétricamente alrededor del 0)
Ejercicio 1Distribuciones
Ej1_distribuciones.r
> pnorm(1.48,0,1)[1] 0.9305634
> pnorm(1.9,0,1,lower.tail=F)[1] 0.02871656
> pnorm(1.9,0,1)-pnorm(1.48,0,1)[1] 0.04072006
> pnorm(1.9,0,1)-pnorm(-1.48,0,1)[1] 0.9018468
> qnorm(0.4,0,1)[1] -0.2533471 > qnorm(0.1,0,1)
[1] -1.281552
> qnorm(0.05,0,1)[1] -1.644854> qnorm(0.95,0,1)[1] 1.644854
Ejercicio 2Índice IMECA zona noreste del DF
8-9 Buena9-10 Buena10-11 Buena11-12 Regular12-13 Regular13-14 Mala14-15 Mala15-16 Regular16-17 Regular17-18 Regular
imeca.txt
Fuente:Dir. De Monitoreo AtmosféricoSecretaria de Medio Ambiente GDF5 de junio de 2012
setwd("D:/curso_PUMA_IIMAS")imeca<-read.table("imeca.txt",sep=",",quote="\"",head=T)attach(imeca)valores<-table(Calidad)barplot(valores,col=c("blue","yellowgreen","brown"),main="Calidad del Aire zona noreste del DF")pctje<-valores/sum(valores)
Ejercicio 3Reforestación urbana en el DF,1999
Delegacion Reforestacion SuperficeAlvaro Obregon 170 96.17Azcapozalco 104 33.66Benito Juarez 229 26.33Coyoacan 70 54.4Cuajimalpa 15 74.58Cuautemoc 357 32.4Gustavo A. Madero 183 94.07Iztacalco 321 23.3Iztapalapa 117 117Magdalena Contreras 106 74.58Miguel Hidalgo 313 46.99Milpa Alta 14 228.41Tlahuac 63 85.34Tlalpan 171 340.07Venustiano Carranza 42 33.4Xochimilco 252 118
datos_reforestacion_df.xlsx
Fuente: Secretaria del Medio Ambiente, GDF, 2000
summary(Reforestacion) Min. 1st Qu. Median Mean 3rd Qu. Max. 14.00 68.25 143.50 157.90 234.80 357.00
sd(Reforestacion)110.671cv0.70072633178926es27.6677412567416
Ejercicio 4Composición de Suelos
Ca K Mg Na1 0.4626 0.1172 0.1119 0.02532 0.0148 0.0616 0.0563 0.01063 0.0184 0.08 0.0891 0.01054 0.0309 0.0698 0.079 0.01265 0.0431 0.0475 0.0569 0.04466 1.6493 0.1737 0.4764 0.03247 0.7555 0.1562 0.383 0.03068 0.1276 0.1844 0.1857 0.01819 9.1687 0.1992 1.3042 0.0218
10 0.2919 0.1826 0.388 0.036611 0.8167 0.1288 0.0012 0.018712 1.8544 0.1219 0.2427 0.013813 0.4856 0.1092 0.0119 0.028914 4.4653 0.1212 0.5535 0.028115 2.9269 0.1498 0.1131 0.0295
datos_suelo.xlsx
mean(datos)Ca K Mg Na 1.5407800 0.1268733 0.2701933 0.0241400
sd(datos) Ca K Mg Na 2.46275245 0.04758910 0.33550117 0.01011807
cov(datos) Ca K Mg Na
Ca 6.065149626 0.057156272 0.7316968163 0.0008993380 K 0.057156272 0.002264722 0.0100180234 0.0000969790 Mg 0.731696816 0.010018023 0.1125610321 0.0004152596 Na 0.000899338 0.000096979 0.0004152596 0.0001023754
cor(datos) Ca K Mg Na
Ca 1.00000000 0.4876808 0.8855567 0.03609145 K 0.48768082 1.0000000 0.6274521 0.20140599 Mg 0.88555670 0.6274521 1.0000000 0.12232853 Na 0.03609145 0.2014060 0.1223285 1.00000000
Ejercicio 5Emisiones totales de fuentes móviles en la Zona Metropolitana de la Ciudad de México 1998
Fuente: Secretaria del Medio Ambiente, GDF, 1999
datos_moviles.xlsx
Fuente HC CO NOx PM10 SO2
particulares 75676 873287 35874 759 2374
taxis 18798 210368 8642 183 572
microbuses 17837 194698 6150 68 271
pick-up 2958 33974 1249 18 73
camion_carga 59882 653643 20646 229 909
autobuses 2416 5977 7330 763 161
carga_mayor_2e 20679 51153 62733 6530 1163
mean(datos) HC CO Nox PM10 SO228320.857 289014.286 20374.857 1221.429 789.000
sd(datos)HC CO Nox PM10 SO228319.4441 339397.2961 22024.5409 2361.0313 805.0832
cov(datos)HC CO Nox PM10 SO2
HC 801990915 9422074226 250450010 -4901465.8 19652056.3 CO 9422074226 115190524605 1652276129 -211270863.6 217857652.8Nox 250450010 1652276129 485080403 45835453.2 12040235.7 PM10 -4901466 -211270864 45835453 5574469.0 517621.3 SO2 19652056 217857653 12040236 517621.3 648159.0
cor(datos)HC CO NOx PM10 SO2
HC 1.0000000 0.9802874 0.4015406 -0.0733060 0.8619508 CO 0.9802874 1.0000000 0.2210382 -0.2636510 0.7973035 Nox 0.4015406 0.2210382 1.0000000 0.8814403 0.6790275 PM10 -0.0733060 -0.2636510 0.8814403 1.0000000 0.2723138 SO2 0.8619508 0.7973035 0.6790275 0.2723138 1.0000000