Displaying 20 results from an estimated 9000 matches similar to: "Non-central distributions"
2002 Aug 24
2
Density of non-central t distribution
Hello Everybody,
Is there a function available for calculating the density of a
non-central t distribution? (dt does not accept the ncp option).
Thanks in advance!
--
Wolfgang Viechtbauer
"Are you not thinking what I am not thinking?"
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read
2005 Apr 26
1
Finding 'ncp' for t
Hi Folks,
I'm looking for a neat procedure for the following:
Given t0 such that
pt(t0,df,ncp=0) = alpha (given)
find ncp0 such that for given beta
pt(t0,df,ncp=ncp0) = (1 - beta)
(In other words, what's the ncp such that you get power (1-beta)
to detect it, using a 1-sided test with size alpha when ncp = 0?)
In the past I've done the groping by hand, but this time it
needs
2006 Jul 01
1
noncentral F-distributed random numbers (PR#9055)
Full_Name: Long Qu
Version: 2.3.1
OS: Windows XP
Submission from: (NULL) (64.113.93.235)
The QQ-plot of two versions of simulating noncentral F-distributed random
numbers has quite different scales:
> qqplot(rf(1000,2,15,3),qf(runif(1000),2,15,3))
The rf() function reads:
> rf
function (n, df1, df2, ncp = 0)
{
if (ncp == 0)
.Internal(rf(n, df1, df2))
else rchisq(n, df1,
2001 May 18
1
Non-Central t
In the help file for the non-central t, the following appears:
ncp: non-centrality parameter delta; currently `ncp <= 37.62'.
I assume that this means the ncp cannot exceed 37.62. Is this
still the case and is there any plans to increase this restriction?
Thanks!
Jeff
Jeff Morris
Design Support
Clinical Chemistry R&D
Ortho-Clinical Diagnostics
email: jmorris6 at ocdus.jnj.com
2000 Nov 28
1
non-centrality parameter in pf() (PR#752)
Bug Description:
Problem with the function pf() when the non-centrality
parameter is large. Here is a sample command. You should
see a smooth line from 0 to about 55, and then the values
of pf() go crazy from 55 to 100.
############################
ncp <- seq(0,100,length=200)
plot(ncp,pf(5,7,2,ncp=ncp))
############################
Version:
platform = i686-pc-linux-gnu
arch = i686
os =
2002 Dec 17
4
Quick tip please!
I have two CSV files (exported from Excel), say file1 and file2.
The have the same number of rows, and each has several columns,
with names on the first line; and some of the columns in file1
are repeated in file2.
Using the "foreign" package, I can read these in separately
to dataframes say d1 and d2 with
> d1<-read.csv("file1")
>
2007 Jul 16
5
formula(CO2)
The formula attribute of the builtin CO2 dataset seems a bit strange:
> formula(CO2)
Plant ~ Type + Treatment + conc + uptake
What is one supposed to do with that? Certainly its not suitable for
input to lm and none of the examples in ?CO2 use the above.
2008 Jun 14
1
qt with ncp>37.62
help(qt) states that:
"ncp non-centrality parameter delta; currently except for rt(), only for
abs(ncp) <= 37.62"
so I would expect that calling qt with non-centrality parameter exceeding
37.62 should fail, instead e.g. calling
> mapply(function(x) qt(p = 0.9, df = 55, ncp = x),35:45)
gives:
[1] 40.21448 41.35293 42.49164 43.68862 44.82945 45.97048 47.11170 48.25310
[9]
2003 Jul 16
2
Density function for non-central t distribution
Hi,
I've written some C code for density evaluation of the non-central t distribution. It works
with the R-1.7.1 source code if placed in the src/nmath directory and after appropriate
changes are made to Makefiles, to the dt function in src/library/base/R/distn.R etc.
I haven't read a lot of R source code so it may need some R-ification, but I've tried to use the
dt.c file as a
2003 Apr 13
2
Peculiarity in non-central qchisq for ncp > 294.92 ...
Hello all,
Here's my query:
Running R 1.6.2 on FreeBSD 5.0, and on WinXP, and I find that the
following hangs the process:
dchisq(alpha=0.01, df=1, ncp=295)
but it does work for ncp < 294.92.
Is this general?
Best wishes to all,
Andrew
Andrew Robinson Ph: 208 885 7115
Department of Forest Resources Fa: 208 885 6226
University of Idaho E : andrewr at uidaho.edu
PO
2005 May 11
7
bug in modulus operator %% (PR#7852)
The following can't be right,
first rw2010:
> 1 %% 0.001
[1] 0.001
Then rw2001:
> 1 %% 0.001
[1] -2.081668e-17
>
and the last seems about right.
--
Kjetil Halvorsen.
Peace is the most effective weapon of mass construction.
-- Mahdi Elmandjra
--
No virus found in this outgoing message.
Checked by AVG Anti-Virus.
2002 Dec 10
3
clogit and general conditional logistic regression
Can someone clarify what I cannot make out from the
documentation?
The function 'clogit' in the 'survival' package is
described as performing a "conditional logistic regression".
Its return value is stated to be "an object of class clogit
which is a wrapper for a coxph object."
This suggests that its usefulness is confined to the sort of
data which arise in
2004 Dec 26
2
Prosodic/phonetic analysis with R
Hi Folks,
I'm interested in looking (in a beginner's amateurish way)
at prosodic/phonetic analysis of recorded speech.
In particular I would want to use R to formulate and
evaluate specific models.
So I would like to ask R people for their recommendations
for a program which would
a) Take as input a sound file in one of the common formats
(".wav", ".au")
b)
2007 Jun 15
3
[OT] 'gv' and fractional points
Hi Folks,
This is off-topic R-wise, but it may be close to
the heart of many R-users, so I think it may be
the best place to ask!
Users of 'gv' (the "front end" to ghostscript) will
be aware of the little window which gives you the
x-y coordinates (in points = 1/72 inch) of the position
of the "cross-hair" mouse cursor. These coordinates
are those of the corresponding
2003 Jun 25
2
Execution of R code
Greetings Folks,
When R code (as entered or read from a courced file) is executed,
is it interpreted from the input form every time having once been
read in, or do subsequent invocations use an "intermediate"
(pre-interpreted) form?
Or, putting it another way, is the execution of R code faster
second time time round (and later) because the pre-interpretation
has already been done once
2005 Apr 21
9
Using R to illustrate the Central Limit Theorem
Dear All
I am totally new to R and I would like to know whether R is able and
appropriate to illustrate to my students the Central Limit Theorem,
using for instance 100 independent variables with uniform distribution
and showing that their sum is a variable with an approximated normal
distribution.
Thanks in advance,
Paul
2007 May 01
2
Simulation using parts of density function
Hi
My simulation with the followin R code works perfectly:
sim <- replicate(999, sum(exp(rgamma(rpois(1,2000), scale = 0.5, shape = 12))))
But now I do not want to have values in object "sim" exceeding 5'000'000, that means that I am just using the beginning of densitiy function gamma x < 15.4. Is there a possibility to modify my code in an easy way?
Thanks for any help!
2003 Oct 09
4
Automatic re-looping after error
Hi Folks,
I'm seeking advice about how to resume an outer loop following
failure of a function within the loop (which issues an error
message).
Essentially, I'm repeating (simulating) a process which involves
random sampling, EM and MCMC. I'm walking on very edge of rather
thin ice -- data rather thinly spread over many parameters!
Occasionally some component of the loop fails, with
2004 May 04
2
Superposing data on boxplot
Hi folks,
I have a vaiable Y and an associated factor Z at several (13)
levels.
boxplot(Y~Z)
produces a nice array of boxplots, one for each level of Z,
and each duly labaelled with its level of Z.
I would like to superpose on each boxplot the actual data
points which it represents, i.e. do something conceptually
(though not in real R) expressed as
points(Y~Z)
or
points(Z,Y)
It can
2005 Feb 26
2
R-help Time Series
Hi Folks,
While I was browsing in the R-help archives yesterday,
I got curious about the time series of the sizes of
the monthly archives in MB.
This turned out to have an unexpected feature or two,
which I leave to readers to explore for themselves.
I'm now wondering at what point in time we might expect
to be receiving 1000MB/month (30+MB/day). It's not that
far away, it seems, but