Re: [R] how to fit gauss beam?

From: Thomas Hunger <hto_at_arcor.de>
Date: Mon 04 Sep 2006 - 21:46:48 GMT

Reply to self:
> I am having a hard time fitting a gauss beam using R. In
> gnutplot I did something like
>
> $ w(z) = w0 * sqrt(1+(z/z0)**2)
> $ fit w(z) 'before_eom.txt' using 1:2 via w0, z0
>
> to obtain w0 and z0. Now I want to do the same in R. I
> tried a linear model like this (r = radius, z =

This works fine:

data <- read.table ("nach_eom.tab", header=T)

M <- 1.1
nls (major ~ M*w0*sqrt(1+(d/z0)^2),

     data = data, 
     start = list(w0 = 460, z0=0.02))

______________________________________________
R-help@stat.math.ethz.ch 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. Received on Tue Sep 05 07:50:14 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Thu 07 Sep 2006 - 07:51:18 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.