Displaying 1 result from an estimated 1 matches for "caudalesb".
2013 May 31
1
consulta
...den]
zeros<-c()
for(i in 1:(length(tiempo)-1)){
if((valor[i]-qmed)*(valor[i+1]-qmed)<0){
zeros<-c(zeros,tiempo[i]+(valor[i]-qmed)*(tiempo[i+1]-tiempo[i])/(valor[i+1]
-valor[i]))
}
}
zeros
}
tiemposb<-c(Tiempos,zeros(Tiempos,Caudales,qmed))
Caudalesb<-c(Caudales,rep(0,length(tiemposb)-length(Caudales)))
orden<-order(tiemposb)
tiemposb<-tiemposb[orden]
Caudalesb<-Caudalesb[orden]
positivos<-ifelse(Caudalesb<=0,0,Caudalesb)
negativos<-ifelse(Caudalesb>=0,0,Caudalesb)
plot(tiemposb,Caudalesb,type="l")
points(t...