similar to: Broken x-axis?

I am trying to read and write database tables that have a date field in them. I am constructing a data.frame, then using dbWriteTable to create the table and dbReadTable to read it. >datatbl<-data.frame( dates, trialnum, fooddel, ethdel, trialtime, trialtype, deliveries, food, ethanol, fcumrec, dcumrec, rrf, rrd) >ifelse(startdate==filelist,

Hi! I use simtest fonction of multcomp package to compile a Dunnett's test. I have 10 treatments and one control group, so i create a matrix with: m<-matrix(0,10,11) m[1,1]<--1 m[1,2]<-1 m[2,1]<--1 m[2,3]<-1 m[3,1]<--1 m[3,4]<-1 m[4,1]<--1 m[4,5]<-1 m[5,1]<--1 m[5,6]<-1 m[6,1]<--1 m[6,7]<-1 m[7,1]<--1 m[7,8]<-1 m[8,1]<--1 m[8,9]<-1

Dear list, I tried to use the nparcomp to run some post hoc non-parametric comparison and got and error. Error in uniroot(pfct, interval = interval) : f() values at end points not of opposite sign Appreciate any comments. the command line: >nparcomp(Ulceration~Group,data=test,type='Dunnett',control='Non-treated') Jun

Dear R-experts, I need to replace the values of a vector(tx) with a word ('Vehicle') when the value of the vector contains the word 'vehicle'. Sometimes, the value could be 'MCT vehicle', or 'control-vehicle', etc. I tried gsub like this treatment<-gsub('vehicle','Vehicle', tx, ignore.case=T) But then I end up with values like 'MCT

I am working on creating an R package for doing fire department analysis and am trying to create a function that can display emergency incident densities. The following code sort of does the trick, but I need a display that shows the number of incidents per square mile. I believe the code below shows incidents per square unit (in this case, degrees lat/long). To solve this problem, I believe

Hello, I would like to use serialized symbols as values for select in forms. eg. class Vehicle < ActiveRecord::Base @@statuses = {:wait => ''No yet'', :ok => ''Run run run''} validates_inclusion_of status, :in => @@statutes.keys Serialize status, Symbol def self.statuses @@statuses end end Then, I can use Vehicle.statuses.invert to

Hello all! I''m having an issue with Observers, and I''m hoping someone knows the answer. Here is the definition: app/models/vehicle_observer.rb > class VehicleObserver < ActiveRecord::Observer > def after_save(vehicle) > breakpoint > if vehicle.has_new_mil_indication? > UserMailer.deliver_mil_notice(vehicle) > end > end > end

Hello there, I''m developing an application which has models like this: class Vehicle < ActiveRecord::Base belongs_to :internal_dim, :class_name => "Dimension" def before_create() idim = Dimension.create self.internal_dim_id = idim.id #self.create_internal_dim # Doesn''t work end end class Dimension < ActiveRecord::Base

Hello, Hierarchical factors are a very common data structure. For instance, one might have municipalities within states within countries within continents. Other examples include occupational codes, biological species, software types (R within statistical software within analytical software), etc. Such data structures commonly use hierarchical coding systems. For example, the 2007 North

I updated my source tree and after putting malloc.h back things look a lot better. Still getting a large number of JIT-only failures. There is now only one assertion failure: /home/cfe/x86/llvm-gcc/bin/g++ -I/usr/home/llvm/obj/projects/llvm-test/MultiSource/Benchmarks/Prolangs-C++/city -I/usr/home/llvm/obj/projects/llvm-test/../../../projects/llvm-test/MultiSource/Benchmarks/Prolangs-C++/city

OEM software means no CD/DVD, no packing case, no booklets and no overhead cost! So OEM is synonym for lowest price. Buy directly from the manufacturer, pay for software ONLY and save 75-90%! Check discounts and special offers! Find software for home and office! TOP ITEMS Adobe Photoshop CS2 V9.0 $69 Macromedia Flash Prof 8 $49 Adobe Premiere 2.0 $59 MS Office

OEM software means no CD/DVD, no packing case, no booklets and no overhead cost! So OEM is synonym for lowest price. Buy directly from the manufacturer, pay for software ONLY and save 75-90%! Check discounts and special offers! Find software for home and office! TOP ITEMS Adobe Photoshop CS2 V9.0 $69 Macromedia Flash Prof 8 $49 Adobe Premiere 2.0 $59 MS Office

Let's say I've run Anova(lm(y~a*b)) and found the a:b interaction to be significant. Now I'm interested in which specific level combinations of a and b significantly differ from the control group. Can I use the t-tests from summary(lm(y~a*b)) to answer that question? I saw no mention of multiple comparison in the documentation for summary.lm, so am I right in assuming I need to

Hi, I used subset command, like this: grupoP = subset(dados, grupos=='P', select=c(mortos, vivos, doses, percevejos, p)) and the variables in select option are numeric. I tried fit a model with command: ajuste.logit = glm(cbind(mortos,percevejos)~log10(doses), family=binomial(logit), data=grupoP) and the output is: > ajuste.logit = glm(cbind(mortos,percevejos)~log10(doses),

On Wed, 17 Nov 2004, Jeff Cohen wrote: > I updated my source tree and after putting malloc.h back things look a lot better. Still getting a large number of JIT-only failures. There is now only one assertion failure: > > /home/cfe/x86/llvm-gcc/bin/g++ -I/usr/home/llvm/obj/projects/llvm-test/MultiSource/Benchmarks/Prolangs-C++/city

I am trying to reproduce some old SAS PROC MIXED code using R and nlme. The data consists of emission readings from vehicles and fuel properties. All variables are real numbers except "study" and "vehicle", which are character. Unfortunately, since the data are confidential, I'm unable to provide an example. The original SAS v6.12 code is provided below:

Hello, I have an application that access two different DBMS - Firebird and Postgres. So far, that''s ok. The problem is, depending on the customer logged, a different database (Firebird) must be selected. To make the things clearer: Postgres: 1 database shared by all customers with these tables: - user - vehicle - driver Firebird: several databases with different names and identical

Dear Dennis, Your question assumes that people know both SAS PROC MIXED and R nlme. Only a limited number of people do. Add the mathematical formulation of the model. That will increase the number of people that can help you. Adding the number of levels in each categorical variable and the number of observation per group is useful too. Best regards, ir. Thierry Onkelinx Instituut voor natuur-

Hi all, I asked this before the holiday, didn't get any response. So would like to resend the message, hope to get any fresh attention. Since this is not purely lme technical question, so I also cc-ed R general mailing list, hope to get some suggestions from there as well. I asked some questions on how to analyze longitudinal study with only 2 time points (baseline and a follow-up)

Hi, all I am trying to use the label_wrap_gen function in this website. https://github.com/hadley/ggplot2/wiki/labeller I tried to make a long name like this Light and heavy good vehicles (diesel) -\nGVX f2 = facet_grid(vehicle ~ ., labeller=label_wrap_gen(width=15)) eventually, I got something like this in my label... *Light and heavy good vehicles (diesel) - GVX* I suppose the