search for: calcuate

Hi, folks, Please let me address the problem by the following codes: first=c('u','b','e','k','j','c','u','f','c','e') second=c('usa','Brazil','England','Korea','Japan','China','usa','France','China','England') third=1:10

Hi, I want to analyse DNA sequence data (mtDNA) in R as in calculate Fst, Heterozygosity and such summary statistics. Package Adagenet converts the DNA sequence into retaining only retaining the polymorphic sites and then calcuates Fst.. but is there any other way to do this? I mean analyse the DNA sequence as it is.. and calculate the statistics? Thanks! [[alternative HTML version deleted]]

Hi Firstly , happy Christmas to R-Help! Secondly, I wonder if anyone can help me with the following query: I am trying to reproduce some explicit probability calculations performed in APPL (a Maple extension for computational probability). For instance, in APPL, to compute the probability that the sum of 10 iid uniform variables [0,1] will be between 4 and 6, (i..e Pr( 4 < \sum_{i=1}^{10}X_i

...surements via appendinng the 6 measurements each time to the current set? I couldn't find anything in Intro to R on appending. cheers, Dave ps - please respond directly to afshar at miami.edu creatine.function.new = function(delta.0.Y.0, gamma, comp.LIS.frm, comp.CAND.frm) { ## function to calcuate the delta.i, i.e. the percent ## leftover ## gamma = rate of Cr going into bucket, e.g., mg/hr ## delta.0.Y.0 = product of delta.0 and Y.0 at baseline ## Y.1 = delta.0.Y.0 + gamma delta = numeric(6) delta.patient = numeric(24) delta.patient.comb = numeric(144) ## for (k in 1:24) { ## each patient...

Dear R community, I am using 6 variables to test for an effect (by linear regression). These 6 variables are strongly correlated among each other and I would like to find out the number of independent test that I perform in this calcuation. For this I calculated a matrix of correlation coefficients between the variables (see below). But to find the rank of the table in R is not the right

Hello, I am using tapply to pull out data by the day of week and then perform functions (e.g. mean). I would like to have the number of values used for the calcuation for the functions, sorted by each day of week. A number of entries in any given column are NAs. I have tried the following code and simple variants with no luck. for (i in 1:length(a[1,])){ x<-tapply(a[,i],a[,1],mean,

Dear expeRts, is it possible to carry out calculations between different foreach() calls? As for nested loops, you want to carry out calcuations not depending on the inner loop only once and not for each iteration of the innermost loop. Cheers, Marius library(foreach) foreach(i=1:3) %:% foreach(j=1:2) %do% { i <- i+1 print(paste(i,j)) } foreach(i=1:3) %:% i

...aving..). When I search the Actions table by populating "Date Closed" in my Find view with a date range (ex, 1/1/2007...1/31/2007) i return no records. When I just use one value (ex, 1/1/2007) i return no records. When I view my list of actions, the "Date Closed" is correctly calcuated, and record DO exist that meet this criteria. When I search on the native "Date Created" field (auto entered based on when the record is created) it works fine. When I search on a related field that is Text, it works fine. When I search on a related calculated field that is text, it wo...

Hi, From the file HowToUseJIT.cpp, I learned how to insert some calcuation IR code like Add/Sub/Mul etc by using Instruction *Add = BinaryOperator::createAdd(One, ArgX, "addresult", BB); By following this way, it works well when I insert some IR code whose operand is integer type like IntTy, however, when I tried to insert similar thing whose operands are float point, I got the

...a coefficients of a gene signature consisting of xx genes and the gene expression is also provided of those xx genes. If I try to use "coxph()" function in R using just the gene expression data alone, the beta coefficients and coxph$linear.predictors will change and I need to use the pre-calcuated linear predictor not re-computed using coxph() function. The reason is I need to compute a quantity that uses as it's input the coxph() output but I need this output to be pre-calculated beta-coefficients and linear.predictor. Any one can show me how to do this in R? Thanks a lot. [[altern...

hi,r project I am a new user who has some questions about r tool,I have some original data file in netcdf format which is a famous geometry data file format, The file contains very large weather information from every zone of the world,I want to extract only a part of them to calcuate, like the information in a city or the information between exact latitute range and longtitude range, can r tool help me to do that or in some point r tool is a statitcal software which focus on mathmatical caculation and the user have to translate the special data files into some data structure l...

I tried to use gpuCor function in the gputools package of R to calcuate the pairwise correlations of a matrix of 40,000 columns. Becuase there would be memory issues if I use the whole matrix at a time, I splitted the matrix into submatrix of 10,000 columns and then calculate the pairwise correlation of different submatrices. There are altogether 4 submatrices, so I n...

I need to calcuate the cumulative probability for the Natural Exponential Family - Hyperbolic secant distribution with a parameter theta between -pi/2 and pi/2. The integration should be between 0 and 1 as it is a probability. The function "integrate" works fine when the absolute value of theta is not to...

Hello, I can't find a package which calculates Kendall's tau-c. There is the package Kendall, but it only calcuates Kendall's tau-b. Here is the example from ttp://www2.chass.ncsu.edu/garson/pa765/assocordinal.htm. cityriots <- data.frame(citysize=c(1,1,2,2,3,3), riotsize=c(1,2,1,2,1,2), weight=c(4,2,2,3,0,4)) cityriots <- data.frame(lapply(cityriots,function(x) rep(x,cityriots$weight))) xtabs(~ rio...

...N-p)*(N-p-2)) + N*(p+1)/((N-p)*(N-p-2)) CAIC <- -2*cond.loglik + 2*K.corr return(CAIC) } Reference: http://r.789695.n4.nabble.com/Using-Conditional-AIC-with-lmer-td847899.html However, in this function the function hatTrace of the lme4 package was used to calcuate rho (line: rho <- hatTrace(model) ) Sadly, R tells me that this function could not be found. A search on the internet revealed that this function was dropped out of the new package release of lme4. Can anyone help me with this? Is there an alternative function that gives me the same result...

