Integrate(dnorm) with different mean and standard deviation help

11 messages
Open this post in threaded view
|

Integrate(dnorm) with different mean and standard deviation help

 I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow.   > integrate(dnorm, -1.96, 1.96) 0.9500042 with absolute error < 1e-11 > mean = 2.0 > sd = 1.0 > integrate(dnorm, -1.96, 1.96) 0.9500042 with absolute error < 1e-11 > integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing Calls: integrate -> match.fun -> dnorm Execution halted   How do I change the built in mean=0 and standard deviation=1 for dnorm using integrate?   Thanks,   Frank Chicago, IL             [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 Hi, Check this link (https://stat.ethz.ch/pipermail/r-help/2010-February/227902.html). Hope it helps. A.K. ----- Original Message ----- From: FJ M <[hidden email]> To: R <[hidden email]> Cc: Sent: Monday, July 23, 2012 10:23 PM Subject: [R] Integrate(dnorm) with different mean and standard deviation help I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow. > integrate(dnorm, -1.96, 1.96) 0.9500042 with absolute error < 1e-11 > mean = 2.0 > sd = 1.0 > integrate(dnorm, -1.96, 1.96) 0.9500042 with absolute error < 1e-11 > integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing Calls: integrate -> match.fun -> dnorm Execution halted How do I change the built in mean=0 and standard deviation=1 for dnorm using integrate? Thanks, Frank Chicago, IL                               [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code. ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 In reply to this post by FJ M Hello, Maybe the following could help:  > f <- function(x) dnorm(x, mean=2, sd=1)  > integrate(f, -1.96, 1.96) 0.4840091 with absolute error < 1.4e-12 HTH Regards. Le 24/07/2012 11:23, FJ M a écrit : > > I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow. > > >> integrate(dnorm, -1.96, 1.96) > 0.9500042 with absolute error < 1e-11 >> mean = 2.0 >> sd = 1.0 >> integrate(dnorm, -1.96, 1.96) > 0.9500042 with absolute error < 1e-11 >> integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) > Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing > Calls: integrate -> match.fun -> dnorm > Execution halted > > How do I change the built in mean=0 and standard deviation=1 for dnorm using integrate? > > Thanks, > > Frank > Chicago, IL   > [[alternative HTML version deleted]] > > ______________________________________________ > [hidden email] 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. > ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 In reply to this post by FJ M Try > integrate(dnorm, mean = 2, sd = 1, -1.96, 1.96) 0.4840091 with absolute error < 1.4e-12 HTH, Jorge.- On Mon, Jul 23, 2012 at 10:23 PM, FJ M <> wrote: > > I'm trying to provide different parameters to the integrate function for > various probability functions. I'm using dnorm as the simplest example > here. For instance integrate(dnorm, -1.96, 1.96) produces the correct > answer for a normal distribution with mean 0 and standard deviation 1. I've > tried two ways to use mean=2.0 and standard deviation 1, but with no luck. > The examples follow. > > > > integrate(dnorm, -1.96, 1.96) > 0.9500042 with absolute error < 1e-11 > > mean = 2.0 > > sd = 1.0 > > integrate(dnorm, -1.96, 1.96) > 0.9500042 with absolute error < 1e-11 > > integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) > Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing > Calls: integrate -> match.fun -> dnorm > Execution halted > > How do I change the built in mean=0 and standard deviation=1 for dnorm > using integrate? > > Thanks, > > Frank > Chicago, IL >         [[alternative HTML version deleted]] > > ______________________________________________ > [hidden email] 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. >         [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 In reply to this post by Pascal Oettli-2 On 2012-07-23 19:48, Pascal Oettli wrote: > Hello, > > Maybe the following could help: > >   > f <- function(x) dnorm(x, mean=2, sd=1) >   > integrate(f, -1.96, 1.96) > 0.4840091 with absolute error < 1.4e-12 Or you could note the '...' argument indicated on the help page:    integrate(dnorm, lower = -1.96, upper = 1.96,              mean = 2, sd = 1) Peter Ehlers > > HTH > > Regards. > > > Le 24/07/2012 11:23, FJ M a écrit : >> >> I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow. >> >> >>> integrate(dnorm, -1.96, 1.96) >> 0.9500042 with absolute error < 1e-11 >>> mean = 2.0 >>> sd = 1.0 >>> integrate(dnorm, -1.96, 1.96) >> 0.9500042 with absolute error < 1e-11 >>> integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) >> Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing >> Calls: integrate -> match.fun -> dnorm >> Execution halted >> >> How do I change the built in mean=0 and standard deviation=1 for dnorm using integrate? >> >> Thanks, >> >> Frank >> Chicago, IL   >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> [hidden email] 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. >> > > ______________________________________________ > [hidden email] 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. > ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 In reply to this post by FJ M integrate(dnorm, -1.96, 1.96, mean=2, sd=1) Read the help for integrate!  It tells you that the integrate function has a "..." argument which consists of "additional arguments to be passed to f".      cheers,          Rolf Turner On 24/07/12 14:23, FJ M wrote: > I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow. >   > >> integrate(dnorm, -1.96, 1.96) > 0.9500042 with absolute error < 1e-11 >> mean = 2.0 >> sd = 1.0 >> integrate(dnorm, -1.96, 1.96) > 0.9500042 with absolute error < 1e-11 >> integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) > Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing > Calls: integrate -> match.fun -> dnorm > Execution halted >   > How do I change the built in mean=0 and standard deviation=1 for dnorm using integrate? ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 In reply to this post by Peter Ehlers Hello, Yes, I should learn to read. Regards Le 24/07/2012 11:54, Peter Ehlers a écrit : > On 2012-07-23 19:48, Pascal Oettli wrote: >> Hello, >> >> Maybe the following could help: >> >>   > f <- function(x) dnorm(x, mean=2, sd=1) >>   > integrate(f, -1.96, 1.96) >> 0.4840091 with absolute error < 1.4e-12 > > Or you could note the '...' argument indicated on the help page: > >    integrate(dnorm, lower = -1.96, upper = 1.96, >              mean = 2, sd = 1) > > Peter Ehlers > >> >> HTH >> >> Regards. >> >> >> Le 24/07/2012 11:23, FJ M a écrit : >>> >>> I'm trying to provide different parameters to the integrate function >>> for various probability functions. I'm using dnorm as the simplest >>> example here. For instance integrate(dnorm, -1.96, 1.96) produces the >>> correct answer for a normal distribution with mean 0 and standard >>> deviation 1. I've tried two ways to use mean=2.0 and standard >>> deviation 1, but with no luck. The examples follow. >>> >>> >>>> integrate(dnorm, -1.96, 1.96) >>> 0.9500042 with absolute error < 1e-11 >>>> mean = 2.0 >>>> sd = 1.0 >>>> integrate(dnorm, -1.96, 1.96) >>> 0.9500042 with absolute error < 1e-11 >>>> integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) >>> Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing >>> Calls: integrate -> match.fun -> dnorm >>> Execution halted >>> >>> How do I change the built in mean=0 and standard deviation=1 for >>> dnorm using integrate? >>> >>> Thanks, >>> >>> Frank >>> Chicago, IL >>>     [[alternative HTML version deleted]] >>> >>> ______________________________________________ >>> [hidden email] 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. >>> >> >> ______________________________________________ >> [hidden email] 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. >> > > ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|

Re: Integrate(dnorm) with different mean and standard deviation help

 In reply to this post by Peter Ehlers It would be a useful additon to the help page to add   integrate(dnorm, lower = -1.96, upper = 1.96, mean = 2, sd = 1) as an example.   Thanks,   Frank Chicago   > Date: Mon, 23 Jul 2012 19:54:45 -0700 > From: [hidden email] > To: [hidden email] > CC: [hidden email]; [hidden email] > Subject: Re: [R] Integrate(dnorm) with different mean and standard deviation help > > On 2012-07-23 19:48, Pascal Oettli wrote: > > Hello, > > > > Maybe the following could help: > > > > > f <- function(x) dnorm(x, mean=2, sd=1) > > > integrate(f, -1.96, 1.96) > > 0.4840091 with absolute error < 1.4e-12 > > Or you could note the '...' argument indicated on the help page: > > integrate(dnorm, lower = -1.96, upper = 1.96, > mean = 2, sd = 1) > > Peter Ehlers > > > > > HTH > > > > Regards. > > > > > > Le 24/07/2012 11:23, FJ M a écrit : > >> > >> I'm trying to provide different parameters to the integrate function for various probability functions. I'm using dnorm as the simplest example here. For instance integrate(dnorm, -1.96, 1.96) produces the correct answer for a normal distribution with mean 0 and standard deviation 1. I've tried two ways to use mean=2.0 and standard deviation 1, but with no luck. The examples follow. > >> > >> > >>> integrate(dnorm, -1.96, 1.96) > >> 0.9500042 with absolute error < 1e-11 > >>> mean = 2.0 > >>> sd = 1.0 > >>> integrate(dnorm, -1.96, 1.96) > >> 0.9500042 with absolute error < 1e-11 > >>> integrate(dnorm(mean=2.0,sd=1.0), -1.96, 1.96) > >> Error in .Internal(dnorm(x, mean, sd, log)) : 'x' is missing > >> Calls: integrate -> match.fun -> dnorm > >> Execution halted > >> > >> How do I change the built in mean=0 and standard deviation=1 for dnorm using integrate? > >> > >> Thanks, > >> > >> Frank > >> Chicago, IL > >> [[alternative HTML version deleted]] > >> > >> ______________________________________________ > >> [hidden email] 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. > >> > > > > ______________________________________________ > > [hidden email] 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. > > >               [[alternative HTML version deleted]] ______________________________________________ [hidden email] mailing list https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.
Open this post in threaded view
|