thr3ads.net - R help - [R] double exponential regression R [Apr 2013]

If this information is useful, please help other people find it:
Share via:

catalin roibu

2013-Apr-21 19:34 UTC

[R] double exponential regression R

Hello all!
I have a problem with a double exponential equation.
This are my data's> structure(list(proc = c(1870.52067384719,
766.789388745793, 358.701545859122,  237.113777545511, 43.2726259059654,
148.985133316262, 92.6242882655781,  88.4521557193262, 56.6404686159112,
27.0374477259404, 34.3347291080268,  18.3226992991316, 15.2196612445747,
5.31600719692165, 16.7015717397302,  16.3923389973684, 24.2702542054496,
21.247247993673, 18.3070717608672,  2.8811892177331, 3.18018869564679,
8.74204132937479, 7.11596966047229 ), cls = c(0.25, 0.5, 0.75, 1, 1, 1.5,
2, 2.5, 3, 3.5, 4, 4.5,  5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10)),
.Names = c("proc",  "cls"), row.names = c("0.25",
"0.5", "0.75", "1", "11",
"1.5",  "2", "2.5", "3",
"3.5", "4", "4.5", "5", "5.5",
"6", "6.5", "7",
 "7.5", "8", "8.5", "9",
"9.5", "10"), class = "data.frame")
I want to compute a double exponential equation like this:
proc=a*exp(b*class)+c*exp(d*class)

Thank you very much for your help!

Best regards!

CR
-- 
---
Catalin-Constantin ROIBU
Lecturer PhD, Forestry engineer
Forestry Faculty of Suceava
Str. Universitatii no. 13, Suceava, 720229, Romania
office phone     +4 0230 52 29 78, ext. 531
mobile phone   +4 0745 53 18 01
                       +4 0766 71 76 58
FAX:                +4 0230 52 16 64
silvic.usv.ro

	[[alternative HTML version deleted]]

Joshua Wiley

2013-Apr-21 20:09 UTC

head link

[R] double exponential regression R

Hi CR,

Assuming your data are stored in "d", perhaps something like this:

m <- nls(log(proc) ~ k + a*(b^cls), start = list(a = 2, b = .25, k 0), data =
d)
summary(m)
# gives
## Formula: log(proc) ~ k + a * (b^cls)
## Parameters:
##   Estimate Std. Error t value Pr(>|t|)
## a  5.24617    0.50905  10.306 1.90e-09 ***
## b  0.69640    0.07718   9.023 1.73e-08 ***
## k  1.93055    0.42331   4.561  0.00019 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05
'.' 0.1 ' ' 1
## Residual standard error: 0.6967 on 20 degrees of freedom
## Number of iterations to convergence: 8
## Achieved convergence tolerance: 4.268e-06

and you can examine the fit:

plot(fitted(m), log(d$proc))

which looks not too bad (see attached)

Cheers,

Joshua


On Sun, Apr 21, 2013 at 12:34 PM, catalin roibu <catalinroibu at
gmail.com> wrote:> Hello all!
> I have a problem with a double exponential equation.
> This are my data's> structure(list(proc = c(1870.52067384719,
> 766.789388745793, 358.701545859122,  237.113777545511, 43.2726259059654,
> 148.985133316262, 92.6242882655781,  88.4521557193262, 56.6404686159112,
> 27.0374477259404, 34.3347291080268,  18.3226992991316, 15.2196612445747,
> 5.31600719692165, 16.7015717397302,  16.3923389973684, 24.2702542054496,
> 21.247247993673, 18.3070717608672,  2.8811892177331, 3.18018869564679,
> 8.74204132937479, 7.11596966047229 ), cls = c(0.25, 0.5, 0.75, 1, 1, 1.5,
> 2, 2.5, 3, 3.5, 4, 4.5,  5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10)),
> .Names = c("proc",  "cls"), row.names =
c("0.25", "0.5", "0.75", "1",
"11",
> "1.5",  "2", "2.5", "3",
"3.5", "4", "4.5", "5", "5.5",
"6", "6.5", "7",
>  "7.5", "8", "8.5", "9",
"9.5", "10"), class = "data.frame")
> I want to compute a double exponential equation like this:
> proc=a*exp(b*class)+c*exp(d*class)
>
> Thank you very much for your help!
>
> Best regards!
>
> CR
> --
> ---
> Catalin-Constantin ROIBU
> Lecturer PhD, Forestry engineer
> Forestry Faculty of Suceava
> Str. Universitatii no. 13, Suceava, 720229, Romania
> office phone     +4 0230 52 29 78, ext. 531
> mobile phone   +4 0745 53 18 01
>                        +4 0766 71 76 58
> FAX:                +4 0230 52 16 64
> silvic.usv.ro
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.


