Re: [R] Self-starting nonlinear power law function

Best regards, Tiago Henriques -- View this message in context: http://r.789695.n4.nabble.com/Self-starting-nonlinear-power-law-function-tp4681626p4712967.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] How to double integrate a function in R

I would like to express my gratitude for the great help given by David and Hans regarding my last query. Thank you very much for your time, people. All the best, Tiago --- Hello, R users! I am trying to double integrate the following expression: # expression (1/(2*pi))*exp(-y2/2)*sqrt((y1

[R] How to double integrate a function in R

(x[2]-x[1])))} adaptIntegrate(fun, lower = c(0,0), upper =c(5, 6), tol=1e-8) However, I don't know how to constrain the integration so that y2>y1>0. Any ideas? Tiago __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

Re: [R] How to compute a P-value for a complex mixture of chi-squared distributions in R

Thank you very much, Rui and Peter, for you detailed and helpful tips! It worked like a charm! I would spend more two weeks (or more) to figure out that by myself. Cheers! Tiago > Hello, > > No, nothing wrong. (I feel silly for not having noticed it.) In fact not > only it's

[R] How to compute a P-value for a complex mixture of chi-squared distributions in R

is not efficient. Are you aware of smart approach? All the best, Tiago __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, m

[R] How to compute a P-value for a complex mixture of chi-squared distributions in R

is not efficient. Are you aware of smart approach? All the best, Tiago __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, m

[R] How to create a correct matrix in R

Hello Rlisters! In my codes, I need to import a matrix: v <- read.table("/home/tiago/matrix.txt", header=FALSE) v<-as.matrix(v) v V1 V2 V3 V4 V5 V6 [1,] 1. -0.89847480 -0.73929292 -0.99055335 -0.0451446

[R] Creating cluster with dbscan from bigmemory

ever dealt with one or both of these problems, please let me know how you achieved it. Thanks. Tiago Cunha -- Tiago Daniel Sá Cunha Faculdade de Engenharia da Universidade do Porto Skype: tiagodscunha Facebook: http://www.facebook.com/TiagoDSCunha Email alternativo 1: ei08...@fe.up.pt Email

[R] aligning image with bar plot in the same plot

there is an easier way to do it. Thank you so much! Tiago # code ## matrix to plot matI <- matrix(c(rep(1:3, 100), rep(4:6, 100), rep(7:9, 100)), nrow=3, byrow=T) ## plot with three panels layout(matrix(1:3, nrow=1), widths=c(1,1,1)) # margins mar1 <- c(1,1,1,1) par(ma

[R] help with time data - R weird behaviour across machines

9:57" "2011-03-27 00:59:58" "2011-03-27 01:00:01" [4] "2011-03-27 01:00:02" > is.na(todelete) [1] FALSE FALSE FALSE FALSE > as.numeric(todelete) [1] 1301180397 1301180398 1301180401 1301180402 any suggestions about what might be going on here would be greatly a

Re: [R] How can I access information stored after I run a command in R?

I would like to thank Justin and Matthias for their very helpful on my query (see it below). All the best, Tiago Dear all, Supposed I run the following command: ### #install.packages("Rassoc", dependencies=TRUE) library("Rassoc") ca=

[R] How can I access information stored after I run a command in R?

) However, instead of typing the value above, I would like to replace it by the macro (scalar, local) that has the accurate p-value. thanks in advance for your help. Tiago __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

Re: [R] How to compute the P-value for a mixture of chi^2 distributions

Thanks for your very constructive and helpful tips, Peter and Ted! -- View this message in context: http://r.789695.n4.nabble.com/How-to-compute-the-P-value-for-a-mixture-of-chi-2-distributions-tp3591276p3595373.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] How to compute the P-value for a mixture of chi^2 distributions

thank you very much for your advice, Thomas! Extremely helpful! Cheers! Tiago -- View this message in context: http://r.789695.n4.nabble.com/How-to-compute-the-P-value-for-a-mixture-of-chi-2-distributions-tp3591276p3591365.html Sent from the R help mailing list archive at Nabble.com

[R] How to compute the P-value for a mixture of chi^2 distributions

Hello, The test I am working on has an asymptotic 0.5*chi2(1)+0.5*chi2(2) distribution, where numbers inside parenthesis stand for the degrees of freedom. Is is possible to compute quickly in R the cumulative distribution of that distribution? Thanks in advance. Tiago -- View this message in

[R] Error when creating 2 R instances on the same context.

v); // create an embedded R instance and I'm using R 2.10 on a Ubuntu OS. Any help will be aprecciated even pointing at wich documentation I can find help. Thanks in advance, Tiago Zortea __

Re: [R] issue with L-BFGS-B in optim (optim just hangs)

story! Best, Tiago Ravi Varadhan wrote: Hi Tiago, It is hard for me to speculate without knowing more about your problem. Here is what I would suggest, assuming your problem specification and its computer implementation are correct: (1) You may try to terminate the algorithm by specifying a

[R] issue with L-BFGS-B in optim (optim just hangs)

that the value is the boundary defined in the function call. As anyone seen this behavior before? Is it me missing something or is this some bug in "method=L-BFGS-B"? Any suggestions on how to deal with it? Many thanks for any useful f

[R] infer haplotypes phasing trios tdthap

nd it) Thank you very much, Tiago Magalhães example (ped file with pedigrees) 9 100 102 101 1 2 1 1 2 1 2 2 1 2 9 101 0 0 2 1 1 1 2 1 2 2 2 2 9 102 0 0 1 1 2 1 2 1 2 2 1 1 data out: hap.transmit(example) pedidfathermother 9 100102 101 f.tr.1f.tr.2f

[R] Flushing (Reset) 'last.warning'

Dear members of the mailing list, I want to fetch warnings() from a series of prop.test calls. I want to get "none" if no warning is issued, and "warning" if there is a problem. I have looked (and relooked) at options(warn) and warning(), warnings() and 'last.warning' but to no avail. I am a biol