search for: lms

...f the intercepts in one array). I figured I'd put the model objects into an array, and then (somehow) I could easily create corresponding arrays of intercepts or residuals or whatever, but I can't the parameters back out. Right now I've stored the model objects in a 2-D array: > lms.ASP <- array(list(), c(3,4)) Then I fill the array element-by-element: > surf105.lm. ASP <- lm(ASP ~ time) > lms.ASP[1,1] <- list(surf105.lm.ASP) Something is successfully being stored in the array: > test <- lms.tx.ASP[1,1] > test [[1]] Call: lm(formula = ASP ~ time...

...- c(4,17,12,12,17,3,9,22,13,20) condition <- c(0.72, 0.72, 0.67, 0.73, 0.76, 0.65, 0.68, 0.78, 0.73, 0.71) test <- data.frame(individual, day, condition) ind.id <- unique(test$individual) ind.list <- lapply(1:length(ind.id), function(i){ subset(test, test$individual==ind.id[i])}) lms <- lapply(ind.list, lm, formula=condition~day) ldply(lms, function(x) x$coefficients) here I can display coefficients, here I need to write code for r2.squared I tried with following script summary(ldply(lms, function(x) )$r.squared, $p value) but it did not work. can any one help me? thanks...

Hello, I'm trying to create birth centile charts from about 90000 birth weights using the package lmsqreg. As this is not a CRAN package I have problems finding any examples, which explain how to estimate the LMS parameters needed for the equivalent degrees of freedom. (I only have some basic understanding of statistics, I roughly understand the formulas, but then I don't know how to recreate t...

