similar to: Date and selection

Dear R experts, I have plot the PDF of t distribution with df = 74. curve(dt(x,df=74),from=-4, to=4) how can I shade the area under curve (for example, col="red") from t=+- 1.996? Thank you. Regards, CH -- CH Chan Research Assistant - KWH http://www.macgrass.com

Dear all, I coded a function called u.glm u.glm <- function (x,ahi,age,bmiz,gender) { library(nortest) lil.rslt <- lillie.test(x) if (lil.rslt$p.value >0.05) { cat("Logtrans=0, lillie=",lil.rslt$p.value,"\n") xmodel<-glm(x~ahi+age+bmiz+as.factor(gender)) summary(xmodel) confint(xmodel) } else { cat("Logtrans=1,

Dear R helpers, I have a question regarding wordpress and R. I have asked this question in Wordpress support ( http://wordpress.org/support/topic/368312 ) but there is no answer so far, maybe the R community can have a better answer. I followed the excellent tutorial on r-statistics-blog ( http://www.r-statistics.com/2010/02/r-syntax-highlight-on-your-blog-a-wordpress-plugin/ ) to enable syntax

Dear R programmers, I can only config. proxy under Mac OS X terminal and launch R under Terminal by Terminal: export http_proxy="http://un:pw at proxy.com:port" Under R: >chooseCRANmirror(graphics=FALSE) >update.packages() I don't know how to config this in R for Mac OS X Aqua GUI. I checked the relevant document but no specific steps for Mac OS X. Regards, Ch Chan --

Hello, is there an up-to-date reference for how many people use R? I'm giving an R demo and want to cite wonderful R usage stats. How many people use it (or download it)? How often is R used in peer-reviewed pubs, etc. Is there any whiz-bang citation that says something like "R is great and developed by the best minds in statistical computing." Any thoughts? -Andy

Hi everybody, somehow i dont get the shapiro wilk test for normality. i just can?t find what the H0 is . i tried : shapiro.test(rnorm(5000)) Shapiro-Wilk normality test data: rnorm(5000) W = 0.9997, p-value = 0.6205 If normality is the H0, the test says it?s probably not normal, doesn ?t it ? 5000 is the biggest n allowed by the test... are there any other test ? ( i know qqnorm

Hi all, I am trying to figure out the formula used by R's Spearman rho (using cor(method="spearman")) because I can't seem to get the same value as by calculating "by hand". Perhaps I'm using "cor" wrong, but I don't know where. Basically, I am running these commands: > y=read.table(file="tmp",header=TRUE,sep="\t") >

Dear list, I need to read a big txt file (around 130Mb; 23800 rows and 49 columns) for downstream clustering analysis. I first used "Tumor <- read.table("Tumor.txt",header = TRUE,sep = "\t")" but it took a long time and failed. However, it had no problem if I just put data of 3 columns. Is there any way which can load this big file? Thanks for any suggestions!

Dear colleagues, I'm conducting a meta-analysis of studies evaluating adherence of HIV-positive drug users into AIDS treatment, therefore I'm looking for some advice and syntax suggestion for running the meta-regression using proportions, not the usual OR/RR frequently used on RCT studies. Have already searched already several handbooks, R-manuals, mailing lists, professors, but... not

Dear R Gurus, I have a problem related to plot. For example, I have two variables, pre and post. pre <- c(1,2,3,4,5) post <- c(2,5,7,2,3) How can I plot a line graph similar to this one? http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1847566&rendertype=figure&id=F1 Would you please provide me a keyword so that I can search. I don't know what is the name of this

In a discussion on IRC, Micheal was complaining about the semantics of StringSet, specifically how the iterators work. DannyB mentioned that we have a Trie class, so I decided to check it out. Lo and behold, I couldn't find _any_ uses of it, and when I deleted it from my system and did a clean rebuild, everything was fine. Looking at the log, the last substantial change to this class was in

Hi all, I need to put two different axis to one plot. On the top of the plot, I need to put one axis, with increments from left side to the right side; then at the bottom of the same plot, I need to put another axis, with increments from right side to the left side and showing a different unit. How do I do that? By the way, is there a "hold" command for plotting? If I first plot a

i am running asterisk-0.4.0 (Micheal asked with version am running)

Delete it. =] Thanks. On Tue, Sep 4, 2012 at 9:53 PM, Marshall Clow <mclow.lists at gmail.com> wrote: > In a discussion on IRC, Micheal was complaining about the semantics of > StringSet, specifically how the iterators work. > > DannyB mentioned that we have a Trie class, so I decided to check it out. > Lo and behold, I couldn't find _any_ uses of it, and when I deleted

Dear R-Helpers: I'm having a problem getting dates into the correct format. I have a data frame, which is based on a .csv file that I imported into R via read.table. R has converted my date variables to factors; when I use the as.Date command, most of the values are converted "correctly" (and by this I guess I mean converted "as I wish them to be") but some have not been.

Dear R help, Apologies for the less than informative subject line. I will do my best to describe my problem. Consider the following matrix: mdat <- matrix(c(1,0,1,1,1,0), nrow = 2, ncol=3, byrow=TRUE, dimnames = list(c("T1", "T2"), c("sp.1", "sp.2", "sp.3"))) mdat In my actual data I have time

Some what In-Reply-To: <1323092879.2629413.1445740496536.JavaMail.yahoo at mail.yahoo.com> On Sun, Oct 25, 2015 at 02:34:56AM +0000, Patrick Masotta via Syslinux wrote: > Micheal Brown > > > > iPXE's HTTP implementation is definitely not slow. On a Gigabit LAN, > > you should get the full 1000Mbps speed for HTTP downloads: a 200MB test > > file should

Hello Micheal, this reminds me of something, that I experienced in the past. Why would the server! complain "Unknown CA"? To test inspect the communication with wireshark and look if the client sends a cert; or: $ echo "a001 LOGOUT" | openssl s_client -msg -connect your.server:993 and grep for "CertificateRequest". Do you have a certificate configured in your

I'm having some trouble figuring out latest.start in the tframe package. The results seem a bit cloudy. Sample session: > library(ts) > library(tframe) Loading required package: syskern > t1<-ts(c(1,2,3,4,5),start=1) > t2<-ts(c(1,2,3,4,5),start=2) > earliest.start(tbind(t1,t2)) [1] 1 1 > earliest.start(tbind(t2,t1)) [1] 1 1 > latest.start(tbind(t2,t1)) [1] 1 1

Dear R Users, I would like to know is there any package to create a plot like this? http://dl.dropbox.com/u/5409929/cs1160521f01.gif X axis is categorical. And the positions of the points are corresponding to the frequency. (similar to violinplot) Thank you. Regards, CH -- CH Chan