similar to: How do I plot geoms in parallel in ggplot

Hello R-List I've run into a minor problem with ggplot2. In particular, I get an error message when trying to plot a horizonal line through zero. Here is what I do: plt.df <- subset(xfm.df, select=c(stage, subj, tX)) plt <- ggplot(pp.df, aes(x=stage, y=tX, group=subj)) + geom_line() plt plt <- plt + geom_hline(intercept=0) plt The first plot yields a most lovely graph.

All, How does one extract the level-1 variance from a model fit via lmer()? In the code below the level-2 variance component may be obtained via subscripting, but what about the level-1 variance, viz., the 3.215072 term? (actually this term squared) Didn't see anything in the archives on this. Cheers, David > fm <- lmer( dv ~ time.num*drug + (1 | Patient.new), data=dat.new )

All, I'm trying to plot an all.effects() object, as shown in the help for all.effects and also Crawley's R book (p.178, 2007). The data has a repeated measures structure, but I'm using all.effects for the simple lm() fit here. Below is a reproducible example that yields the error message. fm.ex = lm(dv ~ time.num*drug*X, data = dat.new) fm.effects = all.effects(fm.ex, xlevels =

Hi ! I am a new user and quite confused by R-indexing. Make a list and get the attributes lst <- list(x = 1:3, y = 4:6, z = 7:9) attributes(lst) This returns: $names [1] "x" "y" "z" I can easily do: nm <-names(lst) or nm <-attr(lst,"names") which both return the assigned names of the named list 'lst', but why then this doesn't

Hello On Wed, 2008-01-23 at 00:10 +0300, Nikita Danilov wrote: > Hello, > > below is a tentative list of tasks into which WBC effort can be > sub-divided. I also provided a less exact list for the EPOCH component, > and an incomplete list for the STL component. > > WBC tasks are estimated in lines-of-code with the total of (9100 + 3000) > LOC, where LOC is a non-comment,

hi guys i have CentOS-5.1 with samba-3.0.25b-0.el5.4 and i took this workin well smb.conf from slackware12 ( commands to be done befor this ): # useradd wbc -c /home/wbc -g users# smbpasswd -a wbc# mkdir /export# chmod u+rwx,g+rx,o+rx /export# chown wbc.users /export the conf file:[global]workgroup = WORKGROUPnetbios name = sanshirosecurity = share [share]comment = Datapath = /exportforce user

Hello, First, ggplot2 is great! Second, sorry for the basic nature of my question. I have data in the form of a melt table like this: Subject Day Ab Variable Value 1 1 Yes A 3 1 3 Yes A 5 1 5 Yes A 7 2 1 No A 2 2 3 No A 4 2 5 No A 6 3 1 Yes A 1 3 3 Yes A 3 3 5 Yes A 5 4 1 No A 4 4 3 No A 6 4 5 No A 8 1 1 Yes B 3 1 3 Yes B 5 1 5 Yes B 7 2 1 No B 2 2 3 No B 4 2 5 No B 6 3 1 Yes B 1 3 3 Yes B 3 3 5

I am using smbldap-tool from IDEALX on my Fedora Core 1 samba 3.0.5 PDC. I have been following their guide step by step and I am getting the following error message as I try to get SID. I have looked in google group, but I couldn't exactly find anybody who had the same problem as I do. Please help me if you can! [root@pdc smbldap-tools]# net getlocalsid [2004/06/17 10:53:57, 0]

Hello all, I am re-posting my previous question with a simpler, more transparent, commented code. I have been ramming my head against this problem, and I wondered if anyone could lend a hand. I want to make parallel a bootstrap of a linear mixed model on my 8-core mac. Below is the process that I want to make parallel (namely, the boot.out<-boot(dat.res,boot.fun, R = nboot) command).

I know that this has been asked before in other variations but I just can't seem to figure out my particular application from previous posts. My apologies if I have missed the answer to this question somewhere in the archives. I have indeed looked. I am running Ubuntu 11.04, with R 2.12.1 and ESS+Emacs. For journal formatting requirements, I need to italicize all the greek letters in any

Dear R user, I'm relative new and need some help / advise from you. I have organised my data in several frames, mainly because the data came from different sources. One common variable in all data frames is the day post transplantation. I would like to combine the data frame to one big data frame which contains all the data, however no duplicate rows should occour (only one row for one

hi, now i did some tests with the "auth methods" option. "auth methods = sam, sam_ignoredomain, unix, winbind, wbc, trustdomain, ntdomain, guest" -> RDP works -> smbclient -k does not work "auth methods = sam, sam_ignoredomain, unix, winbind, wbc, trustdomain, ntdomain, guest, samba4" -> RDP hangs an the welcome-scrren -> smbclient -k works fine

Hi, I have sort of a newbie question. I've seriously put a lot of effort into how to handle simple replicates in a repeated ANOVA design, but haven't had much luck. I really liked reading "Notes on the use of R for psychology experiments and questionnaires", by Jonathan Baron and Yuelin Li ( http://www.psych.upenn.edu/~baron/rpsych/rpsych.html ) but still didn't run across

Hi , I am new to R. I am trying to run a simple R script as shown below: aov.R ------ data1<-c(49,47,46,47,48,47,41,46,43,47,46,45,48,46,47,45,49,44,44,45,42,45,45,40 ,49,46,47,45,49,45,41,43,44,46,45,40,45,43,44,45,48,46,40,45,40,45,47,40) matrix(data1, ncol= 4, dimnames = list(paste("subj", 1:12), c("Shape1.Color1", "Shape2.Color1", "Shape1.Color2",

quantreg package is used. *fit1 results are* Call: rq(formula = op ~ inp1 + inp2 + inp3 + inp4 + inp5 + inp6 + inp7 + inp8 + inp9, tau = 0.15, data = wbc) Coefficients: (Intercept) inp1 inp2 inp3 inp4 inp5 -0.191528450 0.005276347 0.021414032 0.016034803 0.007510343 0.005276347 inp6 inp7 inp8 inp9 0.058708544

I am extremely puzzled by this behavior in R. I have a data frame called Trials in which I have results from an experiment. I am trying to do a subjects analysis, but getting weird results. Each row has 1 trial in it, which includes a column for the subject number I get the list of subject numbers like so: > Subj=unique(sort(Trials$Subj)) Then I loop over them. But I get strange results. As

Hi all, We would like to propose a new Constant type in LLVM for representing entries in the Procedure Linkage Table (PLT). The PLT is a data structure used for dispatching position-independent function calls to appropriate functions where the address of the function is not known statically. Right now, if a call is made to a function, it may be lowered to a direct call to the function itself or

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,

Hi, I would like to represent the treatment course of a subject with an arrow representing passage through time and various ticks, arrows, or bars representing treatments done to the patient along the arrow. I would like to generate this sort of schematic for each patient in a database taking information from the database such as time of treatment and type of treatment. Does anybody know how I

Thanks... I agree that the problem was explained in the documentation but I can't find a way to have axis() working even manipulating par("plt") or with graphics.reset = TRUE: - adding graphics.reset=TRUE does not allow axis() to be shown; - I see that par()$plt is involved but it is the not sufficient to explain why axis() works because if it is changed by hand, axes are not