similar to: Non-Central t

Do you know what I really, really miss in R that I took fopr granted in my XLisp-Stat programming days??? The ability to easily create menus and dialog boxes. I could write little applications and give them out to clients. It was wonderful. It would at least allow each person to begin to build one's own menu driven version. Just a thought!! Jeff > -----Original Message----- >

I am trying to run, in R -Windows95, the example in MASS2, section 17.2. When I run predplot() on cush.lda I get > predplot(cush.lda, "LDA") Error: couldn't find function "max.col" I get the plot OK, but within predict(), it balks on max.col A more general question about using VR. Is there an easy way to use library on VR in general or do I have to move each into

Hello All, Has anyone written a function for the distribution function of a *doubly* non-central F-distribution? I looked through the archives, but didn't find anything. Thanks! Wolfgang

Greetings all, I'm trying to figure out how to calculate the inverse CDF (i.e. a quantile) for a non-central F distribution. I could put together a quick numerical solver routine using the CDF, but I wonder if there's a function that I've missed that would be more efficient? Thank-you, Andrew Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885

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]

Dear all, I am attempting to model some one-dimensional data using a mixture model of non-central Student's t distributions. However, I haven't been able to find any R package that provides this functionality. Could there be a way to "manipulate" the EM algorithms from the mixdist or mixtools package to fit the model, or do you have any other suggestions? If anyone could help

> Date: Mon, 11 Jan 1999 15:23:10 +0100 > From: Kjetil Halvorsen <khal at alumni.uv.es> > To: "Morris, Jeffrey (JCD-RL)" <JMorris at OCDUS.JNJ.COM> > CC: "'R-Help'" <r-help at stat.math.ethz.ch> > Subject: [R] R: Regression with t errors? > > Has anybody implemented/knows about ways to do in R regression > based on

I'm running Guido's latest version on NT 4.0 and when I try to copy a plot to the clipboard, I get an information window that says "Impossible to Open the Clipboard". Has anyone else seen this? BTW, it works fine on my laptop running Windows95. Thanks! Jeff jmorris at ocdus.jnj.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

I have downloaded an installed the xgobi library and can plot(laser), but when I try xgobi(laser), I get these messages in the MSDOS window... Out of environment space Out of environment space Bad command or filename and in the RGUI window it comes back with... D:\RW0901/library/xgobi/scripts/xgobi.bat -vtitle 'laser' C:/WINDOWS/TEMP/xgobi24464 Any ideas what I need to do? I am

Is there a easy way to put confidence bands on a regression plot or do I have to explictly construct the values and send the lines to the plot? Thanks! Jeff jmorris at ocdus.jnj.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

I am submitting the command: calcium<-read.table("calcium.txt",header=TRUE) and with Rw0630 and Bdr0631 I get an application error and R crashes. But the command works fine with rjune. I am running under NT4.0 Anyone have a similar problem? Jeff jmorris at ocdus.jnj.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi, I think I have identified a inaccuracy in dchisq when the non-centrality parameter is non-zero and large. Here's a little test: sys.dchisq.test <- function(N = 100000,mean = 0) { z <- rnorm(N,mean = mean, sd = 1) x <- z^2 xmin <- min(x) xmax <- max(x) br <- seq(xmin,xmax,length = 101) dbr <- br[2]-br[1] hist(x,br) p <- dchisq(br,df = 1,ncp =

I have a dataset that I am running non-linear regression on via the following code: Hill <- function(E0,Em,C50,g,C){ # # Hill is the hill interaction function. # # E0 Represents the minimum interaction Effect # # Em Represents the Maximum Interaction Effect # # C50 represents the concentration at which 50% of the effect occurs. # # gamma represents the cooperativity of the

What in the world am I missing?? > x<-rnorm(20) > mean(x) [1] -0.2272851 > results<-boot(x,mean,R=5) > results[2] $t [,1] [1,] -0.2294562 [2,] -0.2294562 [3,] -0.2294562 [4,] -0.2294562 [5,] -0.2294562 Jeff Morris Ortho-Clinical Diagnostics A Johnson & Johnson Co. Rochester, NY Tel: (585) 453-5794 [[alternative HTML version deleted]]

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

Dear Rusers, I tried to implement a function comparing mean scores between one subject (the patient) and a group a control subjects. The function returns attended results, but I also obtained the following warning : Warning message: the condition has length > 1 and only the first element will be used in: if (substr(fp, 1, 1) == "<") fp else paste("=", fp) Maybe

Hello: I am brand new to R, though I've used S-Plus for some time. I've developed some specialized plotting functions that I want to distribute and would like to provide users with a dialogbox for entering input arguments. Has anyone developed functions in R for displaying dialogboxes that could accept multiple inputs? I found winDialog and winDialogString, but these only accept and

Dear Professor I have read your questions in the website on Doubly non-central F-distribution. I am looking for source code for evaluating this function now. I have matlab code but it's only accurate up to the 4th decimal point. Dataplot is more accurate, but it is not user friendly as you only can evaluate one function at a time. So I would like to know do you have found any R code for

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 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]]