My goal is simple: calcuate GC content of each sequence in a list of nucleotide sequences. I have figured out how to vectorize, but all my attempts at memoization failed. Can you show me how to properly memoize my function? There is a StackOverflow post on the subject of memoization, but it does not help me: http://stackove...

Dear all - I have followed the thread the reply to which was lead by Thomas Lumley about using pweibull to generate fitted survival curves for survreg models. http://tolstoy.newcastle.edu.au/R/help/04/11/7766.html Using the lung data set, data(lung) lung.wbs <- survreg( Surv(time, status)~ 1, data=lung, dist='weibull') curve(pweibull(x, scale=exp(coef(lung.wbs)),

Hi all, I have 2 questions: 1)How do I calculate the mean on an imported txt file? I've imported the file below and that's what it looks like imported. How do I then calcuate the mean, median, or mode on the column LeafArea using the desktop R package? Any help would be greatly appreciated!! Thanks, Nat LeafType Leaflets LeafArea ShapeRatio LeafWeight LeafThickness 1 1 3 0.12 0.12 0.21 0.00 2 1 3 0.17...

Hi, In the code below I create a small data.frame (dat) and then cut it into different groups using CutList. The lists in CutList allow to me choose whatever columns I want from dat and allow me to cut it into any number of groups by changing the lists. It seems to work OK but when I'm done I have a variable number of data frames what I need to do further operations on and I don't know

...wess to smooth a scatter plot, xx<-lowess(xinput,f=.04) #defaults for other args followed by turnpoints(xx$y) #defaults for other args I plot the smoothed result as well as turnpoints (using yy$tppos) on top of raw data plot. Result is exactly as expected, graphically. For another purpose, I calcuate the difference between turnpoints (representing time intervals between turnspoints in my applicaiton), e.g. aduration[j]<-yy$tppos[i+1]-yytppos[i] This also appears to work as expected, HOWEVER, a typical list of such differences looks like: 22,33,22,11,33,44,33,33,11,33,44,66,33,22,............