similar to: limit digits in text()?

Dear experts, I just noticed that I get different results conducting type III anova using drop1 or the Anova command from the car package. I suppose I made a mistake and hope you can offer me some help. I have no idea where I got wrong and would be very grateful for explaination as R is new terrain for me. If I run the commands in line, they produce the same results. But if I run them in

Hi all, I have an easy data set. It has three columns: Subject, Condition, dprime. A small excerpt follows, in order to illustrate: Subject Condition dprime HY s 3.725846 CM s 2.877658 EH s 5 HY st 2.783553 CM st 2.633955 EH st 5 I want to visualize this. What I thought of was having dprime on the y-axis (scale 0-5), Subject on the x-axis, and then two lines plotted to show the dprime value for

This example, x<-1:10 y1<-x y2<-x+2 plot(x,y1,pch=21,ylim=range(c(y1,y2)),type="o",lty="solid",bg="white") points(x,y2,pch=22,type="o",lty="dashed",bg="white") legend(7,3,legend=c("y1","y2"),pch=c(21,22),lty=c("solid","dashed"),bg="white") works if I use pch = one of 21-25. But

I use AfterStep on Linux. I made an R.xpm file (attached) and put this icon in the "wharf". To invoke R I click on the R icon. I did not find an offical R.xpm, and so I am submitting this one for your collective consideration.... The idea behind the icon is that it looks like a scatterplot, which is "statistical". My main beefs with my R.xpm are that it doesn't have a

Hi, I'm trying to calculate the value of the variable, dp, below, in the argument to the integral of dnorm(x-dp) * pnorm(x)^(m-1). This corresponds to the estimate of the sensitivity of an observer in an m-alternative forced choice experiment, given the probability of a correct response, Pc, a Gaussian assumption for the noise and no bias. The function that I wrote below gives me an error:

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

I was wondering if R has an equivalent of #ifdef. I want the code to check if a function has been sourced in yet, and if not, then source it in. (I know of replace(), which works for packages/libraries) Here is what I have come up with. my.objects<-objects() if(my.objects[my.objects=="thingy2"]!="thingy2") source("~/papers/speed/junk.r") junk.r contains the

Guys, I wrote to the finance mailing list earlier with my questions but was directed here. Sorry for the repeat. --------------- library(quantmod) .... now <- Sys.time() midnight <- strptime() # <---- I want to make this a static variable that will be equal to 12:00:00 am but I dont know what to put here. I keep getting NA for everything I do if(now == midnight) {

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

I do: x<-rnorm(1000) quantile(x,c(.025,.975)) 2% 98% -1.844753 1.931762 Since I want to find a 95% confidence interval, I take the .025 and .975 quantiles. HOWEVER R says I have the 2% (not 2.5%) and 98% (not 97.5%) points. Is it just rounding the printed 2% and 98%, or is it REALLY finding .02 and .98 points instead of .025 and .975? Thanks for any help. Bill Simpson

I do: x<-rnorm(1000) quantile(x,c(.025,.975)) 2% 98% -1.844753 1.931762 Since I want to find a 95% confidence interval, I take the .025 and .975 quantiles. HOWEVER R says I have the 2% (not 2.5%) and 98% (not 97.5%) points. Is it just rounding the printed 2% and 98%, or is it REALLY finding .02 and .98 points instead of .025 and .975? Thanks for any help. Bill Simpson

> On 8 Sep 2017, at 15:51 , Martin M?ller Skarbiniks Pedersen <traxplayer at gmail.com> wrote: > > On 8 September 2017 at 14:37, peter dalgaard <pdalgd at gmail.com> wrote: >> >> >>> On 8 Sep 2017, at 14:03 , peter dalgaard <pdalgd at gmail.com> wrote: >>> >>> x <- scan("~/Downloads/digits.txt") >>> x <-

On 8 September 2017 at 14:37, peter dalgaard <pdalgd at gmail.com> wrote: > > > > On 8 Sep 2017, at 14:03 , peter dalgaard <pdalgd at gmail.com> wrote: > > > > x <- scan("~/Downloads/digits.txt") > > x <- x[-seq(1,220000,11)] > > ...and, come to think of it, if you really want the 1000000 random digits: > > xx <-

> On 8 Sep 2017, at 14:03 , peter dalgaard <pdalgd at gmail.com> wrote: > > x <- scan("~/Downloads/digits.txt") > x <- x[-seq(1,220000,11)] ...and, come to think of it, if you really want the 1000000 random digits: xx <- c(outer(x,10^(0:4), "%/%")) %% 10 -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg

Simplest version that I can think of: x <- scan("~/Downloads/digits.txt") x <- x[-seq(1,220000,11)] length(x) # 200000 hist(x) Now, because it's Friday: How does one work out the theoretical distribution of the following table? > table(table(factor(x,levels=0:99999))) 0 1 2 3 4 5 6 7 8 9 10 11 13497 27113 27010 18116 9122

Howdy! I have a problem that's similar to one discussed earlier on the list: When not in xterm all the backspace and cursor-movements and stuff are working well. Only when I start R in an xterm window backspace appears as ^H and so on. What can I do? Ragnar -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I've gotten my CDR working the way I like, but I am looking to customize it a bit. I have set up an IVR menu, which works great. I would like to be able to capture the prompted DTMF digits pressed by callers, to my CDR database but I don't see any AGI or Asterisk commands that allow one to customize the CDR contents. Am I thinking about this on the wrong track? If someone calls sales

================================================= This is an automated summary of the status of the R-bugs repository. Note that this may be neither complete nor perfectly correct at any given instance: Not all bugs are reported, and some reported bugs may have been fixed, but the repository not yet updated. Some bug fixes are difficult to verify because they pertain to specific hardware or

I did some serach around. It seems that ROC curve computation is not supported on R. Anyone has some leads? Thanks. Jason ===== Jason G. Liao, Ph.D. Division of Biometrics UMDNJ School of Public Health 335 George Street, Suite 2200 New Brunswick, NJ 08903-2688 phone (732) 235-9748, fax (732) 235-9777 http://www.geocities.com/jg_liao __________________________________________________ Great