similar to: metafor - interpretion of QM in mixed-effects model with factor moderator

I have a need to call and pass arguments to nlme() from within another function. I use R version 1.8. I have found an apparent way to make this work, but I would appreciate some comments on whether this fix is really appropriate, or there is another way to do it that does not involve changing the source code. I don't have enough experience to start changing the sorurce code of a library

I apologize for posting on this question again, but unfortunately, I don't have and can't get access to MASS for at least three weeks. I have found some code on the web however which implements the prediction error algorithm in cv.glm. http://www.bioconductor.org/workshops/NGFN03/modelsel-exercise.pdf Now I've tried to adapt it to my purposes, but since I'm not deeply familiar

Hi, I'm trying to fit a multinomial logistic regression to my data which consists of 5 discrete variables (scales 1:10) and 1000 observations. I get the following error: Error in `row.names<-.data.frame`(`*tmp*`, value = c("NA.NA", "NA.NA", : duplicate 'row.names' are not allowed In addition: Warning message: non-unique value when setting

*I have the following function:* /plot_mi_time = function(mdata, miname) { mdata2 = mdata[row.names(hakat) == miname, ] print(mdata2) xcoords <- c(1,1,2,2,3,3,4,4,5,5,6,6) plot(c(xcoords), mdata2, xaxt="n", ylab="Expression", xlab="Time(h)", , main=miname) axis(1, at=xcoords,

Hi, I want to control the plots in the output of the xyplot(). It is easier to explain it through an example: #------------------------------------------------------------- library(lattice); # months months <- c("Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep",

Hi, I have the following function that i use to plot graphs. plot_mi_time = function(mdata, miname) { mdata2 = mdata[row.names(mir_test3) == miname, ] # print(mdata2) xcoords <- c(1,1,2,2,3,3) plot(xcoords, mdata2, xaxt="n", ylab="Expression", xlab="Time(h)", , main=miname) axis(1, at=xcoords,

Hi all, I need to replace missing values in a matrix by 10 % of the lowest available value in the matrix. I've got a function I've used earlier to replace negative values by the lowest value, in a data frame, but I'm not sure how to modify it... nonNeg = as.data.frame(apply(orig.df, 2, function(col) # Change negative values to a small value, close to zero { min.val =

This post is NOT a question, but an answer. For readers please disregard all earlier posts by myself about this question. I'm posting for two reasons. First to say thanks, especially to Dimitris, for suggesting the use of errorest in the ipred library. Second, so that the solution to this problem is in the archives in case it gets asked again. If one wants to run a k-fold cross-validation

Hi I was trying to import a dataset from access database into R My code looks like below >library(RODBC) >setwd("//STATESFPSOP1/docs/CRMPO/DEPT/Retail PD") >channel<-odbcConnectAccess("MortgagePD.mdb") >Mdata<-sqlFetch(channel, "MortgagePD") But I got error messages saying Error: cannot allocate vector of size 27.9 Mb In addition: Warning

Hi everybody I have a large number of Excel speadsheets that I want to merge into one R dataframe to process. I can read them one by one with read.xls but I really need a function to read a whole directory at once and merge the columns with the same name into a single dataframe. here is what I have thus far (on Mac OS X 10.5.2. for R 2.8.0): # Define the directory where the data is

Hi everyone, I have some data about a market research which I want to arrange in one plot for easy viewing, the data looks something like: Product Color StoreA StoreB StoreC StoreD Price ProdA R NA 4.33 2 4.33 35 G NA 4.33 2 4.33 35 B

Hi, I've been trying to encode a live audio input from the microphone on iOS device using opus. Uncompressed audio recording works fine with http://theamazingaudioengine.com/ Then, when I tried to do encoding, I'm stuck at figuring out why the buffer is empty: static int ec_write_byte(ec_enc *_this,unsigned _value){ if(_this->offs+_this->end_offs>=_this->storage)return

am trying to produce two dot plot figures in ggplot2. So far the first one (p) in the program below is working fine. However when I want to move to a faceted plot (p1) I seem to lose my ordering or, more likely, I'm just getting an ordering I am not expecting and I always have trouble understanding ordering in R! What I would like is the 2011 panel to be ordered in descending order as is

Hi, Sorry for my ignorance, but how can I see if a function is 'compiled code' or plain R? E.g. the daisy function from the cluster package.

Hi everyone, I'm trying to encode/decode with speex, when I do not, the audio is loud and clear, but when I encode/decode to test audio quality, I get a really poor audio quality and a robotic sound. Here's my init audio method : #define AUDIO_QUALITY 10 - (void) initAudio { try { //SPEEX CONFIG speex_bits_init(&bits_in);

Dear R list., I am running a script to get a compact letter display. library(lme4) library(multcomp) library(gplots) ####### Mixed Effects Model ######### data <- read.table("AJmix.txt",header=TRUE, sep="\t") attach(data) y<-cbind(positive,negative) treatment<-factor(treatment) mouse<-factor(mouse) data$obs<-1:nrow(data) names(data) detach(data) attach(data)

Hello all I have a mail store on a volume replica 3 with no arbiter. A while ago the disk of one of the bricks failed and I was several days late to notice it. When I did, I removed that brick from the volume, replaced the failed disk, updated the OS on that machine from el8 to el9 and gluster on all three nodes from 10.3 to 11.1, added back the brick and started a heal. Things appeared to work

Dear all, I need some help on plotting multiple boxplots on one figure. I have three matrix A, B and C. Each of them is a 1000 by 10 matrix. The 10 columns of all three matrix correspond to the 10 values of the same parameter, say k=1, ..., 10. I want to make a plot where x axis represents different values of k. For each k value, I want to plot three boxplots, one on top of another. For

Dear All Can you please do me favor and let me know how can I compare two files but not in line-by-line basis on my CentOS server ? I mean say row#1 in file1 has the same data as say row#5 in file2 , but the comm compares them in line-by-line basis that is not intended . It seems that the diff cannot do the job as well . Thank you in advance -------------- next part -------------- An HTML

I am using the following function to plot columns from a CSV-file: plot_mi_time = function(mdata, miname) { mdata2 = mdata[row.names(mir_test) == miname, ] # print(mdata2) plot(c(1:3), mdata2, xaxt="n", ylab="Expression", xlab="Time(h)", , main=miname, pch=16) axis(1, at=c(1,2,3),labels=c("a","b","c")) } No I have a CSV file