similar to: How to shadow 'power' area?

Dear R useRs, I have a vector with positive and negative numbers: A=c(0,1,2,3,0,4,5) Now if i-th element in vector A is > 0, then i-th element in vector B is a+1 else i-th element in vector b=a (or 0) vector A: 0 1 2 3 0 4 5 vector B: 0 2 3 4 0 5 6 What's the right way to do this. I still have some problems with for and if statements... Cheers, Andrej

Dear useRs, I have two different length vectors: one column (1...m) and one row vector (1...n): 20 40 20 60 5 4 2 Now I have to calculate ratios between column vector elements and each row vector elements: 4 5 10 8 10 20 4 5 20 15 12 30 Thank's in advance for any suggestions, Andrej

Hello all, I have a question about basic statistics. Given a PDF value of 0.328161, how can I find out the value of -0.625 in R? It is like reversing the dnorm function but I do not know how to do it in R. > pdf.xb <- dnorm(-0.625) > pdf.xb [1] 0.328161 > qnorm(pdf.xb) [1] -0.444997 > pnorm(pdf.xb) [1] 0.628605 Many thanks, Edwin -- View this message in context:

Dear R useRs, I have big (23000 rows), vertical bar delimited file: e.g. A00001|Text a,Text b, Text c|345 A00002|Text bla|456 ... .. . Try using A <- read.table('filename.txt', header=FALSE,sep='\|') process stop at line 11975 with warning message: number of items read is not a multiple of the number of columns I have no problems with processing similar file, which is

Hi, Having a weird problem with the integrate function. I have a function which calculates a loss density: I'd like to integrate it to get the distribution. The loss density function is: lossdensity<-function(p,Beta,R=0.4){ # the second derivative of the PDF # p is the default probability of the pool at which we are evaluating the lossdensity # Beta is the correlation with the market

Hi, Am not sure if my code itself is correct. Here's what am trying to do: Minimize integration of a function of gaussian distributed variable 'x' over the interval qnorm(0.999) to Inf by changing value of parameter 'mu'. mu is the shift in mean of 'x'. Code: # x follows gaussian distribution # fx2 to be minimized by changing values of mu # integration to be done over

Hi, This problems has bothered me for the lase couple of hours. > 1e-100==0 [1] FALSE > (1-1e-100)==1 [1] TRUE How can I tell R that 1-1e-100 does not equal to 1, actually, I found out that > (1-1e-16)==1 [1] FALSE > (1-1e-17)==1 [1] TRUE The reason I care about this is that I was try to use qnorm() in my code, for example, > qnorm(1e-100) [1] -21.27345 and if I want to

Dear all, Does anybody have an idea or suggestion how to construct (plot) 4-dimensional hypercube in R. Thanks in advance for any pointers. Regards, Andrej

Dear R useRs, is there any simple, build in function to match specific regular expression in data file and write it to a vector. I have the following text file: *NEW RECORD *ID-001 *AB-text *NEW RECORD *ID-002 *AB-text etc. Now I have to match all ID fields and print them to a vector: 001 002 etc. I know that this is very simple with Perl or R-Perl interface, but if possible, I want to do

Hello, could You please help: I am looking for a way to formulate test accuracy measures such as test sensitivity, specificity, predictive values, and correct classification rate using p/d/qnorm. The tests' primary values follow a bimodal distribution, which is modelled by a mixture of two normal distributions: p * dnorm ((x - u1) / s1) / s1 + (1 - p) * dnorm ((x - u2) / s2) / s2)

Is there a function to compute the derivative of the probit (qnorm) function in R, or in any of the packages? Thanks, -Andrew [[alternative HTML version deleted]]

Dear All I am trying to find a uniroot of a function within another function (see example) but I am getting an error message (f()values at end points not of opposite sign). I was wondering if you would be able to advise how redefine my function so that I can find the solution. In short my first function calculates the intergrale which is function of "t" , I need to find the uniroot of

dear all, here's a couple of questions that puzzled me in these last hours: ##### issue 1 qnorm(1-10e-100)!=qnorm(10e-100) qnorm(1-1e-10) == -qnorm(1e-10) # turns on to be FALSE. Ok I'm not a computer scientist but, # but I had a look at the R inferno so I write: all.equal(qnorm(1-1e-10) , -qnorm(1e-10)) # which turns TRUE, as one would expect, but all.equal(qnorm(1-1e-100) ,

Dear all, I try to install the latest R version using checkinstall (v. 1.6.2) on Ubuntu 11.10. After solving all the dependencies (installed using apt-get build-dep r-base) checkinstall fails to build and install R package with the following error (the same commands build and install R-2.13.2 on the same machine without any problem): wget

Dear list, has anyone written a package/function in R for computing a point- biserial resp. biserial correlation? Thanks in advance Bernd

Here's another example of my plotmath whipping boy, the Normal distribution. A colleague asks for a Normal plotted above a series of axes that represent various other distributions (T, etc). I want to use vectors of equations in plotmath to do this, but have run into trouble. Now I've isolated the problem down to a relatively small piece of working example code (below). If you would

Full_Name: Xiong Guanglei Version: 1.8.1 OS: Linux Submission from: (NULL) (202.38.103.50) qnorm(x) when x>1.0

Hi, In R, qnorm is the inversion function of pnorm. (c.d.f) But there is no inversion function for dnorm. (p.d.f). Is there any easy (and quick) way to compute the inversion function of p.d.f in R ? Thanks ...

Dear all, I'm working on ability estimates using Rasch model. Using the "ltm" package, the procedure is quite simple: ## Factor Scores for the Rasch model fit <- rasch(LSAT) factor.scores(fit) What about Partial Credit Model (PCM)? For PCM I use PCM function from eRm package. Is there any similar function like factor.scores to estimate ability scores using PCM model? Best,

(1) ...need to speed up a monte-carlo sampling...any suggestions about how I can get R to use all 8 cores of a mac pro would be most useful and very appreciated... (2) spent the last few hours trying to get pnmath to compile under os- x 10.5.4... using gcc version 4.2.1 (Apple Inc. build 5553) as downloaded from CRAN, xcode 3.0... ...xcode 3.1 installed over top of above after