similar to: Help: find the index of the minimum of entries

Dear All, I have a text file, tab delimited, called "sample.txt",as follows: ID_REF 382 GC_Score Theta R B_Allele_Freq Log_R_Ratio 200003 BB 0.9101527 0.9734979 0.8788951 1 0 200006 AB 0.6003323 0.4385073 2.033364 0.4850979 0.01553433 I have explored various options of the command: read.table, with one as:

Dear All, As always, I appreciate all your help. I would like to know the easiest way to convert each of the homogeneous elements of a numeric list into a matrix. Each element of this list is also a list such that when displayed, looks like a 2-by-3 matrix , I would like to convert each of them into a matrix, without changing the double index of each entry. Suppose: a<-

Dear All, I would like to know how to sort a vector of numeric values such that we know the original index of each ordered component. Say, we have c <- c(1,4,3,2) csort <- sort(c,descreasing=FALSE) With a few components of c, we can manually find out: csort[1] = 1 = c[1], ie, the original index of csort[1] is 1, csort[2] =2 =c[4], ie, the original index of csort[2] is 4. When length(c) is

Dear All, I would like to know, beside writing a function and then apply it to a list, or using a for loop, whether there is a one-line command to do the following. Suppose we have a list, each of whose components are numeric: > a <- vector("list",3) > a[[1]] <-c(1,2) > a[[2]] <-c(3,4) > a[[3]] <- c(5,6) > a [[1]] [1] 1 2 [[2]] [1] 3 4 [[3]] [1] 5 6 Target:

Dear All, I would like to ask a question about how to locate the indices after comparing two numeric vectors. Say, I have, A <- c(1,2,3,4) and A <- c(0,3,1,5) and we cmpare: > idx <- A<B So: > idx [1] FALSE TRUE FALSE TRUE Question 1: I would like to get a vector that stores the indices for entries of A such that this relationship is true. In this case, a vector idx1

Dear All, I would like to define a function: f(x,y,z) with three arguments x,y,z, such that: given values for x,y, f(x,y,z) is still a function of z and that I am still allowed to find the root in terms of z when x,y are given. For example: f(x,y,z) = x+y + (x^2-z), given x=1,y=3, f(1,3,z)= 1+3+1-z is a function of z, and then I can use R to find the root z=5. Thank you. -Chee [[alternative

Dear All, I would like to generate m positive real numbers c_i, I=1,...,m, such that (1) c_1 + c_2 + ... + c_m=m, (1) after being ordered into c_1 >= c_2 >= .... >=c_m>0, we have that c_m is of the same order of m^(-1/8), when m is sufficiently large. Thanks, -Chee [[alternative HTML version deleted]]

Dear All, For function: rmvnorm{mvtnorm} in (library mvtnorm, not splus2R), if I generate 2 bivariate normal samples as follows: > rmvnorm(2,mean=rep(0,2),sigma=diag(2)) [,1] [,2] [1,] 2.0749459 1.4932752 [2,] -0.9886333 0.3832266 Where is the first sample, it is stored in the first row or the first column? Does this function store samples row-wise or column-wise? Thank

Dear All, I would like to ask a question on how to do overlay plots in each subgraph of xyplot. 1. I did simulations for m=1000, 2500, 5000, 10000, as the sample sizes. 2. for each sample size value m, 4 graphs are generated; each graph contains overlayed comparisons between 4 methods, 3. now I want put them into a 4-by-4 plot by xyplot, i.e., 4 sample size values, each of which has 4 plots.

Dear All, I would like to ask a question on how to incorporate into an R script help information for the user. I vaguely recall that I saw some instructions on an R manual, but am not able to figure them out. Hereunder is the basic setting: 1. I finished writing an R script, my_script.r, that is a function (the function is named "my_func") for a statistical procedure. This means that

Dear All, In Storey's estimator of the proportion of true nulls, the estimator depends on the tuning parameter lamda. Suppose now that an estimator of this proportion has been obtained by the qvalue package, what is the lamda that corresponds to the estimate? How to get this lamda? Thanks, -Chee [[alternative HTML version deleted]]

Dear All, Suppose Z_i, i=1,...,m are marginally identically distributed as a two normal mixture p0*N(0,1) + (1-p0) *N( miu_i, 1) where miu_i are identically distributed according to a mixture and I have generated Z_i one by one . Now suppose these m random variables are jointly m-dimensional normal with correlation matrix M= (m_ij). How to proceed next or how to start correctly ? Question:

Dear All, I would like to know how to do the following: 1. suppose I have x values from the ordered from 0, 0.5, 1, and would like to label these three points on the x-axis. 2. However, R labels them as 0.0, 0.5, 1.0. But I wan5 them to be 0, .5, 1, since the former way uses limited space of a multi-subgrap plot by adding extra zeros Thank you, Chee [[alternative HTML version deleted]]

Dear All, I would like to seek for help on this issue: 1. I set par(mfrow=c(2,2)), hoping to plot 4 subgraphs in a whole graph 2. Each subgraph has its own x,y axes and each has x-axis label and y-axis label 3. moreover, subgraphs in the left column of the whole graph are all 3D, and have z axes and labels for z axes 4. subgraphs in the right column of the whole graph are all 2D 5. In each

Hi, After reading some of the posting in this list, I came across this posting from: From: Peter Dalgaard <p.dalgaard_at_biostat.ku.dk> Date: Fri 29 Oct 2004 - 08:02:40 EST Replying to Re: [R] Errors during make check He described a problem similar to mine, that build 64-bit R (I am building R version 2.0.1) with sunperf library gives a core dump during make check. So I configured my

Dear All, I would like to report a possible mistake in the generic function rzigp in ZIGP package, and seek for help. Yesterday I tried the "rzigp" function, in this way: 1. generate poisson means from a pareto distribuion, 2. generate dispersion for each poisson from gamma distribution 3. generate ratios between poisson means from uniform 4. generate dispersed poisson counts by

Just installed samba-3.0.0beta1 and trying to use swat. When clicked on help on a particular paramenter (any one), a new window pops up and supposingly it will load smb.conf.5.html and then jumps to the label of the parameter. But now it doesn't work. smb.conf.5.html is no longer complete. Checked under docs/docbook/smbdotconf and found that smb.conf.5.xml supposingly should include files

I am looking for a way to do multiresponse parameter estimation for a problem that is quite similar in form to the alpha-pinene example in Bates and Watts (1988), p. 147, where the model is a set of ODEs. I am specifically interested in obtaining the parameter covariance matrix, as well. Which package(s) in R would I use to do such a thing? Mark Schaffer Dept. of Chemical Engineering

Hello, Could anybody tell me how to find the index for a value which I define in an array? ie. c<-matrix(1:9,1,9) #A 1x9 matrix filled with numbers 1 - 9 in order I want to know the index for the value 5. Thanks in advance, George Chen

Hi: setup: Samba-3.0alpha24 with OpenLDAP 2.1.20 client: XP smbclient test with an LDAP user works But I cannot join the domain PDC (with root, uidNumber = 0): 1. if I set passwd backend = ldapsam_compat, I got: Access Denied 2. if I set passwd backend = ldapsam_compat, ldapsam, I got No mapping between account names and security IDs was done. The machine account was created on the fly