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]]
Hola¡ Ya se cual fue mi error,no fue R... Gracias, un saludo 2014-07-10 17:43 GMT+02:00 Marta valdes lopez <martavaldes85@gmail.com>:> 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]]