search for: coh

Hello everybody I'm trying to write a simple version of matlab's "quiver". The idea is that I have fluid with velocity defined on a grid. I have a matrix of x-components of velocity and a matrix of y-components and I want to see the overall flow pattern. (I work with 2D fluid mechanics problems). My first-stab function is below: quiver <- function(u,v,scale=1) # first

...Settings\Suzette\Desktop\MyFile\MyProgram.tex" Any advice would be much appreciated, Suzette Suzette Blanchard, Ph.D. Assistant Professor, Dept. of Biostatistics City of Hope National Medical Center and Beckman Research Institute 1500 East Duarte Road, Duarte, CA 91010-3000 sblanchard@coh.org <mailto:sblanchard@coh.org> ph: (626) 256-4673 ext: 64446 "EMF <COH.org>" made the following annotations. ------------------------------------------------------------------------------ SECURITY/CONFIDENTIALITY WARNING: This message and any atta...{{dropped}}

...unction and need to find a good guide source about writing R functions. If you know a good book or website for guidance of R function, please let me know. Thanks, Ding From: Bert Gunter [mailto:bgunter.4567 at gmail.com] Sent: Monday, January 15, 2018 10:11 AM To: Ding, Yuan Chun <ycding at coh.org> Cc: Richard M. Heiberger <rmh at temple.edu>; r-help at r-project.org Subject: Re: [R] consolidate three function into one That is certainly OK, but you can also just use print(ggsurvplot(...)) as your final statement. out <- RFS( ...) would then return the ggsurvplot object *and*...

I'm using spec for multivariate time series, and have a question about "coh". This is probably a very silly question, but I'm a bit stumped. >From the documentation: Column i + (j - 1) * (j - 2)/2 of coh contains the squared coherency between columns i and j of x, where i > j. What's in column 1 of coh? Cheers Jason -- Indigo Industrial Control...

