similar to: Mathematical annotation axis in lattice

Dear list, I am trying to write a package for simulating meioses in R. We defined a class 'haplotype' which contains the basic units of our simulation: setClass("haplotype",representation(snp = "numeric",qtl = "list", hID = "numeric",phID0 = "numeric",phID1 = "numeric"),

Dear list, I am trying to develop a package. I used the function package.skeleton to make the directory tree of the package and then build and compiled the package as described (I hope). Now, I would like to increase the number of functions in the package without overwriting the existing package directory tree, so I put these new functions in the directory R of the tree. When I then build and

Hello list, I am making graphics for an article which I want to publish. The article is about several methods (to calculate breeding values of individuals) applied in several genetic scenarios (scen1 in the example) and using data from two sources (scen 2 in my example). I want to specify the ylim of my plot and have relation = 'free' for the yaxis but I would to avoid plotting the axis

Dear forum, I'm trying to install RMySQL but I'm having problems in loading it, here is the message that I'm getting in R install.packages("RMySQL", dependencies=TRUE) trying URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.5/RMySQL_0.6-0.zip' Content type 'application/zip' length 391364 bytes opened URL downloaded 382Kb package 'RMySQL'

Hello, I have to add "Age (bar(x)=14.3) as a title on a chart. I am unable to get this to working. I have tried bquote, substitute and expression, but they are only doing a part of the job. new<- c(14.3, 18.5, 18.1, 17.7, 18, 15.9, 19.6, 17.3, 17.8, 17.5, 15.4, 16.3, 15, 17.1, 17.1, 16.4, 15.2, 16.7, 16.7, 16.9, 14.5, 16.6, 15.8, 15.2, 16.2, 15.6, 15, 17.1, 16.7, 15.6, 15, 15.8, 16.8,

Dear all,   I was using older version of R (installed early). I install new version of R (R.2.9.0) but i could not find package "xlsReadWrite" to read Excel file. Is there any alternatives to read Excel file ? I hope it should be. I need help although it is minor querry. Thanks in advance.   Sincerely, Ram Kumar Basnet Wageningen University, The Netherlands   [[alternative HTML

Dear all, I am trying to make barplot from matrix with "beside=FALSE". I have 165 labels in my bar plot. When i plot the label horizontally below the barplot, then all labels are not appear in the plot due to lack of space. I used the following function: barplot(t(bar), main="Barplot based on LCMS", ylab="RF membership Probability",   

Good morning, I am trying to recode the values of a column in a data frame and, instead of using a long series of if/else statements I would like to use a switch (as it is, for instance, in the C programming language). I have a column of this type: AT BB B1 B1 CC CC T6 B1 CC BB TT AT AT AT TT BB ... and I need to change this coding (if it's AT then it should be A1, if it's BB then

Dear R users,   Does somebody know the way to import many files by a single command in R ? I have 50 files in a directory and now, i am importing the files repeatedly (one by one). If there is a way to import all files at a time, it makes much more easy and save times too. Thanks in advance.     Sincerely, Ram Kumar Basent Wageningen University, the Netherlands   [[alternative HTML

Dear R users,   I am not so used to this R software. I have to use the package " Rgraphviz" but found some problem in the installation process. I download this package and store in R library but i am not getting this package in R installation list. I made review in google search net and use the following command:   ### source("http://bioconductor.org/biocLite.R")

Hello, I am a Masters student and I am working on my thesis modelling smallholder farms using a program in R. I have modified the original code and I am having some issues with lists that I cannot figure out. Originally, I had list file defining lists such as: Param, Crop1, Crop1, Soil, etc. (ex. Param <- list() ). Their subsets were listed as Crop1$CContent for example and there was quite

Hello, I have a list with n numerical components of different length (3, 4 or 5 values in each component of the list); I need to export this as a text file where each component of the list will be a row and where missing values should fill in the blanks due to the different lengths of the components of the list. I think that as a first step I should convert my list to a data frame, but this is

Hi   I think this may be simple task but creating trouble to me. In my calculation, i have got large number of rows (5546) and column 182 but i am not able to see all the rows. May be somebody have idea how i can see all the rows ? I will be greatful if any body help me. Thanks in advance   Sincerely, Ram Kumar Basnet Graduate student Wageningen University Netherlands   [[alternative

Hi R users, I want to calculate the partial correlation (first or second order) and their corresponding significance test (p-value). I used > library(corpcor) and "cor2pcor" function to calculate partial correlation but could not calculate "p-value". If some one knows the R package or function to calculate corresponding "p-value", It will be great help. Thanks  

Dear all,   I did ANOVA for many response variables (Var1, Var2, ....Var75000), and i got the result of p-value like below. Now, I want to select those predictors, which have pvalue less than or equal to 0.05 for each response variable. For example, X1, X2, X3, X4, X5 and X6 in case of Var1, and similarly, X1, X2.......X5 in case of Var2, only X1 in case of Var3 and none of the predictors in case

Dear R users,   It may be very simple but it is being difficult for me. I have two vectors with some common string. And, i want to combine into a vector in such a way that it includes string from both vectors and make a unique.   For example:   x <- paste(rep("A",5),1:5,sep = ".") x [1] "A.1" "A.2" "A.3" "A.4" "A.5" y <-

Hi, I am trying to produce a plot with an xlabel that reads (x >= 5%) with the >= turned into the correct epression. I can do this up to the percentage sign by specifiing xlab=expression(x>=5). Whatever I do to include the % sign as well doesn't seem to work. xlab=bquote(x>=5.("%")) almost works but includes brackets. Anyonw know how to solve this one Dan PS I am

Good morning, I am merging two datasets and I would like to save the non-matching rows in a separate file. The problem is how to retrieve the non-matching rows in R. Example: DATASET A code nomi A1 Franco A2 Mario A3 Andrea A4 Sandro A5 Luca DATASET B code book A1 Guerra e Pace A1 Storia di Roma A2 La coscienza di Zeno A4

Hello, How can I read a xlsx file using xlsx package? Thanks Albert [[alternative HTML version deleted]]

Dear R users,   I am new user for elastic net. I am trying to use elasticnet library. I have marker data with 359 markers and 168 samples, and response is metabolites. I am trying to do regression between a metabolite and markers.  But i am getting the following error:   > en<-enet(marker,as.numeric(vio),lambda=0.5,normalize=FALSE,intercept=TRUE) Error in one %*% x : requires numeric