-- 
Joshua Wiley
Ph.D. Student, Health Psychology
University of California, Los Angeles
http://joshuawiley.com/
Senior Analyst - Elkhart Group Ltd.
http://elkhartgroup.com
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fit.png
Type: image/png
Size: 5534 bytes
Desc: not available
URL:
<https://stat.ethz.ch/pipermail/r-help/attachments/20130421/62979957/attachment.png>

Jose Iparraguirre

2013-Apr-21 20:20 UTC

head link

[R] double exponential regression R

Dear Catalin,

I think the package dcemriS4 has this capability.
Hope it helps,
Jose
________________________________________
From: r-help-bounces at r-project.org [r-help-bounces at r-project.org] On
Behalf Of catalin roibu [catalinroibu at gmail.com]
Sent: 21 April 2013 20:34
To: r-help at r-project.org
Subject: [R] double exponential regression R

Hello all!
I have a problem with a double exponential equation.
This are my data's> structure(list(proc = c(1870.52067384719,
766.789388745793, 358.701545859122,  237.113777545511, 43.2726259059654,
148.985133316262, 92.6242882655781,  88.4521557193262, 56.6404686159112,
27.0374477259404, 34.3347291080268,  18.3226992991316, 15.2196612445747,
5.31600719692165, 16.7015717397302,  16.3923389973684, 24.2702542054496,
21.247247993673, 18.3070717608672,  2.8811892177331, 3.18018869564679,
8.74204132937479, 7.11596966047229 ), cls = c(0.25, 0.5, 0.75, 1, 1, 1.5,
2, 2.5, 3, 3.5, 4, 4.5,  5, 5.5, 6, 6.5, 7, 7.5, 8, 8.5, 9, 9.5, 10)),
.Names = c("proc",  "cls"), row.names = c("0.25",
"0.5", "0.75", "1", "11",
"1.5",  "2", "2.5", "3",
"3.5", "4", "4.5", "5", "5.5",
"6", "6.5", "7",
 "7.5", "8", "8.5", "9",
"9.5", "10"), class = "data.frame")
I want to compute a double exponential equation like this:
proc=a*exp(b*class)+c*exp(d*class)

Thank you very much for your help!

Best regards!

CR
--
---
Catalin-Constantin ROIBU
Lecturer PhD, Forestry engineer
Forestry Faculty of Suceava
Str. Universitatii no. 13, Suceava, 720229, Romania
office phone     +4 0230 52 29 78, ext. 531
mobile phone   +4 0745 53 18 01
                       +4 0766 71 76 58
FAX:                +4 0230 52 16 64
silvic.usv.ro

        [[alternative HTML version deleted]]

______________________________________________
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.

Please donate to the Syria Crisis Appeal by text or online:

To donate ?5 by mobile, text SYRIA to 70800.  To donate online, please visit 

http://www.ageinternational.org.uk/syria

Over one million refugees are desperately in need of water, food, healthcare,
warm clothing,
blankets and shelter; Age International urgently needs your support to help
affected older refugees.

Age International is a subsidiary charity of Age UK and a member of the
Disasters Emergency Committee (DEC).
The DEC launches and co-ordinates national fundraising appeals for public
donations on behalf of its member agencies.

Texts cost ?5 plus one standard rate message.  Age International will receive a
minimum of ?4.96.
More info at ageinternational.org.uk/SyriaTerms

-------------------------------
Age UK is a registered charity and company limited by guarantee, (registered
charity number 1128267, registered company number 6825798).
Registered office: Tavis House, 1-6 Tavistock Square, London WC1H 9NA.

For the purposes of promoting Age UK Insurance, Age UK is an Appointed
Representative of Age UK Enterprises Limited, Age UK is an Introducer
Appointed Representative of JLT Benefit Solutions Limited and Simplyhealth
Access for the purposes of introducing potential annuity and health
cash plans customers respectively.  Age UK Enterprises Limited, JLT Benefit
Solutions Limited and Simplyhealth Access are all authorised and
regulated by the Financial Services Authority. 
------------------------------

This email and any files transmitted with it are confidential and intended
solely for the use of the individual or entity to whom they are
addressed. If you receive a message in error, please advise the sender and
delete immediately.

Except where this email is sent in the usual course of our business, any
opinions expressed in this email are those of the author and do not
necessarily reflect the opinions of Age UK or its subsidiaries and associated
companies. Age UK monitors all e-mail transmissions passing
through its network and may block or modify mails which are deemed to be
unsuitable.

Age Concern England (charity number 261794) and Help the Aged (charity number
272786) and their trading and other associated companies merged
on 1st April 2009.  Together they have formed the Age UK Group, dedicated to
improving the lives of people in later life.  The three national
Age Concerns in Scotland, Northern Ireland and Wales have also merged with Help
the Aged in these nations to form three registered charities:
Age Scotland, Age NI, Age Cymru.

Possibly Parallel Threads

Search for more possibly parallel threads

R help - Apr 2013 - double exponential regression R

[R] double exponential regression R

[R] double exponential regression R

[R] double exponential regression R

Possibly Parallel Threads