...73, 26.417, 24.868, 29.895, 24.200, 23.215, 21.862, 22.274, 23.830, 25.144, 22.430, 21.785, 22.380, 23.927, 33.443, 24.859, 22.686, 21.789, 22.041, 21.033, 21.005, 25.865, 26.290, 22.932, 21.313, 20.769, 21.393) and I would like to fit the following model: mod1 <- lqs(y ~ x1 + x2, method="lms", nsamp="exact") mod1$coefficients (Intercept) x1 x2 35.5293489 0.4422742 -1.2944534 mod1$bestone [1] 12 17 27 Now, instead of using the formula, I want to provide the design matrix and the response, then: X <- cbind(1, x1, x2) mod2 <- lqs.default(X, y,...

Hi, I have been using the lqs function with method='lms'. However the results I get are a little different from the results noted by Rousseeuw & Leroy (Robust Regression and Outlier Detection) and I was wondering how to use these results for outlier detection. I'm using the stackloss dataset, for which the original Rousseeuw et al. program...

Everything is ok on Firefox, IE, and iPad. Thanks, Erin On Tue, Jan 8, 2013 at 7:58 PM, Erin Hodgess <erinm.hodgess at gmail.com> wrote: > Dear R People: > > Has anyone used R2HTML in web files that were on the Blackboard LMS, please? > > I'm starting to do these but wanted to know if there were any > potential pitfalls. > > Thanks, > Erin > > > -- > Erin Hodgess > Associate Professor > Department of Computer and Mathematical Sciences > University of Houston - Downtown > mai...

...the random seed for. Is it a starting point for iterations? Or do different results occur because of the estimation doesn't converge? I tried data "phones" from package MASS. You find this example as well in the MASS-book on page 160. > lqs(calls~year, data=phones, method="lms") Call: lqs.formula(formula = calls ~ year, data = phones, method = "lms") Coefficients: (Intercept) year -55.947 1.155 Scale estimates 0.9377 0.9095 > > lqs(calls~year, data=phones, method="S") Call: lqs.formula(formula = calls ~ year, dat...

Dear R People: Has anyone used R2HTML in web files that were on the Blackboard LMS, please? I'm starting to do these but wanted to know if there were any potential pitfalls. Thanks, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: erinm.hodgess at gmail.com

...et?s say I have 10 models, each named m1,m2,m3..., and I would like to summarize them automatically and simultaneously - e.g., to extract parameter estimates later on from all models; how can I do that? I have tried: x=1:10 #this creates some example data y=rnorm(10) m1=lm(x~y) m2=lm(x~1) sum.lms=function(x)summary(paste("m",x,sep="")) sum.lms(1:2) but the paste() function always returns "m1" "m2" in quotation marks instead of m1 m2 Do you know a way around this? I guess I would need to use expression() or substitute() but I don?t know how to const...

...b/pdb/mysql.so:mysql mysql:mysql host = localhost mysql:mysql user = samba mysql:mysql password = ambas mysql:mysql database = samba root@auth /usr/src/installs/samba-3.0.0/examples/pdb/mysql# It booted up correctly but when i try to had a user i get: root@auth /usr/local/samba/bin# ./smbpasswd -a lms New SMB password: Retype new SMB password: No builtin nor plugin backend for plugin found Loading plugin:/usr/local/samba/lib/pdb/mysql.so:mysql failed! Failed initialise SAM_ACCOUNT for user lms. Failed to modify password entry for user lms root@auth /usr/local/samba/bin# Can anybody tell me how...

Hi all, We all know that mismatch between clocks of ADCs of far-end voice and near-end voice is not allowed in a time-domain or frequency-domain LMS based AEC system. It means that capture and render audio streams must be synchronized to a same sample rate. However, I found that this restriction is removed in microsoft AEC from Windows XP SP1. Anyone knows how microsoft AEC do it? This technology is much helpful for us to implement AEC in commo...

...We are collection data of Age, Weight and Length. The data are used to produce diagnosis-specific Growth charts (like the CDC Growth Charts: http://www.cdc.gov/nchs/data/nhanes/growthcharts/set1clinical/cj41l017.pdf) Data-example is in the UsingR package as data(kid.weights). In the CDC-papers a LMS approch is described. Cole has programmed the methology (I guess) in a <R 2.0.0 package called lmsqreg: """ > library(lmsqreg) Error in library(lmsqreg) : lmsqreg' is not a valid package -- installed < 2.0.0?) """ Where to go forward in R. Thanks for any...

...ning data), and I'left 10 as test data. Training data and test data are stored in an excell. training <- read.xls("C:/...../training.xls") , the same for test. That is: v1 v2 ... v15 When I type str(training) and str(test), both sets have the same names The resulting model is lms <- lm(vd ~ log(v1) + fv2+ fv5+ fv7 ) -fvi means they were turned into factors- plms<- predict(lms, new=test ,interval="prediction") Error at model.frame.default(Terms, newdata, na.action = na.action, xlev = object$xlevels) : length of the variables are different (found for ...

...rts with aa* and the last one which appears is "zmonteiro" (our last user). Now, after some more tests, i found out that this is only happening on winXP professional SP2 and win2003 server. I took some screenshots so anybody that is willing to help me can check out my problem: 1? http://lms.ispgaya.pt/goodies/sambaXP/erro_XP_cant_find.JPG 2? http://lms.ispgaya.pt/goodies/sambaXP/jaam_windows_2000.jpg on the first image we can see a logged in user on my domain (jaam) and, allthough c:\documents and settings\jaam\Desktop got automatically assigned ISPGAYA\jaam as its owner, that user c...

...odels, IRLS and vector smoothing. It is a large R package that fits many models and distributions, usually by maximum likelihood estimation. See below for a partial listing, or better, http://www.stat.auckland.ac.nz/~yee/VGAM/doc/VGAMrefcard.pdf. Some broad categories are categorical data analysis, LMS quantile regression, extreme value models, univariate and continuous distributions. The package is under continual development and definitely not complete yet. Until then it should probably be considered 'alpha software'. Nothing is cast in iron yet, so please expect lots of changes if you...

Hi everybody, We have FLAC supported for input/output in REAPER (http://www.reaper.fm), and the problem is that when writing 24 bit FLAC files, the data isn't compressed (i.e. the FLAC is slightly larger than if it was writing to .WAV). The files play back fine, however, and 16 bit mode works great. We're using flac-1.1.2, built on win32 with MSVC6 w/ SP5 + VCPP, and NASM version 0.98.39

Hi, trying to do a robudt regression of a two-way linear model, I keep getting the following error: > lqs(obs ~ y + s -1,method="lms", contrasts=list(s=("contr.sum"))) Error: lqs failed: all the samples were singular Robust regression with M-estimators works (also regular least square fits, of course): rlm.formula(formula = obs ~ y + s - 1, method = "M", contrasts = list(s = ("contr.sum")))...

Hi All, Many Thanks to Underwood for her excellent review of our big trouble which prevent LMS-based AEC algorithms to be used in most computer. Maybe it can be summaried as follows: 1. Different sample rate of sampling and rendering does exists in most low-cost soundcards (In my experiments over more than 20 soundcards, the differences range from 0.5Hz to more than 50Hz when sample rate is...

Hello, I have a double matrix that I want to represent in a line chart. Although I have seen some examples I still don't manage to get it. My data is this (a double matrix called mymatrix) : Blogs Wikis Redes Etiq. SPC LMS Menor de 30 57.14 28.57 14.29 28.57 57.14 28.57 de 31 a 40 63.83 61.70 29.79 17.02 59.57 70.21 de 41 a 50 72.64 70.75 47.17 20.75 55.66 75.47 Mayor de 51 62.07 58.62 48.28 17.24 50.00 67.24 As I understand from http://www.statmethods.net/graphs/line.html this example (see bottom of the page),...

...is cancelled than the one without resampling. But there is still echo can be heared. After all, frequency step of sample rate conversion is limited, mismatch is still exist after resampling. Someone told me that capture and render codec have different clock generator which shift independently. And LMS algorithm is very sensitive to the difference between sample rates. Sincerely Maoquan At 2011-04-12 21:46:26?"Shridhar, Vasant" <vasant.shridhar at harman.com> wrote: I would imagine that it is handle through basic asynchronous sample rate conversion. There is a lot of literat...