similar to: Time-Ordered Clustering

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'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 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, 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

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, 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 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

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

Hello, Can anyone explain why the following message appears for the function model.tables, where se=T? In the V&R MASS text, p.285 the se=T option works for a split plot example that seems similar to my operation. But the model.tables documentation, in the Arguments section for "se", states "should standard errors be computed?". "Warning: Warning in

>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 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

Hello, Does anyone know whether a Windows version of the Ratings package will be available? Thank you, Paul Prew CONFIDENTIALITY NOTICE: \ This e-mail communication an...{{dropped:11}}

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

Hello, I've been enjoying using the "Mixture and Hidden Markov Models in R" by Visser & Speekenbrink to learn how to apply these analyses to my own data using depmixS4. I currently have a fitted 4-state mixture model with three emissions variables and one binomial covariate (HS). I am trying to compute confidence intervals using the following code, where fmms4s is the model:

On Tue, 30 May 2023 17:43:31 +0000 Heather Lucas <hlucas2 at lsu.edu> wrote: > Hello, > > I've been enjoying using the "Mixture and Hidden Markov Models in R" > by Visser & Speekenbrink to learn how to apply these analyses to my > own data using depmixS4. > > I currently have a fitted 4-state mixture model with three emissions > variables and one

Hi All, Does someone have an idea of how to cleverly convert a categorical timeseries into a transition matrix? Ie, I have something like: x<- c(1,1,2,1,1,2,2,2,1,2), And I want a matrix with counts and/or probabilities: > tr <- matrix(c(2,3,2,2),2,2) > tr [,1] [,2] [1,] 2 2 [2,] 3 2 Meaning that there are two transitions from 1 to 1, two from 1 to 2, three from 2 to 1

Dear useRs, I'm trying to create a barplot like so: x=matrix(1:10,2,5) barplot(x,leg=c("left","right"),besid=T) The legend is placed in default position topright, however the data are plotted there too. I tried controlling the legend position by adding x="topleft" but this results in an error that x matches multiple formal arguments. Leaving out the legend

Dear List, We are using RHmm to cluster data through HMM. We would like to restrict the transition matrix of HMM, to get hierarchical connections between clusters. But, RHmm doesn't seem to support these restrictions. Can any one suggest a library to do that. Thanks, Kishor

Dear R users & Experts, This is just a curiousity, I was wondering why the dominant eigenvetor and eigenvalue of the following matrix is given as the third. I guess this could complicate automatic selection procedures. 0 0 0 0 0 5 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 Please

Hi All, I'm building a package using C/C++ and Fortran code which usually runs fine. However I do get occasional bus errors around the time of exiting one of the C functions. Where do I need to be looking to solve this problem? Do bus errors stem from unmapped memory exceptions? How can I find out if that is what is happening while running R? Any hints greatly appreciated. Best, ingmar