Fit power curve that goes through specified coordinates in r -


i have following x-y pairs:

x <- seq(.3, .7, by=0.01) y <- .5 * (x**-2.5) + runif(n = 41, min = -0.25, max = 0.25)

i find parameters , b of curve form y-hat = * (x ** b), such sum of squared error between y , y-hat minimized.

furthermore, curve must pass through coordinates (0.7, 1.0). need solve a, since b = -log(a)/log(0.7).

i can in excel solver, writing own algorithm solve in r beyond current ability. i'm guessing built in solution exists anyway. looked powerlaw package there doesn't appear way force curve through specified point.

set.seed(101) x <- seq(.3, .7, by=0.01) y <- .5 * (x**-2.5) + runif(n = 41, min = -0.25, max = 0.25)

you need reasonable starting point: can 1 log-log regression.

m1 <- lm(log(y)~log(x))  ## hope data positive a_est  <- exp(coef(m1)[1])

now you're ready estimate:

n1 <- nls(y ~ a*x^(-log(a)/log(0.7)),start=list(a=a_est)) <- coef(n1) b <- -log(a)/log(0.7) plot(x,y) curve(a*x^b,add=true,col="red")

Comments

Post a Comment

Popular posts from this blog

cakephp - simple blog with croogo -

i not problem, made simple blog in croogo, problem is not adding or modifying or delete, because arabic language ?? function add in controller `public function admin_add() { $this->set('title_for_layout', __('add part')); if (!empty($this->request->data)) { $this->part->create(); if ($this->part->save($this->request->data)) { $this->session->setflash(__('the part has been saved'), 'default', array('class' => 'success')); $this->redirect(array('action' => 'index'));//, $this->part->id)); } else { $this->session->setflash(__('the part not saved. please, try again.'), 'default', array('class' => 'error')); } } }` => model `class part extends appmodel { public $name = 'part'; var $belongst...
Read more

How to group boxplot outliers in gnuplot -

i have large set of data points. try plot them boxplot, of outliers exact same value , represented on line beside each other. found how set horizontal distance between outliers in gnuplot boxplot , doesn't much, apparently not possible. is possible group outliers together, print 1 point , print number in brackets beside indicate how many points there are? think make more readable in graph. for information, have 3 boxplots 1 x value , times 6 in 1 graph. using gnuplot 5 , played around pointsize, doesn't reduce distance anymore. hope can help! edit: set terminal pdf set output 'dat.pdf' file0 = 'dat1.dat' file1 = 'dat2.dat' file2 = 'dat3.dat' set pointsize 0.2 set notitle set xlabel 'x' set ylabel 'y' header = system('head -1 '.file0); n = words(header) set xtics ('' 1) set [i=1:n] xtics add (word(header, i) i) set style data boxplot plot file0 using (1-0.25):1:(0.2) boxplot lw 2 lc rgb '#8b0000...
Read more

bash - Performing variable substitution in a string -

i have string contains use of variable, , i'd substitute variable's value string. right best have is #!/bin/bash foo='$name; echo bar' name="the name" expanded="$(eval echo "$foo")" echo "$expanded" which has obvious defects: prints the name bar while i'd print the name; echo bar instead of eval can bash's regex matching bash's string replacement: foo='$name; echo bar' name="the name" [[ "$foo" =~ \$([[:alnum:]]+) ]] && s="${!bash_rematch[1]}" && expanded="${foo/\$${bash_rematch[1]}/$s}" echo "$expanded" name; echo bar
Read more