similar to: modelling problem

I fit the model fit<-lm(thresh~cond*Ne) where thresh is the reponse cond is a factor with levels a, b, and c Ne is a continuous indep var I think of this full model as having three lines: thresh as a function of Ne for each condition. Thus we have slopea, slopeb, slopec, inta, intb, intc. lm output my params ------------------------- (Intercept) inta condb intb - inta condc

I wanted to fit a nonlinear model using nls to two groups. Here is a toy example. library(nlme) x<-rep(1:5,2) cond<-c(rep("a",5),rep("b",5)) cond<-as.factor(cond) y[cond=="a"]<-1/x +2 y[cond=="b"]<-2/x+3 df<-data.frame(x,y,cond) fit<-nls(y~c(a1,b1)[cond]/x+c(a2,b2)[cond],data=df, start=list(a1=1,b1=2,a2=2,b2=3)) I based above on an

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

Full_Name: Bill Simpson Version: 1.2.3 OS: linux Submission from: (NULL) (193.62.250.209) I think I have found an accuracy problem with sum() > a<-1:1000 > b<-1:1000 > > crossprod(a,a)*crossprod(b,b) [,1] [1,] 1.114448e+17 > drop(.Last.value) [1] 1.114448e+17 > #this should give same answer as above > sum(a*a)*sum(b*b) [1] -652010736 #seems to be due to

There were some questions about hist() a couple of days ago which triggered this post. My question/suggestion is about the y-axis in hist. There are reasons to prefer making the y-axis density=relative frequency/bin width. One reason is that the height of the plot does not depend on the bin width; another is that if your histogram is in density then you can easily superimpose a smooth theoretical

I'll be even more tangent. Those interested in Ox, see http://www.de.ufpe.br/~cribari/ox.pdf Cheers, Francisco. Date: Fri, 6 Apr 2001 09:34:19 +0100 (BST) From: Bill Simpson <wsi at gcal.ac.uk> Subject: Re: [R] Using Gauss with R This is a tangent to your question. The economist Jurgen Doornik has written a language called Ox: http://www.nuff.ox.ac.uk/Users/Doornik/doc/ox/ox.htm

I was wondering if anyone out there has created a series of images (e.g. jpegs) using R, then animated them on a web page. I have looked around a bit and it seems fairly complicated. One approach that seems sensible is to use a java script that displays the list of images. I have only used html on web pages before so I don't know how to use java... Ideally the final result is easily ported to

Hallo! I have run a EEG experiment and got the following data group: 1 2 1 2 1 2 1 2 ... as factor, 2 levels between subjects fixed effect (patient vs control) subj: 1 2 ... 14 1 2 ... 14 as factor 7 patients 7 control random effect condition: 1 1 ... 2 2 ... 1 1 ... 2 2 as factor, 2 levels within subjects, ie every subject worked on every cond fixed effect roi: 1 ... 2 ... 3 ... 4 ... as factor,

Full_Name: Bill Simpson Version: 65.1 OS: linux Submission from: (NULL) (193.62.250.209) Try this code as separate lines entered interactively, then try doit() doit<-function() { x<-seq(1,10) y1<-x+rnorm(10,1,1.5) y2<-x+rnorm(10,1,1) t.test(x,y1,paired=TRUE) t.test(x,y2,paired=TRUE) } doit() apparently executes only the last of a series of t.test()s. Maybe this is no bug,

I have not really used lm before and I was hoping for some help on a simple problem. Here is a toy version of the problem I want to solve. y x grp -.9 1 a -.8 2 a -.7 3 a -.7 1.5 b -.5 2.5 b -.3 3.5 b -.19 2.7 c -.11 3.7 c -.41 4.7 c I want to fit a model that has one y-intercept and three slopes, one for

Full_Name: Bill Simpson Version: 0.64.1 OS: linux Submission from: (NULL) (193.62.250.209) Two problems with dev.print help. 1. Under Examples: dev.print(width=6, height=6, horizontal=F) ^ FALSE 2. dev.print help does not describe any of its parameters (e.g. height or width). Not described on the other dev.xxx pages either

If there are no objections, I'll submit a patch for review that sets the default value of "-enable-cond-stores-vec" to "true". Thanks! -- Matt On Wed, Dec 14, 2016 at 12:55 PM, Michael Kuperstein via llvm-dev < llvm-dev at lists.llvm.org> wrote: > I haven't verified what Matt described is what actually happens, but > assuming it is - that is a known

I just installed 99.0a, and I can't set mar and mgp: > par(mar(2,2,.1,.1)) Error in par(mar(2, 2, .1, .1)) : couldn't find function "mar" I had the same happen with mgp. Are these bugs or a problem in my installation (everything else seems OK)? Is R core going to kill mar and mgp? If so, how should I fix the huge margins and huge space between axis title and axis? Try this

Full_Name: Bill Simpson Version: 65.1 OS: Linux Redhat 6.1 Submission from: (NULL) (193.62.250.209) Try this: #default plot symbols and lettering are too small, need to scale up par(cex=2,mex=.7) par(mar=c(5,5,1,1)) spdiff<-c(1,2,3,4,5) dpdet<-c(1,2,3,4,5)/10 dpsp<-c(2,3,4,5,6)/10 dpdir<-c(2,4,6,8,10)/10 plot(spdiff,dpdet,pch=15,ylim=c(0,max(dpdet,dpsp,dpdir)),xlab="Speed

The problem is illustrated below: line<-function(x,y) { fit<-lm(y~x) fit } func<-function() { x<-c(1,2,3,4) y<-c(1.1,2,3,4) for(i in 1:2) { fit<-line(x,y) summary(fit) } } Why don't I get two printed summaries when I run func()? How to make it work? (In reality I will be getting different summaries of different data sets using this loop) Thanks very much for any

Thanks Michael and Dibyendu for doing the experimentation and bringing this up to our attention. It might be the case what Matt described here. I will take a look at it. Farhana From: Michael Kuperstein [mailto:mkuper at google.com] Sent: Wednesday, December 14, 2016 9:56 AM To: Das, Dibyendu <Dibyendu.Das at amd.com>; Aleen, Farhana A <farhana.a.aleen at intel.com> Cc: Matthew

I want to display a series of plots: display plot 1 wait for user input (keypress or mouse click) display plot 2 etc How to do it? Thanks for any help. (I have been using locator().) Bill Simpson -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi Matt- Yeah I used a pretty recent llvm (post 3.9) on an x86-64 ( both AMD and Intel ). -dibyendu From: Matthew Simpson [mailto:mssimpso at codeaurora.org] Sent: Wednesday, December 14, 2016 10:03 PM To: Das, Dibyendu <Dibyendu.Das at amd.com> Cc: Michael Kuperstein <mkuper at google.com>; llvm-dev at lists.llvm.org Subject: Re: [llvm-dev] Enabling scalarized conditional stores in

I have a data file like this: std cf hit miss fa cr 920 980 40 15 14 31 950 1010 24 23 23 30 1190 1250 26 21 27 26 1010 1070 33 10 28 29 1040 1100 35 10 11 44 I use read.table to read it in. My problem is that read.table makes std and cf into factors. I want them just

Hi Dibyendu, Are you using a recent compiler? What architecture are you targeting? The target will determine whether the vectorizer thinks vectorization is profitable without having to manually force the vector width. For example, top-of-trunk vectorizes your snippet with "clang -O2 -mllvm -enable-cond-stores-vec" and "--target=aarch64-unknown-linux-gnu". However, with