Hello,
I am deriving near real-time liner relationships based on 5-min
precipitation data, sometimes the non-qced data result in a slope of NA. I
am trying to read the coefficient (in this example x) to see if it is equal
to NA, if it is equal to NA assign it a value of 1. I am having trouble
with the if statement not recognizing the coefficient or "NA" value in the
test.
Any thoughts, I supplied a really basic example below.
Thank you,
Doug
#####################
###### Example ######
#####################
> x=c(1,1,1)
> y=c(1,1,1)
> fit <- lm(y~x)
> fit
Call:
lm(formula = y ~ x)
Coefficients:
(Intercept) x
1 NA
> coef <- coef(fit)
> fit$coef[[2]]
[1] NA
> if("fit$coef[[2]]" == "NA") {.cw = 1}
***************************************
Douglas M. Hultstrand
Ph.D. Candidate Earth Sciences
Watershed Science Program
Colorado State University
http://www.cnr.colostate.edu/~dmhultst/
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