Greetings, I would like to estimate a spectral coherence between two timeseries. The stats : spectrum() returns a coh matrix which estimates coherence (squared). A basic test which from which i expect near-zero coherence: x = rnorm(500) y = rnorm(500) xts = ts(x, frequency = 10) yts = ts(y, frequency = 10) gxy = spectrum( cb...

...ght find useful. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Mon, Jan 15, 2018 at 9:57 AM, Ding, Yuan Chun <ycding at coh.org> wrote: > Hi Richard, > > Thank you so much!! I understand the problem now, I assign a name to the > "ggsurvplot" object and then add print(fig) at bottom of function > definition, now figure gets printed on screen. > > Ding > > # function to generate...

...nt this type of analysis. I tried all the examples which generated a series of fancy plots. But I need to work on the numerical results. I have read that the following info is available through Cross Spectra analysis: *Cross-periodogram, Cross-Density, Quadrature-density, Cross-amplitude, Squared Coherency, Gain, and Phase Shift* I went through a couple of the two-series (bivariate) cross-spectrum analysis examples with R. I also printed out the attributes of the analysis (see the following). I cannot quite match the above quantities with the attributes/features output of cross-spectra analysis...

...", legend.labs = legend.labs, lty=1, lwd=3) print(fig) } -----Original Message----- From: Richard M. Heiberger [mailto:rmh at temple.edu] Sent: Sunday, January 14, 2018 1:34 PM To: Ding, Yuan Chun <ycding at coh.org> Cc: Bert Gunter <bgunter.4567 at gmail.com>; r-help at r-project.org Subject: Re: [R] consolidate three function into one FAQ 7.22 You must print a ggplot object, for example with print(m52.2cluster) For the FAQ, run the line system.file("../../doc/FAQ") in R on your co...

...See data (from R) below Readable: "Df" "Sum Sq" "Mean Sq" "F value" "Pr(>F)" "pre" 1 0.00593 0.00593936 0.7450563 0.401636958677004 "coh" 1 0.04311 0.04311302 5.4082639 0.0344751749542619 "Residuals" 15 0.11957 0.00797169 NA NA Original: "Df" "Sum Sq" "Mean Sq" "F value" "Pr(>F)&q...

...'ve scaled down the other labels using labelcex whch leaves the disparity in sizes that I'm trying to remove. I'd like something like this but with the x axis lables the same size as the y. > mpop=c(50,75,80,100,110,90,85,60,40,20) > fpop=c(50,75,80,100,110,90,85,60,50,30) > coh=c("0-9","10-19","20-29","30-39","40-49","50-59","60-69","70-79","80-89","90-99") > pyramid.plot(mpop,fpop,labels=coh,gap=15,labelcex=0.7,top.labels=c("males","age group","...

I know that City of Heroes/Villains doesn't work in wine right now and I couldn't find much on the wine website that looked like anyone was actively working on it, but I am very interested in playing this game in Linux. I have zero programming skills, but I do have some free time and would be more than happy to test stuff out and give feedback and/or to hunt through forums and google

...icking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Mon, Mar 5, 2018 at 2:03 PM, David Winsemius <dwinsemius at comcast.net<mailto:dwinsemius at comcast.net>> wrote: > On Mar 5, 2018, at 8:52 AM, Ding, Yuan Chun <ycding at coh.org<mailto:ycding at coh.org>> wrote: > > Hi Bert, > > I am very sorry to bother you again. > > For the following question, as you suggested, I posted it in both Biostars website and stackexchange website, so far no reply. > > I really hope that you can do me a grea...

...eople keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Mon, Mar 5, 2018 at 2:03 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Mar 5, 2018, at 8:52 AM, Ding, Yuan Chun <ycding at coh.org> wrote: > > > > Hi Bert, > > > > I am very sorry to bother you again. > > > > For the following question, as you suggested, I posted it in both > Biostars website and stackexchange website, so far no reply. > > > > I really hope that you can...

...(u4/u2^2) df <- df * (N0 / N) ## << since R 1.9.0 bandwidth <- bandwidth * xfreq/N pgram <- pgram[2:(Nspec+1),,, drop=FALSE] spec <- matrix(NA, nrow = Nspec, ncol = nser) for (i in 1L:nser) spec[, i] <- Re(pgram[1L:Nspec, i, i]) if (nser == 1) { coh <- phase <- NULL } else { coh <- phase <- matrix(NA, nrow = Nspec, ncol = nser * (nser - 1)/2) for (i in 1L:(nser - 1)) { for (j in (i + 1):nser) { coh[, i + (j - 1) * (j - 2)/2] <- Mod(pgram[, i, j])^2/(spec[, i] *...

...s crash dialog. I tried with different wine versions (atm 1.1.42 and 1.2) under Debian lenny 32bit (Linux 2.6.32) and Ubuntu Lucid 64bit with different fglrx versions. My graphics hardware is a Radeon HD 4200. Via winetrix I installed d3dx9, gecko, ie6, ie7, vcrun2005, and vcrun2008 and copied my CoH installation from Windows 7 RC 64bit (with/without using CoH profile from Win7). Any hints? TIA Martin

...statistical site. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Fri, Mar 2, 2018 at 10:34 AM, Ding, Yuan Chun <ycding at coh.org<mailto:ycding at coh.org>> wrote: Dear R users, I need to analyze data generated from a partial two-by-two factorial design: two levels for drug A (yes, no), two levels for drug B (yes, no); however, data points are available only for three groups, no drugA/no drugB, yes drugA/no dru...

> On Mar 5, 2018, at 8:52 AM, Ding, Yuan Chun <ycding at coh.org> wrote: > > Hi Bert, > > I am very sorry to bother you again. > > For the following question, as you suggested, I posted it in both Biostars website and stackexchange website, so far no reply. > > I really hope that you can do me a great favor to share your points...

...uot;bottom", legend.title = "Tree Cluster", legend.labs = legend.labs, lty=1, lwd=3) } From: Bert Gunter [mailto:bgunter.4567 at gmail.com] Sent: Sunday, January 14, 2018 9:50 AM To: Ding, Yuan Chun <ycding at coh.org> Subject: Re: [R] consolidate three function into one I have not looked at your code *at all*, so I am not sure what exactly you mean by "consolidate." Is this what you are after? f1 <- function(...) { code 1} f2 <- function(...) { code 2} f3 <- function(...) { code 3} #...

...must print a ggplot object, for example with print(m52.2cluster) For the FAQ, run the line system.file("../../doc/FAQ") in R on your computer. Open up the resulting filepath in your favorite editor and scroll down to 7.22 On Sun, Jan 14, 2018 at 4:21 PM, Ding, Yuan Chun <ycding at coh.org> wrote: > Hi Bert, > > I am sorry to bother you on weekend. > > I am still struggling on defining a correct function. > > I first defined the function RFS (see below), then run it by provide the two argument. > > m52.2cluster <-RFS(inputfile =allinfo_m52, N=2 )...

...statistical site. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Fri, Mar 2, 2018 at 10:34 AM, Ding, Yuan Chun <ycding at coh.org> wrote: > Dear R users, > > I need to analyze data generated from a partial two-by-two factorial > design: two levels for drug A (yes, no), two levels for drug B (yes, no); > however, data points are available only for three groups, no drugA/no > drugB, yes drugA/no drugB,...