similar to: quantile function for noncentral f-distribution

Dear R users, I am searching for a general minimisaton function like 'ms' in R. Please give me a hint! -- ------------------------------------------------------------ Dipl. Inform. J. Hedderich Institut f?r Medizinische Informatik und Statistik im Klinikum an der CAU Phone : 0431/5973182 Brunswiker Str. 10

Hello R friends, I am looking for a function to perform the Siegel-Tukey test in R. Can anyvbody help me? Thanks so much! J. Hedderich -- --------------------------------------------------------------- Dipl. Inform. J. Hedderich Institut f??r Medizinische Informatik und Statistik im Klinikum S.H. - Campus Kiel - Phone : 0431/5973182 Brunswiker Str.

Hi, I am new to lattice graphics and have the following question. I have a dataframe d2: cat PointsInTime MeanScore 1 Parent Intake 26.25000 2 Youth Intake 9.75000 3 Worker Intake 20.63636 4 Parent Discharge 24.00000 5 Youth Discharge 15.60000 6 Worker Discharge 23.90909 and to produce a plot I am using the following code print(xyplot(MeanScore ~

The trailing backslash seems to be the culprit. After looking at rw1020\bin\INSTALL, I noticed $pkg =~ s/\/$//; which I assume is to strip the trailing slash from a package name before obtaining the package name and cleaning out its directory I added $pkg =~ s/\\$//; immediately after this line. Rerunning Rcmd INSTALL C:\rw1011\src\library\bqtl\ with that change seemed to work OK

Ok, one last try. R1.2.0, WinNT 4.0: [R is installed in: R_HOME=D:\Programme\R\rw1020] Can anyone give a hint what the following error message (when running Rcmd INSTALL) means and what to do about it? cd D:\Programme\R\rw1020\src\library D:\Programme\R\rw1020\bin\rcmd install testfunctions make: Entering directory `/cygdrive/d/Programme/R/rw1020/src/gnuwin32' dlltool -k --as as --dllname

Dear R users, this is a question for R version 1.2.0 under Windows NT 4.0, regarding the documentation of ones own functions. I have ceated a private library for my functions, which works all very well, apart from that I am not able to create appropriate help files. For illustration, here is a simple example of where I'm stuck: > # first, create a simple function f1: > f1 <-

Dear List: The square of the noncentral t-statistic with noncentrality parameter \delta is a noncentral F with noncentrality parameter \lambda=\delta^2. So, t^2_{\nu,\delta} = F_{1,\nu,\lambda=\delta^2}. Consequently, it should follow that t^2_{1-\alpha/2,\nu,\delta} = f_{1-alpha,1,\vu,\lambda=\delta^2}. However, this is not what is happening with the following code. The central

Hello I could not find whether there is any R-function that implements the inverse of a noncentral Beta. Could someone out there tell me where I can find it? Or how to implement it? Many thanks Ed [[alternative HTML version deleted]]

Hi, Actually I am facing a similar problem. I would like to fit both an ordinary (symmetric) and a non-central t distribution to my (one-dimensional) data (quite some values.. > 1 mio.). For the symmetric one, fitdistr or funInfoFun (using fitdistr) from the qAnalyst package should do the job, and for the non-central one.. am I right to use gamlss(x ~ 1, family=GT()) ? Anyway, I am a little

Hi there: I'd thought these two versions of noncentral t-distribution are essentially the same: > qqplot(rt(1000,df=20,ncp=3),qt(runif(1000),df=20,ncp=3)) But, the scales of the x-axis and the y-axis are quite different according to the QQ-plot. Did I make any mistakes somewhere? Thanks, Long ---------------------------------

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,

For those who are interested, I have made available a R function for noncentral hypergeometric distribution at http://www.geocities.com/jg_liao/software/Hypergeometric/hypergeometric_in_R.txt The paper that describes the algorithm will appear in The American Statistician. The function does not run on S-plus as the R's scoping rule is used. Here is how the function can be used: > n1

Hi, I have written out the log-likelihood function to fit some data I have (called ONES20) to the non-central chi-squared distribution. >library(stats4) >ll<-function(lambda,k){x<-ONES20; 25573*0.5*lambda-25573*log(2)-sum(-x/2)-log((x/lambda)^(0.25*k-0.5))-log(besselI(sqrt(lambda*x),0.5*k-1,expon.scaled=FALSE))} > est<-mle(minuslog=ll,start=list(lambda=0.05,k=0.006))

I have a data with 236 observations. After plotting the histogram, I found that it looks like non-central t distribution. I would like to get MLE for mu and df. I found an example to find MLE for gamma distribution from "fitting distributions with R": library(stats4) ## loading package stats4 ll<-function(lambda,alfa) {n<-200 x<-x.gam

Hello all. It seems that the pf() function when used with noncentral parameters can behave badly at times. I've included some examples below, but what is happening is that with some combinations of df and ncp parameters, regardless of how large the quantile gets, the same probability value is returned. Upon first glance noncentral values greater than 200 may seem large, but they are in

Does anybody have any experience fitting data to the non-central chi-squared or chi-squared distribution? I am trying to fit some data to this distribution but there is error after error. audaces fortuna iuvat --------------------------------- [[alternative HTML version deleted]]

Hi, I would like to know if there is a program written in R to get the CDF (cumulative distribution function) of a bivariate non-central chi-square distribution. Hope someone will reply. Thank you, Rossita M Yunus yunus@usq.edu.au This email (including any attached files) is confidentia...{{dropped:19}}

Full_Name: Mike Meredith Version: 2.3.0 OS: WinXP SP2 Submission from: (NULL) (210.195.228.29) Using dt() with a non-centrality parameter and near-zero values for 'x' results in erratic output. Try this: tst <- c(1e-12, 1e-13, 1e-14, 1e-15, 1e-16, 1e-17, 0) dt(tst,16,1) I get: 0.2381019 0.2385462 0.2296557 0.1851817 0.6288373 3.8163916 (!!) 0.2382217 The 0.238 values are okay,

I have re-built all the compiled packages for R 1.2.0 for Windows, and the set available in CRAN/bin/windows/contrib are now for rw1020. If you are upgrading to rw1020, please replace all the compiled packages you have downloaded by ones from this set. There a few packages available for the first time: RMySQL mlbench netCDF and a few that are not yet running on 1.2.0 on any platform hdf5

I have re-built all the compiled packages for R 1.2.0 for Windows, and the set available in CRAN/bin/windows/contrib are now for rw1020. If you are upgrading to rw1020, please replace all the compiled packages you have downloaded by ones from this set. There a few packages available for the first time: RMySQL mlbench netCDF and a few that are not yet running on 1.2.0 on any platform hdf5