search for: proba

...I get two different curves when I plot the estimated density and the density given by pnbinom. Shouldn't it be the same? This is only the case, I think, if I use the parameter size = 1. I just tried it for size = 2 and the result seemed correct.. Anyways, this is what I tried: #parameters proba<-0.001039302; s<-1; #random numbers with this distribution dom<-rnbinom(10000,s,proba); #first plot plot(ecdf(dom),do.points=FALSE); #second plot x<-seq(0,8700,0.5); lines(x,pnbinom(x,s,proba),col="red");

Hi list, I have problems with the as.numeric function. I have imported probabilities from external data, but they are classified as factors as str() shows. Therefore my goal is to convert the colum from factor to numeric level with keeping the decimals. I have googled the problem for a while now and kept to several advices like http://cran.r-project.org/doc/FAQ/R-FAQ.html...

...newbie with R. I have vectors of the form v = c(1,1,-1,-1,-1,1,1,1,1,-1,1) (longer than this one of course). The elements are only +1 or -1. I would like to calculate : - the frequencies of -1 occurences after 2 consecutives -1 - the frequencies of +1 occurences after 2 consecutives +1 It looks probably something like : Proba( Ut+2=1 / ((Ut+1==1) && (Ut==1))) could someone please give me a little hint about how i should/could begin to proceed ? Thanks (Thanks also to the R creators/contributors, this soft seems really great !)

...choixg <- rep(0,length(conc)) for (i in 1:length(conc)){ choixg[i] <- (k + conc[i])^2/((k+conc[i])^n + (k+1)^n) } return(choixg) } simulation(5,1) plot(conc, choixg, main ="fonction de choix", col= "blue", pch=20, xlab = " concentration", ylab="proba de choisir la gauche") ########################################################## Lassana KOITA Service Technique de l'Aviation Civile (STAC) Direction G??n??rale de l'Aviation Civile (DGAC) Tel: 01 49 56 80 60 Fax: 01 49 56 82 14 http://www.stac.aviation-civile.gouv.fr

...#if no summaries are entered, they are computed if(missing(n)) { n <- length(controls) mean.c <- mean(controls, na.rm=na) sd.c <- sd(controls, na.rm=na) } dl <- n-1 #degrees of freedom of the test #t.test computation t.obs <- (patient-mean.c) / (sd.c*(((n+1)/n)^0.5)) proba.onetailed <- 1-pt(abs(t.obs), df=dl) rar <- pt(t.obs, df=dl) #point estimate of the rarity #confidence intervals computation on the rarity (Crawford & Garthwaite, 2002, Neuropsychologia) c <- (patient-mean.c)/sd.c #finding the non central parameter of t distributions f <- f...

...ta of the second kind, eg. pgbeta2 <- function(quint,b,a,p1,p2) { integrate(function(x) {exp(log(a)+(a*p1-1)*log(x)-(a*p1)*log(b)-log(beta(p1,p2))-(p1+p2)*log(1+(x/b)^a))},0,quint)$value } but I'm facing problems with the quantile function. I tried something like qgbeta2 <- function(proba,b,a,p1,p2) { optimize(function(z) {(proba-pgbeta2(z,b,a,p1,p2))^2},lower=0, upper=10^200) } but it doesn't work. I tried with other non linear optimization command like optim but it is apparently not the solution. Any idea ?

...x <- NULL; y <- NULL simResult <- simulation (5,n) conc <- simResult$conc choixg <- simResult$choixg prin("n"); print(c(k=5, n)) points(conc, log10(1-choixg), main ="fonction de choix", col= n, pch=20, lwd = 3, xlab = " concentration", ylab="proba de choisir la gauche", type="l") } ######################################################################### Lassana KOITA Etude S??curit?? et Exploitation a??roportuaires / Aerodrome Safety & Statistical analysis Service Technique de l'Aviation Civile (STAC) / Civil Avia...

Happy New Year, r-helpers! I am using lowess 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

...conc)) return(list(choixg=choixg, conc = conc)) } #choixg <- simulation(5,2) #mydf<-simulation(5,1) simResult <- simulation (5,5) #plot(mydf$conc, log10(1-mydf$choixg), main ="fonction de choix", col= #"blue", pch=20, #xlab = " concentration", ylab="proba de choisir la gauche") plot(simResult$conc, log10(1-simResult$choixg), main ="fonction de choix", col= "blue", pch=20, lwd = 3, xlab = " concentration", ylab="proba de choisir la gauche") lines(simResult$conc, log10(1-simResult$choixg), col= "red&...

