similar to: Aov: SE's for split plot

Hello, I'm not much of a programmer, and am trying to understand the workings of the function below called "RStatFctn" within this bootstrap procedure. RStatFctn is defined to have two arguments: "x", intended to be a data vector; and "d" intended to be an index (or so it looks to me). Later, "rnormdat" is created to be the data vector. However,

>Hello, > >I'm trying to introduce myself to ggplot2. I'm using syntax from the help file, but pasting in my own data. I don't understand the error message I'm getting. Can anyone clue me in? A Google search of this error statement didn't return anything I could recognize as useful. > >Thanks, Paul > > >> str(d.AD) >'data.frame': 9

Hello, I'm trying to install the package metRology from SourceForge. I save the zip file metRology_0.9-06.tar.gz to my Windows 7 machine, and try to install the package using the RGUI: Packages > Install packages from local zip files > metRology_0.9-06.tar.gz. There's no .zip extension, but R seems to go to work on the installation with a couple warning messages and one error.

Hello, very basic question from a user who is baffled by the workings of computers in general When logging off R, a dialog box asked if I wanted to save my log, I chose yes. Then I noticed that the following message appeared in the Command Window Error in saveLog(currentLogFileName) : unused argument(s) ("C:/Documents and Settings/prewpj/My Documents/Data/Analyses/Healthcare/Hand Care

Hello, I?m wondering if anyone can offer advice on the out-of-memory error I?m getting. I?m using R2.12.2 on Windows XP, Platform: i386-pc-mingw32/i386 (32-bit). I am using the quantreg package, trying to perform a quantile regression on a dataframe that has 11,254 rows and 5 columns. > object.size(subsetAudit.dat) 450832 bytes > str(subsetAudit.dat) 'data.frame': 11253 obs.

Hello, I am not having success in a simple task. Using the effects package, I would like to add reference lines at probability values of 0.1 – 0.6 on a plot of the effects. The plot command works, but following up with an abline command produces the message “plot .new has not been called yet”, and of course the reference lines were not added. Looking through past R help lists, there was a

Hello, I’m trying to follow the syntax of a script from a journal website. In order to create a regression formula used later in the script, the regression matrix must have column names “X1”, “X2”, etc. I have tried to assign these column names to my matrix ScoutRSM.mat using a for loop, but I don’t know how to interpret the error message. Suggestions? Thanks, Paul

Hello, I am looking for a sample size function for samples sizes, to test proportions that are not binomial proportions. The proportions represent a ratio of (final measure) / (baseline measure) on the same experimental unit. Searches using RSeek and such bring multiple hits for binomial proportions, but that doesn't seem to fit my situation. Perhaps there's some standard terminology

Hi All, I would like to use meta.summaries from package rmeta to do a meta-analysis. I have available effect sizes as r's (which could be easily transformed to effect sizes in terms of d's). My problem is that I'm not sure what the se's of these r's should be ... The r-values are themselves computed from F-tests and t-tests for various studies. Are there R-functions that

Full_Name: Kirk Hampel Version: 2.4.1 OS: Windows Submission from: (NULL) (144.53.251.2) Given an AR(p) model, the last p SE's are wrong. The source of the bug is that the C code (ver 2.4.0) assumes *npsi is the length of the psi vector (which is n+p), whilst the predict.ar function in R passes out as.integer(npsi), where npsi <- n-1. Some R code following reproduces the error. Let p=4,

Hello All, Does anyone know of a package that performs constraint-based clusters? Ideally the package could perform "Time-Ordered Clustering", a technique applied in a recent journal article by Runger, Nelson, Harnish (using MS Excel). Quote, "in our specific implementation of constrained clustering, the clustering algorithm remains agglomerative and hierarchical, but observations

dear R wizards: apologies for two queries in one day. I have a long form data set, which identifies about 5,000 regressions, each with about 1,000 observations. unit date y x 1 20060101 <two values> 1 20060102 <two values> ... 5000 20081230 <two values> 5000 20081231 <two values> I need to run such regressions many many times, because they are part of an

Dear all, I would like to characterize the behavior of a certain product over its lifetime. The purpose is to compare the behavior of experimental versions of the product vs. a baseline version. It's known that the behavior of the product follows a roughly trapezoidal pattern: start-up effect where the performance ramps up, then steady state performance, finally a degradation at the end of

I?m entirely stumped on this particular issue and really hoping someone has some advice for me. I am running a covariant model in lm I would like to give the standard errors or the confidence intervals for the fitted lines. I?ve been using the dataset OrangeSprays where I want lines for each level of treatment over the covariant ?colpos?. I?ve been able to calculate intercepts and slopes for

Hello all, I am unsure of how to interpret the output from a Generalized Estimating Equation analysis of an ordinal response. I hope someone can enlighten me. The analysis was done using package 'repolr'. The data consists of a Score on a 3-point scale from 56 Subjects after repeatedly washing their hands with soap. Two soap Products were tested, each panelist washed 10 times = 10

I am trying to make predictions with se's using a nlme (kew11.nlme below). I get an error indicating levels for a factor are not allowed. I have searched and read Rnews, MEMSS, MASS, R-Help, and other lists in Spanish where I found questions similar to mine but not solution. I do not really care about the method used. Any suggestions to obtain predictions with se's from an nlme

Hi John: Thanks! Below is the data using your suggestion. I used "ggplot" to make a graph. I am not too happy with it. I am looking for something simpler and cleaner. Plot is attached. I also tried "lattice" package, but nothing got plotted with "xyplot" command, because it is looking for a numeric variable on x-axis. ggplot(TrialData4, aes(x=Income, y=Percent,

Btw, I think "lattice" graphics will provide a better solution than "ggplot", because it puts appropriate (space saving) markers on the axes and does axes labels well. However, I cannot figure out how to do it in "lattice". On Thu, 6 Jul 2023 at 15:11, Anupam Tyagi <anuptyagi at gmail.com> wrote: > Hi John: > > Thanks! Below is the data using your

Dear all, Myself Vijaykumar Muley working as senior research fellow. By training I am a computational biologist with not a strong knowledge of statistics. I have done some analysis which is explained as follows, I have 10340 (X) profiles of binary vectors with same length(N=845), I will call then "gene profiles" for example... v1 v2 v3 v4.....vN a 1 0 1 0 1 b 0

Hallo Anupam I do not see much difference in ggplot or lattice, they seems to me provide almost identical results when removing theme part from ggplot. library(ggplot2) library(lattice) ggplot(TrialData4, aes(x=Income, y=Percent, group=Measure)) + geom_point() + geom_line() + facet_wrap(~Measure) xyplot(Percent ~ Income | Measure, TrialData4, type = "o", pch = 16, as.table =