similar to: Axis trouble

Dear All, I have built a package which has a data.frame "annotIndex.rda in its "data" folder. I am using this data frame within two functions in my package. Though my package works fine, yet when I do a R CMD check to my package I get the following two notes get.affy.ensembl: no visible binding

Hello, I want to plot 6 line graphs. I have 10 points 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0. At each point say 0.1, I have 6 variables A, B, C, D, E and F. The variables all have values between 0 and 1 (and including 0 and 1). I also want to label the x axis from 0.1 to 1.0 and the y axis from 0.1 to 1.0. My goal is to plot a line graph representing the mean of the variables at each

Dear All, I am trying to build a package for a set of functions. I am able to build the package and its working fine. When I check it with R CMD check I get a following warning : no visible global function definition for ‘biocLite’ I have used biocLite to load a user defined library from within a function if that library is not pre-installed

Hi all. I'm working on an assignment for a psychology class and I am not sure how to adjust the y-axis so it displays the range: 0, 5, 10, 15 The code below is almost ideal: lsd = c(3, 5, 13) mar = c(1, 2, 3) g_range <- range(0, lsd, mar) plot(lsd, type="o", col="blue", ylim=g_range, axes=FALSE, ann=FALSE) axis(1, at=1:3,

Dear All, I am running R in a system with the following configuration *Processor: Intel(R) Xeon(R) CPU X5650 @ 2.67GHz OS: Ubuntu X86_64 10.10 RAM: 24 GB* The R session info is * R version 2.14.1 (2011-12-22) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8

Hi all, I'm new to R and have the following problem: I have a 2 factor design (a has 2 levels, b has 3 levels). I have an object kidney.aov which is an aov(y ~ a*b), and when I ask for model.tables(kidney.avo, se=T) I get the following message along with the table of effects: Design is unbalanced - use se.contrast() for se's but the design is NOT unbalanced... each fator level

Dear All, I am executing my R script from PHP using the shell_exec. $return = shell_exec("/usr/bin/R --slave --args ".$Domain.",".$Gender." < /var/www/html/trends/newTrend.R"); In newTrend.R, I am plotting graph and save in pdf file using the below code. It will create a pdf file with two line. pdf(file="/var/www/html/trends/bmtrend.pdf", height=7.5,

Hi, I am trying to convert the following NLMIXED code to NLME, but am running into problems concerning 'Singularity in backsolve'. As I am new to R/S-Plus, I thought I may be missing something in the NLME code. NLMIXED *********** proc nlmixed data=kidney.kidney; parms delta=0.03 gamma=1.1 b1=-0.003 b2=-1.2 b3=0.09 b4=0.35 b5=-1.43 varu=0.5; eta=b1*age+b2*sex+b3*gn+b4*an+b5*pkn+u;

hey, I can not find a function for the following problem, hopefully you can help me. I have a vactor like this one v = c(NA,NA,TRUE,TRUE,NA,TRUE,NA,TRUE,TRUE,TRUE) and I would like to the TRUE values by the their "local sequence number". This means, the result should look thike this: c(NA,NA,1,2,NA,1,NA,1,2,3) Of course I could solve the problems using a loop, but this would be

I think the code is part of the RTAQ package but is not included in it, as I obtained it from https://r-forge.r-project.org/scm/viewvc.php/pkg/RTAQ/R/HAR_model.R?view=markup&root=blotter&sortby=author&pathrev=1028. It is not my code and I make no claim to other's good work, and apologize if I should even be posting it I am not sure, but in the transform function it allows to

hi, everyone: i have a question on the reshape function. i have the following dataset : gene tissue patient1 patient2 patient3............. _________________________________________________ gene1 breast 10 20 50 gene2 breast 20 40 60 gene3 breast 100 200 300 which i hope to convert to the following format: gene patientID

hello, I have got some trouble with R functions nlm(), nls() or optim() : I would like to fit 3 parameters which must stay in a precise interval. For exemple with nlm() : fn<-function(p) sum((dN-estdata(p[1],p[2],p[3]))^2) out<-nlm(fn, p=c(4, 17, 5), hessian=TRUE,print.level=2) with estdata() a function which returns value to fit with dN (observed data vactor) My problem is that only

Anybody who have a working .emacs (++) for Ubuntu (or Debian) with the newest emacs-rails (v.0.3). I have especially trouble with getting HTML Mode Snippets to work. - Henrik

Hi , I am planning to use an open source .ogg player on my ARM simulator on QEMU ( for ARM Application Baseboard ) . So far I have been able to bring up a Linux kernel 2.6.26.5 on the simulator . Right not I want an open source ogg vorbis player which I can use to call the libraries and play some files . However the catch is that I need the APIs exposed of the player application to me . I am

Dear R users - I have a list of patient identifiers and diagnoses from inpatient admissions. I would like to reorganize the list, presently in a long format to a wide format in reshape, but in the absence of a "time" element, I am uncertain how to do this - any help greatly appreciated. ID Dx A nausea A diabetes A kidney failure A heart attack A fever B fever B

Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ######################################################> kfitm1 Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = "gauss"), data = kidney)

Hi, I need help to perform a Shapiro.test on a data frame, I know that this test works only with vector but I guess there most be a way to permor it on a data frame instead of vactor by vector (i.e. I've got 40 variables to analyze and its kinda annoying to do it one by one) Thanks to anyone that can help me. Gonzalo Quiroga

hey guys, I am trying to add mapping functionality to the example application from the book. The mapping (mapquest''s openAPI) uses javascript to generate maps. So I was thinking that I''ll have a partial view to which i''ll pass the addresses and it''ll plot them. I''ll then put the partial view in the main view. Is this a good way of implementing this

Hi All, I have a data frame with a group of variables named b1, b2, b3, ..., b18. These variables take the value 1, 2 or NA. For each observation, I want to do some computation by looping over the values for the group of variables: b1 to b18. In STATA I would do: forval i=1/18 { --- use b`i' for computation ---- } How can this be done in R? Deepankar

Hi Dear R-users I have a database with 18000 observations and 20 variables. I am running cox regression on five variables and trying to use survfit to plot the survival based on a specific variable without success. Lets say I have the following coxph: >library(survival) >fit <- coxph(Surv(futime, fustat) ~ age + rx, data = ovarian) >fit what I am trying to do is plot a survival