Displaying 1 result from an estimated 1 matches for "ddt2".
Did you mean:
dat2
2012 Mar 19
3
Issue with asin()
...) / 2)
print((Tseuil-meanT)/amp)
alpha <- asin((Tseuil - meanT) / amp)
DD_ifelse3 <- ((1 / pi) * ((meanT - Tseuil) * ((pi/2) - alpha)) +
amp*cos(alpha))
DD <- ifelse(cond1, 0, ifelse(cond2, (meanT - Tseuil), DD_ifelse3))
}
### method (Tmin + Tmax) / 2
else if(method == "DDt2"){
cond1 <- (meanT > Tseuil)
DD <- ifelse(cond1,(meanT - Tseuil),0)
}
else{
stop("\nMethod name is invalid.\nMethods available = DDsin (sinus) or DDt2
(mean)\n")
}
return(DD)
}
BUT! When I try to process random data:
library(reshape2)
library(plyr)
station &l...