Hola !!
Estoy utilizando este script hace tiempo y hasta ahora he ido encontrando
errores que he ido solucionando, pero hoy al utilizarlo otra vez me da el
error de que al crearme las tablas que pido me borra parte de los nombres
de los puertos (Harbour), y ya comprobe en el original y aqui estan todos
bien.Lo que no entiendo es que el script es el mismo y no se porque ahora
me crea ese error.
setwd("C:/Users/mvaldes/Desktop/DB")
library(chron)
library(xlsx)
#Read output of previous routine
#Give FILENAME here
filename<-"univerest_50.csv"
DBxy<-read.csv(filename, sep=";",header=TRUE,dec=",")
DBx<-DBxy
names(DBx)<-c("Boat","DateTime","TimeDiff", "Latitude", "Longitude",
"Course", "Speed", "distNm", "calcSpeed", "calcCourse", "distHb",
"Harbour", "idTrip","vmsAngle", "calcAngle", "vmsLeg", "calcLeg",
"Trip_vmsLeg", "Trip_calcLeg", "lengthTrip", "lengthTrip_vmsLeg",
"lengthTrip_calcLeg","Time", "Date")
#Formatting date and time variables
DBx$Date<-strptime(DBx$Date, "%d-%m-%Y", tz="GMT") ###Monicap
use the same ; other use Y
DBx$Year<-as.POSIXlt(DBx$Date)$year+1900
if(filename!="monicap_50.csv") {DBx$Time<-paste(DBx$Time, ":00",
sep="")} #NOT necessary for Monicap and Univerest_50
DBx$Time<-times(DBx$Time) #Works for Monicap AND UNIVEREST_50 ONLY
DBx$Boat<-gsub("^\\s+|\\s+$", "", DBx$Boat)
#Read file with boat codes and gears
codeBoats<- read.csv("CODES_2002-2010New.csv",
sep=",",header=TRUE) #Laptop
codeBoats$CODIGO<-gsub("^\\s+|\\s+$", "", codeBoats$CODIGO)
#Assigning a Fishing license based on Boat and Year
DBx$gear<-codeBoats$Lic[match(paste(DBx$Boat,DBx$Year),
paste(codeBoats$CODIGO,codeBoats$Year))]
z0<-length(DBx$gear)
######PoleLine analyses / MAY 2013 #######Changes are: harbour with speed
rule (<rest)
#Define GEAR HERE based on list above. Analysing BLL boats
DBx <- DBx[DBx$gear %in% c("PoleLine"),]
z1<-length(DBx$gear)
#defining speed and distance limits
speedFishing<-2.0 #less or equal to, in knots (NEW MAY 2013
reduced from 5.3)
speedHarb<-1.0 #new speed harbour
distHbRule<-3.0 #greater than, in nautical miles nm
minTime<-times(c("05:59:59")) #usual beginning of fishing
operations
maxTime<- times(c("20:59:59")) #usual finishing of fishing
operations
#Selecting Harbour
DBharbour<- na.omit(DBx[DBx$distHb<=distHbRule &
DBx$calcSpeed<=speedHarb,])
DBharbour$State<-"Harbour"
z2<-length(DBharbour$State)
write.table(DBharbour, file =
paste("C:/Users/mvaldes/Desktop/DB/results_",filename, sep=""), append = F,
quote = F, sep = ",", dec = ".",
row.names = F, col.names = T)
#Selecting Steaming
DBsteaming<- na.omit(DBx[DBx$calcSpeed>speedFishing |
(DBx$distHb<=distHbRule & DBx$calcSpeed>speedHarb),])
DBsteaming$State<- "Steaming"
DBsteaming$Harbour<-""
z4<-length(DBsteaming$State)
write.table(DBsteaming, file =
paste("C:/Users/mvaldes/Desktop/DB/results_",filename, sep=""), append = T,
quote = F, sep = ",", dec = ".",
row.names = F, col.names = F)
#Selecting Fishing
DBfishing<- na.omit(DBx[DBx$calcSpeed<=speedFishing &
DBx$distHb>distHbRule & (DBx$Time>minTime & DBx$Time<=maxTime),])
DBfishing$State<-"Fishing"
DBfishing$Harbour<-""
z5<-length(DBfishing$State)
write.table(DBfishing, file =
paste("C:/Users/mvaldes/Desktop/DB/results_",filename, sep=""), append = T,
quote = F, sep = ",", dec = ".",
row.names = F, col.names = F)
#Selecting nigth
DBnight<- na.omit(DBx[(DBx$calcSpeed<=speedFishing &
DBx$distHb>distHbRule &(DBx$Time<=minTime | DBx$Time>maxTime)),])
DBnight$State<-"Night"
DBnight$Harbour<-""
z6<-length(DBnight$State)
write.table(DBnight, file =
paste("C:/Users/mvaldes/Desktop/DB/results_",filename, sep=""), append = T,
quote = F, sep = ",", dec = ".", row.names = F, col.names = F)
Si alguien lo ve y me puede decir lo agradeceria, o simplemente R me vacila
ya no se que pensar...
Muchas gracias por vuestro tiempo, un saludo
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
