similar to: Quantile question

Dear All, I am trying to calculate the Hardy-Weinberg Equilibrium p-value for 42 SNPs. I am using the function HWE.exact from the package "genetics". In order not to do a lot of coding "by hand", I have a for loop that goes through each column (each column is one SNP) and gives me the p.value for HWE.exact. Unfortunately some SNP have reached fixation and HWE.exact requires a

Hi, I need to support port ranges in tc filter rules. I know how to formulate the rule but , I am not able to understand how to calculate the mask value for a perticular range so as to segregate the port values that lie within this range . I got the following sample "tc filter add dev eth1 parent 1:1 protocol ip prio 10 u32 match ip sport 0x1ae0 0x1ff0 flowid 1:10 This rule will match all

Dear All I have problem about calculte variance use R. I use trimm mean by mean(x, trim=.05) and I want to calculate variance with trimm=.05. How I can do this. Many Thanks. Jumlong -- Jumlong Vongprasert Assist, Prof. Institute of Research and Development Ubon Ratchathani Rajabhat University Ubon Ratchathani THAILAND 34000 [[alternative HTML version deleted]]

Dear list, Hello to everybody, I´m interested in finding a package for statistical analysis of binary data, I have a matrix with the following structure: Case1 Case2 Case 3 ....... CaseX Control1 Control2 Control3 ....... ControlY Pep1 1 0 1 1 0 0 0 1 Pep2 1 1 1 1 1 0 0 1 Pep3 0 1 1 1 1 0 0 1 ... Pepz 1 1 1 1 1 0 0 1 I would like to know some R package to calculte similarity matrix or

Hi everyone I have a little problem with R. I built a lineal regression equation using stepAIC function in both directions. Once I get this formula (lineal regression), I would like to save in a txt file data refered to p-value, r-squared, coefficients,...from the choosen model previously using stepAIC and p-value,r-squared, coefficients,...from each predictor. I try to use function

Dear R family I have a following question. Suppose I have a matrix as follows, for instance: tau= 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 1 0 0 0 0 I want to calculate (-m) power of tau, for example, m=893. When I run tau^2, the outcome is just tau. Any help would be appreciated. thanks in advance. best moohwan

Greetings, I am having a problem with DescTools::Quantile (a function computing quantiles from weighted samples): # these sum to one probWeights = c( 0.0043, 0.0062, 0.0087, 0.0119, 0.0157, 0.0204, 0.0257, 0.0315, 0.0378, 0.0441, 0.0501, 0.0556, 0.06, 0.0632, 0.0648, 0.0648, 0.0632, 0.06, 0.0556, 0.0501, 0.0441, 0.0378, 0.0315, 0.0257, 0.0204, 0.0157, 0.0119, 0.0087,

Dear R sages, I used the function rbind to combine a matrix (M) and a vector (Fert) to get a new matrix (A). This was fine. The issue is however, that the new matrix A has as its row names the name of the vector Fert, even though I set teh new vector A to have dimnames=NULL. See short code below fyi. *Fert<-c(0,1,5) * *M <- matrix(0, 2, 3) diag(M) <- c(0.3,0.5) * *A<-

Hi, pchisq -> distribution function dchisq -> density function pval is the area under the curve, to calculte it you use distribution function which is the integral of density function. See: http://www.itl.nist.gov/div898/handbook/eda/section3/eda362.htm http://mathworld.wolfram.com/DistributionFunction.html f(x) density function F(x) distribution function =Pr(X<x)= integral(f(x))

I have a following matrix and wish to define a variable based the variable A=matrix(0,5,5) A[1,]=c(30,20,100,120,90) A[2,]=c(40,30,20,50,100) A[3,]=c(50,50,40,30,30) A[4,]=c(30,20,40,50,50) A[5,]=c(30,50,NA,NA,100) > A [,1] [,2] [,3] [,4] [,5] [1,] 30 20 100 120 90 [2,] 40 30 20 50 100 [3,] 50 50 40 30 30 [4,] 30 20 40 50 50 [5,] 30 50 NA NA

I'm a PhD student and I'm working with covariance function. I'm interested to know if exist some packages in R to calculate and plot the bidimensional Autocovariance Function. the input matrix is a matrix that describe a spatial location over a 2-D space and I want to use it in the same way I can use a time serie in the 1-D acf. Thanks, Nicola.

I'm a PhD student and I'm working with covariance function. I'm interested to know if exist some packages in R to calculate and plot the bidimensional Autocovariance Function. the input matrix is a matrix that describe a spatial location over a 2-D space and I want to use it in the same way I can use a time serie in the 1-D acf. Thanks, Nicola.

Hi, I have a Community (COM) composed of 6 species: A, B, C, D, E & F. The density of my Community is thus (Eq.1): dCOM = dA + dB + dC + dE + dF I would like to calculate and plot a linear regression between the density of each of my species and the density of the whole community (illustrating how the density of each species varies with variations of the whole community). For example, I would

Dear R users, I would like to create a filled.contour plot with colour bar and overlay a map. I need to create several of these representing variables with different spatial extents for comparison. When I try to fix the plot area size and axes limits as follows: par(fig=c(0,1.0,0.3125,1.0) filled.contour(x=seq((320/2)-180.5,(480/2)-180.5,len=nrow(oct94)),

Hello, I saw that Asterisk don't calcultate fine the ANSWEREDTIME. I want that when ANSWEREDTIME =~ 5.6 become 6 and if =~10.3 become 10 because, now, if ANSEREDTIME =~ 15.9, it become 15! it isn't correct How can I have a rounded ANSWEREDTIME ? Where have I to manipulate the sources? thank you -- Francois -------------- next part -------------- A non-text attachment was scrubbed...

Hello Dear, I have a function, like z=f(x,y), and try a surface plot with this function. But, on the reference of "wireframe" requires data option, so I generated x and y, and computed z with them. But, still I have a problem to draw a surface plot. The code and errors are ################################################## mle_beta0=64.43707; mle_beta1=-24365.16; # generating for

I have a following matrix and wish to define a variable based the variable A=matrix(0,5,5) A[1,]=c(30,20,100,120,90) A[2,]=c(40,30,20,50,100) A[3,]=c(50,50,40,30,30) A[4,]=c(30,20,40,50,50) A[5,]=c(30,50,NA,NA,100) > A [,1] [,2] [,3] [,4] [,5] [1,] 30 20 100 120 90 [2,] 40 30 20 50 100 [3,] 50 50 40 30 30 [4,] 30 20 40 50 50 [5,] 30 20 NA NA

It looks like a homework assignment. It also looks like you didn't read the documentation carefully enough. The 'len.out' argument in seq is solely for specifying the length of a sequence. The 'quantile' function omputes the empirical quantile of raw data in the vector 'x' at cumulative probabilit(y)(ies) given in the weights' argument, with interpolation I'm

Hi, I'm confused about how to specify random and fixed factors in an aov() term. I tried to reproduce a textbook example: One fixed factor (Game, 4 levels) and one random factor (Store, 12 levels), response is Points. The random factor Store is nested in Game. I tried > str(kh.df) `data.frame': 48 obs. of 4 variables: $ Subj : Factor w/ 48 levels

Hello. I have n files in a directory: file1, ..., filen. I read them with the following commands: list=scan(file="list",what=list(nom="")) # in the file list, I have all the filenames. n=length(list[[1]]) for (i in 1:n){ aux <- paste("p",i,sep="") assign(aux, as.matrix(read.table(list[[1]][i]))) } R creates the matrices p1,p2,...,pn. I want