I encountered this error/warning (I don't know) and I saw in http://lists.samba.org/archive/samba/2005-May/105124.html that jerry was aware of the problem. Should I worry about it or ignore this message ? Thanks.

...nes in 'Matrix-class.Rd' below and do not overfill the line, while in the second case they do. > library(Matrix) > promptClass("Matrix") A shell of class documentation has been written to the file 'Matrix-class.Rd'. > promptMethods("[<-",file="proba.Rd") A shell of methods documentation has been written to the file 'proba.Rd'. The syntax is very similar in both cases and uses \item{}{}, except that in the case of class the first argument is the function name (printed in bold). The signature is put in the second argument togeth...

...using the --csum-length option. Any value less than or equal to 16 is valid. Note that if you use this option then you run the risk of ending up with an incorrect target file. The risk with a value of 16 is microscopic and can be safely ignored (the universe will proba- bly end before it fails) but with smaller values the risk is higher. Current versions of rsync actually use an adaptive algorithm for the checksum length by default, using a 16 byte file checksum to determine if a 2nd pass is required with a longer bl...

Dear R-users, I have problems with printing Central European (CE) characters (e.g. {\H o} and {\H u}) in graphics made by R. After executing the following code I can see the CE characters with gv on my screen, but I cannot print them on a postscript printer. > postscript(file="proba.ps",enc="ISOLatin2.enc") > plot(0:15,0:15,type="n") > grid(15,15,lty=1) > for(i in c(32:127, 160:255)){x <- i %%16;y <- i%/% 16;points(x,y,pch=i)} There were 50 or more warnings (use warnings() to see the first 50) > warnings() Warning messages: 1: font wi...

...ntext, const std::vector<double>& inX, const std::vector<double>& inF); void RSymbReg(SEXP RPopSize, SEXP RNbrGen, SEXP RNbrPartTournament, SEXP RMaxDepth, SEXP RMinInitDepth, SEXP RMaxInitDepth, SEXP RInitGrowProba, SEXP RCrossoverProba, SEXP RCrossDistribProba, SEXP RMutStdProba, SEXP RMutMaxRegenDepth, SEXP RMutSwapProba, SEXP RMutSwapDistribProba, SEXP RSeed); void RSymbReg(SEXP RPopSize, SEXP RNbrGen, SEXP RNbrPartTournament, SEXP RMaxDepth, SEXP RMinInitDepth, SE...

I am using the rw0901 version of R and I fail in reading large files. Actually, I really dont know how to extend the memory. Does somebody knows? Louis FERRE D?partement de Math?matiques et Informatique Equipe GRIMM |Laboratoire de Stat et Proba Universit? Toulouse Mirail |Universit? Paul Sabatier 5, all?es Antonio Machado |Ade:lferre at cict.fr 31058 Toulouse Cedex | Tel:05/61/50/46/11 | Adel:loferre at univ-tlse2.fr |http://www-sv.cict.fr/lsp/Ferre/index.html -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-....

..."rpanel" ¿Cómo puedo introducir una etiqueta utilizando rpanel? creo un panel de la siguinte forma, pero me gustaría introducir una etiqueta antes de los botones y no veo como puede hacerse anchoVentana<-1024 altoVentana<-600 panel<-rp.control(title="probando rpanel",size=c(anchoVentana-70,altoVentana-70)) rp.button(panel,title=" PRUEBA ", pos=c(150,300,200,100)) intento introducir un botón de la siguiente forma label1<-tklabel(panel,text="proba etiqueta")...

Buenas tardes Les consulto los siguiente, dentro de unos datos hay registros faltantes, estos son NA. Me gustaría poder avisar cuáles son los registros no tomados en cuenta porque son NA. Como no puedo reproducir el error, copio y pego un script R que es la forma que yo razone, la cuál me funciona, y un trozo de la captura de pantalla donde aparece mi error, es decir todo NA. # NA's

...ote server as part of the compilation process. The latency for doing this is rather high, so I want to allow users to do a 'make -j' to parallelize the build, in order to hide the network latency. I guess for now I'll tell users to run 'make -j N' with N < 6 or so (which is probably not a bad idea anyway). But I can't imagine I'll be the last person to pound on ssh like this... ------- You are receiving this mail because: ------- You are the assignee for the bug, or are watching the assignee.

Hi, we use OpenSSH_3.8.1p1 and we would like to change the complete printer Spool communication from Unix r-Commands to ssh "OpenSSH". Sometimes, we have problems, because of very high connection set-ups in extremely short time intervals, more than 10 open sessions. The connection set-up is partly declined with the error message: ssh_exchange_identification: Connection closed by remote

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi folks. Please, someone confirm if i'm right. if i'm wrong, please forgive-me. - -------------------------------------------------------------------------------------------------------------- I've developed a little tool to stress test tcp connections.( Sending syn and answer ack-syn ). that simuates a real tcp connection. (