similar to: Data analysis: normal approximation for binomial

Thank you for the suggestion, Bill. The result is not quite what I would like. Here's sample code for you or anyone else who may be interested: Al1 = c('A','C','C','C') Al2 = c('G','G','G','T') Freq1 = c(0.0078,0.0567,0.9434,0.9908) MAF = c(0.0078,0.0567,0.0566,0.0092) m1 = data.frame(Al1=Al1,

Hello, All! How might one accomplish this using the by() function? m1 is a data frame. # populate column "m1$major_allele" for ( i in 1:length(m1$major_allele)) { if ( m1$Freq1[i] == m1$MAF[i]){ m1$major_allele[i] = m1$Al1[i] } else{ m1$major_allele[i] = m1$Al2[i] } } Jim [[alternative HTML version deleted]]

Hello R useRs, I have a function which returns a list of functions : freq1 <- function(x) { lev <- unique(x[!is.na(x)]) nlev <- length(lev) args <- alist(x=) if (nlev == 1) { body <- c("{", "sum(!is.na(x))", "}") f <- function() {} formals(f) <- as.pairlist(args) body(f) <- parse(text = body) namef <-

Hi, I have loaded a dataset in R : data = label freq1 freq2 news 54 35 fun 37 21 milk 19 7 food 3 3 .... etc And I have a vector flist<-c("fun","food") Now I want to use the vector 'flist' for selecting these values from 'data' so that I get the following dataset : label freq1 freq2 fun 37 21 food

I was wondering how do I actually see what's inside a function, say, density of t distribution, dt()? I know for some, I can type the function name inside R and the code will be displayed. But for dt(), I get > dt function (x, df, ncp, log = FALSE) { if (missing(ncp)) .Internal(dt(x, df, log)) else .Internal(dnt(x, df, ncp, log)) } <environment: namespace:stats> I

Dear masters, I have data that looks like this: #val Freq1 Freq2 0.000 178 202 0.001 4611 5300 0.002 99 112 0.003 26 30 0.004 17 20 0.005 15 20 0.006 11 14 0.007 11 13 0.008 13 13 ...many more lines.. Full data can be found here: http://dpaste.com/173536/plain/ What I intend to do is to have a cumulative graph with "val" as x-axis with "Freq1" & "Freq2" as

Good day, All, Is there any way to maintain the number of decimal places in the type of situation below? I would like to maintain the number of decimal places in 0.667, despite the fact that its column-mates have a fourth decimal place. Thank you for your time. Jim dat.txt contents: MARKER ALLELES FREQ1 RSQR EFFECT2 STDERR CHISQ PVALUE rs6599753 C,T

Hi, I have a grouped data set and would like to calculate weighted proportions for a large number of factor variables within each group member. Rather than using dplyr::count() on each of these factors individually, the idea would be to do it for all factors at once. Does anyone know how this would work? Here is a reproducible example: ############################################################

Dear List Is there a package for leapfrog plotting (Hamiltonian Monte Carlo estimation) in R? I tried the actual "LEAPFrOG" package which doesn't actually give the plot like this one? http://xianblog.files.wordpress.com/2010/09/hamilton.jpg How doe one plot this in R? So, there semi-circle and dots on that semi-circle. I don't think curve() or plot() would produce such plot.

I am wondering if anyone could help me figure out how to turn down the volume on all the dial tones, indications, etc.. and especially the call-waiting CHIRP! I don't want to change the txgain and rxgain because they are working at levels that I would like. However, when voice conversations and voicemail recordings are at good levels then the dial tones, busy tones, etc are way too loud.

I am trying to configure a VegaStream 50 FXO to work with asterisk. The problem that I am having is that the VegaStream does not support incoming registration from asterisk. VegaStream only allows outbound registration. My question is does asterisk allow incoming registration from an FXO? If yes how? Or better yet, has anybody been able to make the VegaStream FXO work with asterisk? According

Dea'R' helpers I have following data - prob = c(0.1, 0.2, 0.3, 0.4, 0.5) frequency = c(100, 75, 45, 30, 25) no_trials = c(10, 8, 6, 4, 2) freq1 = rbinom(frequency[1], no_trials[1], prob[1]) freq2 = rbinom(frequency[2], no_trials[2], prob[2]) freq3 = rbinom(frequency[3], no_trials[3], prob[3]) freq4 = rbinom(frequency[4], no_trials[4], prob[4]) freq5 = rbinom(frequency[5],

Hello, Everyone. Does anyone know how to create a Nightingale’s Rose chart by using R? Hopefully, the graph could be displayed like this: http://mbostock.github.com/protovis/ex/crimea-rose.html Thanks a lot. Kind regards, Henry [[alternative HTML version deleted]]

Hi, First time using the optim(), can someone please tell me what I am doing wrong? The error looks like this Error in .Internal(pnorm(q, mean, sd, lower.tail, log.p)) : 'sd' is missing An example of the error dat = c(20, 19, 9, 8, 7, 4, 3, 2) dat_mu=mean(dat) dat_s=sd(dat) max.func = function(dat, mu, sd) { pnorm(dat, mu, sd) } optim(fn=max.func, dat=dat, par=c(mu=dat_mu,

Hi there, Besides the problem of Voicetronix dialing too early before the carrier gives a dial tone, there also appears to be issues with the audio quality when more than 1 channel is utilized.

Dear R-expert, For histogram function, can we get the table of bin and frequency like in excel, together with the histogram? Therefore, we can check the number of data included. Thank you so much for your attention and help. ____________________________________________________________________________________ [[elided Yahoo spam]]

Hi, I have a couple of observations, a question or two, and perhaps a suggestion related to the plotting of density on the y-axis within the hist() function when freq=FALSE. I was using the function and trying to develop an intuitive understanding of what the density is telling me. After reading through this fairly helpful post:

i have created a package which allows R to play sounds the package is available from http://sunsite.univie.ac.at/~neuwirth/Rcsound/ both as zipped binary for windows and as installable tar.gz preequisite: you have to install csound on yout machine. the necessary links are at http://www.csound.org csound is a very portable program for creating sounds, and using it as out sound creation tool we

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hello, I'm testing out dovecot with postfix. I'm running postfix 2.5.2 with dovecot version 1.1.1 with dovecot sieve plugin version 1.1.5. I have everything compiled and working except for the sieve plugin. I'm now trying to test out the sieve plugin and having no luck in getting a simple auto vacation reply to work. Here is my