similar to: contrast matrix in package multcomp

I'm trying to get the following section of code to work, I think the problem is being caused by the assignment of data to the lm function not evaluating to "train" in the parent environment but I can't seem to figure out how to do this. fitmodel <- function(trial,data) { wrap.lm <- function(formula,data,...) { cat("in wrap lm",NROW(data),"\n");

Hi again, yesterday I mailed this query however I could not see this on the mail list. Therefore, I am reposting it again. I was using package.skeleton() function to create the skeleton of my package in windows. Here is my attempt: rm(list = ls()) setwd("F:/R_PackageBuild") package.skeleton("trial1", namespace = TRUE, code_files = "F:/R_PackageBuild/trial.r") In

i have a function that returns a list containing a variety of variable types i am trying to run the function multiple times and return the output into a variable with a semi-consistent naming pattern i.e., for ten trials i want to return the list into variables trial1,trial2,...trial10 is there a generic way to get this to happen i have a similar process that does the same thing to an external

Hi, I am wondering if anyone can suggest how to test the equality of 2 proportions. The caveat here is that the 2 proportions were calculated from the same number of samples using 2 different tests. So essentially we are comparing 2 accuracy rates from same, say 100, samples. I think this is like a paired test, but don't know if really we need to consider the "paired" nature of the

Hi again, I was using package.skeleton() function to create the skeleton of my package in windows. Here is my attempt: rm(list = ls()) setwd("F:/R_PackageBuild") package.skeleton("trial1", namespace = TRUE, code_files = "F:/R_PackageBuild/trial.r") In the trial.r file, there are 2 objects, one is a function and another is data. Here they are: fn1 <-

Hi every one, I am very very new to R and solicit your kind help I am trying to use my excel files in R. I used the Xlread/write package and able to load the a sample file into console. Then I tried to find the mean of a column vector for example. But I could not proceed Here is the problem rfile <- "Trial1.xls" > read.xls( rfile, colNames = TRUE, sheet = 1, type =

Hello all, I am low down on the learning curve of R but so far I have had little trouble using most of the packages. However, recently I have run into a wall when it comes to bootstrapping a Levene's test (from the car package) and thought you might be able to help. I have not been able to find R examples for the "boot" package where the test statistic specifically uses a

Dear R-users, I want to plot residuals vs fitted for multiple groups with ggplot2. I try this code, but unsuccessful. library("plyr") models<-dlply(dat1,"d",function(df) mod<-lm(y~x,data=df) ggplot(models,aes(.fitted,.resid), color=factor(d))+ geom_hline(yintercept=0,col="white",size=2)+ geom_point()+ geom_smooth(se=F) -- --- Catalin-Constantin ROIBU

Hi, It would be better to give an example. If your dataset is like the one attached: con<-file("Trial1.txt") ?Lines1<- readLines(con) ?close(con) #If the data you wanted to extract is numeric and the header and footer are characters, dat1<-read.table(text=Lines1[-grep("[A-Za-z]",Lines1)],sep="\t",header=FALSE) dat1 #?? V1 V2 V3 V4 V5 #1? 38 43 39 44 45 #2? 39

Dear R family,i am trying to plot and save, simultaneously, about 1000. the name of each plot is contained in "names" file. when i run this loop, i get an error. "Error in plot.new() : Unable to open file 'C:/R/SAVEHERE/myplot_Tak.jpg' for writing". could you please correct the mistake in the loop? >names<-(names(sp)) >for(a in seq_along(names)){ >mypath

Sorry I don't think my first trial worked I am just starting to use samba, and I keep running into the same problem time and again. I have installed SAMBA 2.0.2 correctly and have (i think) a correct smb.conf file. Now when I try to log on to one of my shares it will not let me in. I get the following errors when logging in: [dutchboy@lisa dutchboy]$ smbclient \\\\lisa\\dutchboy doing

Dear R-devel, I've written a numerical solver for SOCPs (second order cone programs) in R, and now I want to move most of the solver code into C for speed. I've written combined R/C packages before, but in this case I need to do matrix operations in my C code. As I have never done that before, I'm trying to write some simple examples to make sure I understand the basics. I am stuck

Hello I have data of body length and body weight of people of different skin colors. I tried to write a code to plot body length and body weight according to the skin colors. (Thanks for Petr's advice so far.) A loop is used but an error shows up in the following code. It says: unexpected '}' in "

Hello, I am using R 1.7.1 on a Redhat Linux machine, version 7.3. The following works fine: setClass("ok", representation( "A" = "matrix", "Cmatrix" = "matrix")) new("ok", "A" = diag(4), "Cmatrix" = diag(4)) But the following doesn't work: setClass("notok", representation(

Hello, happy new year to all. The new Year gave me a new computer with Win98 and promptly I installed R on it. I created a directory R, with subdirectories gcc, perl,bin, helpwrk and rw1040. I got rw1040 from BDR place, the rest were gotter from the internet in the last week of the year, i.e. they are the newest versions. I modified autoexec.bat to put gcc\bin, perl\bin, etc in the path. I put the

Dear r-helpers, I have a question about multiple testing. Here an example that puzzles me: All matrixes and contrast vectors are presented in treatment contrasts. 1. example: library(multcomp) n<-60; sigma<-20 # n = sample size per group # sigma standard deviation of the residuals cov1 <- matrix(c(3/4,-1/2,-1/2,-1/2,1,0,-1/2,0,1), nrow = 3, ncol=3, byrow=TRUE, dimnames =

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

Dear R-help list, I would like to perform multiple comparisons for lme. Can you report to me if my way to is correct or not? Please, note that I am not nor a statistician nor a mathematician, so, some understandings are sometimes quite hard for me. According to the previous helps on the topic in R-help list May 2003 (please, see Torsten Hothorn advices) and books such as Venables &

Dear list, As a reply to my recent mail: > simint and TukeyHSD work for aov objects. > Can someone point me to similar functions for lme objects? Douglas Bates wrote There aren't multiple comparison methods for lme objects because it is not clear how to do multiple comparisons for these. I don't think the theory of multiple comparisons extends easily to lme models. One could

I have run the following logistic regression model: options(contrasts=c("contr.treatment", "contr.poly")) m <- glm(wolf.cross ~ null.cross + feature, family = "binomial") where: wolf.cross = likelihood of wolves crossing a linear feature null.cross = proportion of random paths that crossed a linear feature feature = CATEGORY of linear feature with 5 levels: