R help - Mar 2008 - howto find corresponding values in datasource?

Hi,


i am sure, that this is a noob-question, but i have searched for 
hours without any good result. 

I want to draw a vertical line through the maximum of the first derivation.

Here is a small example. 

--8<-----[mydata.csv]--------------------------------
H    M
1    15
2    22
3    23
4    17
5    10
--8<-----[myquestion.R]------------------------------
mydata <- read.table("mydata.csv", header=TRUE, sep="\t")

attach(mydata)
# make a smooth fit through the points and calculate the first derivation d/dx
myspl <- smooth.spline(H, M, all.knots = FALSE, nknots = 10, spar=0.5) 
mydiff <- predict(myspl, 0:max(H), deriv=1)


# find the maximum peak max(y)
y1=y2= max(mydiff$y)

# now i want to plot a line through the maximum of the derivation
# how can i get the x-coordinate?
# segments(x1,y1,x2,y2)

detach()
--8<-------------------------------------

Thank you,

-- 
Jonas Stein

The x coordinate of the max y value:

x[which.max(mydiff$y)]




Jonas Stein wrote:
> 
> Hi,
> 
> 
> i am sure, that this is a noob-question, but i have searched for 
> hours without any good result. 
> 
> I want to draw a vertical line through the maximum of the first
> derivation.
> 
> Here is a small example. 
> 
> --8<-----[mydata.csv]--------------------------------
> H    M
> 1    15
> 2    22
> 3    23
> 4    17
> 5    10
> --8<-----[myquestion.R]------------------------------
> mydata <- read.table("mydata.csv", header=TRUE,
sep="\t")
> 
> attach(mydata)
> # make a smooth fit through the points and calculate the first derivation
> d/dx
> myspl <- smooth.spline(H, M, all.knots = FALSE, nknots = 10, spar=0.5) 
> mydiff <- predict(myspl, 0:max(H), deriv=1)
> 
> 
> # find the maximum peak max(y)
> y1=y2= max(mydiff$y)
> 
> # now i want to plot a line through the maximum of the derivation
> # how can i get the x-coordinate?
> # segments(x1,y1,x2,y2)
> 
> detach()
> --8<-------------------------------------
> 
> Thank you,
> 
> -- 
> Jonas Stein
> 
> ______________________________________________
> R-help at r-project.org 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.
> 
> 
-- 
View this message in context:
http://www.nabble.com/howto-find-corresponding-values-in-datasource--tp16016130p16023738.html
Sent from the R help mailing list archive at Nabble.com.