similar to: finding ncp for t distribution

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

Dear useRs, I have written a very simple function to compute some probabilities on words (function is below). The function includes a which() statement applied to a vector of characters (word.split): sapply (word.split, function(x) which(letters==x)). This statement worked as expected when used outside the global function : > word <- "hello" > (word.split <-

Dear useRs, I have written a function that implements a Bayesian method to compare a patient's score on two tasks with that of a small control group, as described in Crawford, J. and Garthwaite, P. (2007). Comparison of a single case to a control or normative sample in neuropsychology: Development of a bayesian approach. Cognitive Neuropsychology, 24(4):343?372. The function (see

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 expeRts, I would like to use a oneway_test (from package coin) to test whether two groups differ on various variables. The variables are encoded within a data frame. Unfortunately, I obtained an error, that I don't understand. Could you please help me ? Example: library(coin) y <- as.data.frame(matrix(rnorm(200), ncol=2)) group <- as.factor(unif(100)) lapply(y,

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

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

?pt states that ncp=0 is default, so this is correct > pt(3,3) [1] 0.9711656 > pt(3,3, ncp=0) [1] 0.9711656 but this is not > pt(-Inf,3) [1] 0 > pt(-Inf,3,ncp=0) [1] 0.5 is it? Torsten --please do not edit the information below-- Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu system = i686, linux-gnu status = Under development (unstable) major = 1 minor

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

Hi, I am running cwrsync. 2004/05/11 [93] rsync error: some files could not be transferred (code 23) at /home/lapo/packaging/tmp/rsync-2.5.7/main.c(383) I am uncertain as to what this path statement is about in the above error message. /home/lapo/packaging/tmp/rsync-2.5.7/main.c(383) Thanks, Pat -- Patricia Palumbo DuBois & King, Inc. ppalumbo@dubois-king.com 802-728-4113 | ext 322

Does anybody has a function or idea how to calculate qt with a noncentrality parameter? I'm porting a SAS-macro for equivalence (Wellek, 1994 TT2ST) to R and it seems that I need this to calculate tinv (SAS-command, I hope that this is R's qt,...I get the SAS-macro manual next week). thanks Peter -- P.Malewski, Limmerstr.47, 30451 Hannover, 0511-2135008 At work:

Hi everyone, I tried to use yum to search for ncpfs-2.2.1-1.i386.rpm, but it couldn't find it. I googled and found this package is included with redhat,but i'm not sure if it's included with CentOS 4.2. Any idea where I can find and install it using yum? Thanks JC

Is there any package that allow you to perform "MA plot" like graph for Toray microarray data? Unlike Affymetrix CEL file which contain 2 values (R and G), Torray raw data only contain 1 value. MA-plot is Affymetrix specific which usually available for in (limma) package. P. Dubois [[alternative HTML version deleted]]

Hello helpful R folks, I am simply trying to graph a quarter circle centered at the origin in the first quadrant. When I set the xlim of the plot to the radius of the circle, the plot appears correct. However, I'd like to see a slight extension of the axes beyond the domain of the function itself. When I do this, a portion of the plot seems to be missing by the edge of the domain. Here is

Hi all, I'm having trouble getting access to the pwr. This is on Ubuntu Lucid Lynx, 64 bit. I'm installing pwr via packages.install('pwr'), and loading it via library(pwr), both of which appear successful. Strangely, I never get access to the pwr object in R. I tried installing it to /usr/local/lib/R/site-library and ~/R/x86_64-pc-linux-gnu-library. I also tried installing it

Dear Christian, thanks for writing. I am no longer active in driver development, so I am forwarding your question to the NUT user list. The people on the list are very helpful. Best wishes, -- Peter Christian Dubois wrote: > > Hi there, I found this document on the internet and it was very useful > for me (http://www.mscs.dal.ca/~selinger/ups/backups.html), thank you. > > I

Hey there, I want to plot 5 parabola functions, which happen to be f(x) = 0.25x? + 6,47x -32.6 g(x)=0.99x? -6x -195 j(x)= 0.77x? +14x -495 k(x)=0.001x? + 65x -785 l(x) = 0.9x? -2x -636 in the same graph. Sadly I even do not really understand how to plot just one graph... I found this code in the Internet, which plots a cos-function: >x <- seq( -10, 10, length = 1000) >plot(x,

Hello all, I apologize in advance if this is not the correct list to post these sorts of issues. When attempting to install the latest stable build of R (2.6) or the most recent nightly build from R.research.att.com, the installer disables the ability to actually install R or any of the other bundled frameworks (i.e. the list is unchecked and grayed out). Before I start trying to either coerce the

Dear R Users, I am trying to plot a barchart with a line graph superimposed (using par(new=TRUE)). There are 12 bars and 12 corresponding points for the line graph. This is fine, except that I'm encountering two problems: 1) The position of the points (of the line graph) are not centred on the middle (horizontally) of each corresponding bar. In fact, whilst the first point is located on the

Your response nicely clarifies a question that I've had for a long time, but which I've dealt with by giving each subject a unique label. Unless I'm missing something, both techniques should work as the toy example below gives exactly the same output in all 3 cases below (forgetting about the convergence problem). Would there be a reason to prefer labeling the levels one way or