similar to: saving a regression model to a file

Hello all, I have generated a principal components regression model using the pcr() function from the PLS package (R version 2.12.0). I am getting a "non-conformable arguments" error when I try to use the predict() function on new data, but only when I try to read in the new data from a separate file. More specifically, when my data looks like this #########training data

is there a way to use the actual index value for plotting zoo objects this is the way that the index is set up and a sample range of what I would like to plot 01/01/06 00:00:00 - 01/01/06 23:45:00 { library(zoo) # chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) }} x <- structure(c(15.57, 15.5,

Dear Debian persons I would like enquire about whether the repository deb https://cloud.r-project.org/bin/linux/ubuntu bionic-cran35/ I tried to follow the instructions and I cannot get R3.5.1 - I get R3.4.4??? I have been able to install R3.5.1 on Ubuntu 16.02 using the xenial-cran35/ repository & so not sure if the bionic link is up to date? thanks in advance Regards Alison

Hello all, I am trying to use the micEcon 'insertRow' function inside a function I have written. For example: insert_row_test <- function(m){ insertRow(m,nrow(m)+1,v=0,rName="test") } However, when I try to call the 'insert_row_test' function (after loading the micEcon package), it does not insert a row into the matrix I pass in. When I call the insertRow

Hi all, I have 6 datasets(dataframes Assem_ContigsLen7 through all_ContigsLen12) containing 3 columns (contig_id, contig_length, read_count). Each dataset is composed of 3 types of contigs (assemblies of genomic fragments), 1- all Bacterial fragments, 2 - all Viral fragments, 3 - mixed fragments. I identified the type of contig through a merge with another table with just contig_id and

Hello, I am not only interested in finding out which genes are the most highly up- or down-regulated (which I have done using the linear models and Bayesian statistics in Limma), but I also want to know which genes are consistently highly transcribed (ie. they have a high intensity in the channel of interest eg. Cy5 or Cy3 across the set of experiments). I might have missed a straight forward

Hi, I have a table with data, as below: dput(table): structure(list(Adrian.Cole = c(0L, 0L, 0L, 0L, 0L, 0L), Alison.Wong = c(0L, 0L, 0L, 0L, 0L, 0L), Andrei.Savu = c(0L, 0L, 0L, 0L, 0L, 0L), Bruno.Dumon = c(0L, 0L, 0L, 0L, 0L, 0L), Edward.J..Yoon = c(0L, 0L, 0L, 0L, 0L, 0L), Eugene.Koontz = c(0L, 0L, 0L, 0L, 0L, 0L), Jakob.Homan = c(0L, 0L, 0L, 0L, 0L, 0L), Kelvin.Kakugawa = c(0L, 0L,

Hi, I''m trying to setup my mongrel clusters to start on boot. I have followed the instructions under the "On Boot Initialization Setup" section at http://mongrel.rubyforge.org/docs/mongrel_cluster.html. When I try to start it up, I get the following error: me at mybox.com: sudo /etc/init.d/mongrel_cluster start Starting all mongrel_clusters... !!! Path to log file not

Hello all, I am trying to use cbind to add a named empty column to a Matrix: outputmatrix <- cbind(outputmatrix,kog_id = seq(0,0,0)) The problem I have is that kog_id is a variable that has a value e.g. "KOG1234", but I when I try to use this to name the added column, it is named literally "kog_id" instead of "KOG1234". How can I name a column by passing in a

Hi all, I'm having problems installing plotrix. I tried installing it through install.packages, and from the unix command line, but each time it seems to stall when it is installing the help indices. has anyone had this same problem, is this package still maintained ? any help? thanks > install.packages("plotrix") > > I also tried using the source package > R CMD

Hi all, I have a large table mapping thousands of COGs(groups of genes) to pathways. # Ex COG0001 patha pathb pathc COG0002 pathd pathe COG0003 pathe pathf pathg pathh ## I would like to combine this information into a big list such as below COG2PATHWAY<- list (COG0001 = c ("patha ","pathb ","pathc

Hi folks, With thanks to Alison Keenan (another Alison!) for the voice, Chris Bagnal for converting from 44k wav to sln and finally Terje Elde for debugging my HTML code, the British English files are now ready for download. They can be got from http://www.enicomms.com/cutglassivr/ Thanks and don't forget to practice safe IAX ;-} Mark -- Mark Phillips <g7ltt@g7ltt.com>

CVS HEAD from about 1 week ago. TDM30P and call through Nufone. I was talking and wanted to park the call and move to another phone to pick it up. I hit #701 instead of #700 though -- after a pause, I got a fast busy and the call was gone. When I called the person back, she said that Alison told HER that 701 was an invalid extension. I should have heard that though, not her. If I dial

Ok..this seems like it shouldn't be that hard to figure out but I'm confused anyway. How can I make sure that whatever case a file is created with is what's shown in NT explorer? Let me explain: If I create a file in all uppercase on my unix file system (which is also a samba share) it shows up as a mixed case file within NT explorer, with only first letter capitalized. (Normally

If I want to create a heat map where I have a graphic of a US State with county borders, how do I import and work with the graphic image? The reason I would like to use R is so I can then define color parameters and quickly produce heat maps that correspond to various data in dataframes.   Thanks. Alison [[alternative HTML version deleted]]

I'd like to increase the speed of the Asterisk sound files. Miss Alison talks a bit slow. I can use sox to increase the speed, but then the pitch changes and she starts to sound like a chipmunk. Any audio experts out there know how I can increase the speed a little bit, and change the pitch accordingly so that she sounds ... normal? Thanks Doug.

Dear all, I want to use string split to parse column names, however, I am having some errors that I don't understand. I see a problem when I try to rbind the output from strsplit. please let me know if I'm missing something obvious, thanks, alison here are my commands: >strsplit<-strsplit(as.character(Rumino_Reps_agreeWalign$geneid),"\\.") >

Hello, Using the dse package I have estimated a VAR model using estVARXls(). I can perform forecasts using forecast() with no problems, but when I try to use simulate() with the same model, I get the following error: Error in diag(Cov, p) : 'nrow' or 'ncol' cannot be specified when 'x' is a matrix Can anyone shed some light on the meaning of this error? How can I

Hi, I have a data set that I'd like to run logistic regressions on, using ddply to speed up the computation of many models with different combinations of variables. I would like to run regressions on every unique two-variable combination in a portion of my data set, but I can't quite figure out how to do using ddply. The data set looks like this, with "status" as

Sorry for the newbie question but I''m stumped: I have an app that requires sudo to run mongrel_rails. For my capistrano recipe is do it like this: sudo "mongrel_rails start -e production -p #{application_port} -d -c #{current_path}" Which works famously. Now I have tried several things for my restart task. Here they are: sudo "cd #{current_path} && sudo