j<-c("Gener","Febrer","Març","Abril","Maig","Juny","Juliol","Agost","Setembre","Octubre","Novembre","Desembre") i<-c("1","2","3","4","5","6","7","8","9","10","11","12","13","14","15","16","17","18","19","20","21","22","23","24","25","26","27","28","29","30","31") resultat=matrix(ncol=12,nrow=5,dimnames=list(c("p-valor","Assal. darrer dia","Int. conf. inf.","Int. conf. sup.","Situació dada darrer dia"),j)) Mataro.diari<-matrix(c(1191,1193,1193,1202,1201,1195,1202,1200,1195,1195,1205,1205,1203,1201,1200,1187,1183,1187,1190,1190,1188,1192,1190,1190,1195,1196,1196,1198,1199,1190,1190,1190,1190,1191,1191,1190,1187,1187,1194,1191,1193,1195,1199,1198,1197,1201,1200,1201,1203,1205,1202,1200,1207,1208,1208,1213,1215,1208,1206,NA,NA,NA,1207,1206,1209,1209,1207,1205,1204,1210,1213,1211,1211,1212,1205,1205,1210,1211,1215,1216,1214,1212,1208,1214,1215,1218,1216,1217,1214,1214,1214,1216,1219,1211,1204,1204,1205,1204,1207,1209,1212,1212,1208,1208,1214,1213,1218,1218,1220,1219,1219,1220,1223,1224,1223,1228,1221,1218,1220,1223,1222,1221,1221,NA,1213,1212,1221,1220,1220,1220,1222,1220,1218,1219,1220,1221,1220,1220,1223,1221,1224,1227,1228,1228,1234,1224,1223,1223,1223,1222,1221,1223,1213,1210,1217,1218,1220,1222,1222,1221,1220,1225,1225,1226,1229,1231,1228,1230,1239,1240,1245,1243,1242,1238,1236,1238,1233,1215,1214,1216,1214,1207,1231,1233,1233,NA,1227,1233,1225,1226,1231,1229,1227,1231,1232,1227,1222,1228,1222,1227,1245,1247,1241,1235,1242,1246,1250,1248,1249,1235,1236,1237,1237,1237,1238,1239,1228,1218,1222,1223,1227,1228,1226,1222,1224,1225,1229,1226,1228,1232,1223,1224,1227,1224,1218,1220,1222,1213,1214,1215,1219,1218,1218,1218,1211,1210,1216,1208,1199,1202,1203,1197,1194,1197,1215,1220,1222,1220,1212,1210,1217,1217,1219,1222,1222,1215,1215,1219,1220,1216,1218,1216,1216,1214,1214,1209,1212,1212,NA,1201,1200,1204,1211,1211,1213,1212,1218,1219,1216,1223,1216,1217,1216,1220,1216,1211,1218,1220,1221,1219,1218,1221,1220,1225,1228,1228,1224,1221,1215,1217,1205,1217,1217,1216,1218,1217,1216,1223,1222,1220,1220,1221,1217,1216,1224,1228,1234,1232,1229,1226,1222,1227,1231,1231,1233,1233,1226,1229,1234,1233,NA,1237,1237,1235,1230,1228,1228,1235,1229,1235,1236,1229,1228,1241,1243,1236,1240,1243,1236,1236,1247,1248,1250,1251,1247,1238,1237,1251,1252,1252,1250,1244),byrow=FALSE,ncol=12,dimnames=list(i,j)) last.day<-c(1190,1206,1219,1221,1217,1233,1228,1208,1212,1217,1233,1244) for (t in 1:12) { resultat[1,t]<-round(wilcox.test(Mataro.diari[,t],mu=last.day[t],conf.int=TRUE,na.rm=TRUE)$p.value,digit=4); resultat[2,t]<-last.day[t]; resultat[3,t]<-round(wilcox.test(Mataro.diari[,t],mu=last.day[t],conf.int=TRUE,na.rm=TRUE)$conf.int[1],digit=1); resultat[4,t]<-round(wilcox.test(Mataro.diari[,t],mu=last.day[t],conf.int=TRUE,na.rm=TRUE)$conf.int[2],digit=1); resultat[5,t]<-ifelse(last.day[t]resultat[4,t],"Per sobre la mitjana","En la mitjana")) } resultat