Message-ID: <DC3FB55EE7FEFD409A6FC1EA00D50D0B0F87AF35@NDJSMBX203.ndc.nasa.gov>
Date: 2017-09-25T19:28:28Z
From: Evans, Richard K. (GRC-H000)
Subject: bowed linear approximations
Hello,
Please run the following code snippet and note the resulting plot:
x <- c(10, 50)
y <- c(0.9444483, 0.7680123)
plot(x,y,type="b",log="x")
for(i in 1:50){
xx <- exp(runif(1,log(min(x)),log(max(x)) ))
yy <- approx(x,y,xout=xx, method = "linear")
points(xx,yy$y)
}
notice the "log=x" plot parameter and the resulting "bow" in the linear approximation.
This makes sense when I realized that the plot command is first making the plot and then drawing straight lines between points on a log plot AFTER the plot is generated and that that's why the line is straight. I get that.
.. and it also makes sense that the bowed points are a result of the linear approximations being made BEFORE plotted in a logarithmic plot. I get that..
My goal is to make approximations that lie on the line produced on the plot as shown, so I realize that what I want to do is NOT linear approximations, but maybe "log" approximations?
However, the approximation methods are only "linear" and "constant" .. there isn't a "log" method to approximate with.
So can anyone tell me how to fix the code such that he approximated points DO lie on the line as plotted with the "log=x" plot parameter?
Oh, and they have to be uniformly distributed along the Log=x axis.. I think that's the tricky part.
Any help and/or insight would be greatly appreciated.
Thank you!
-Rich
[[alternative HTML version deleted]]