similar to: No error message no display output

Hi all, I want to make the following: I want to run a linear regression on each column of a matrix "estima" on the correspondent column on the matrix "estima2". You see I want to regress estima[,1] on estima2[,1] this way to all columns.... At the same time I want to make a regression adding each time a new observation. You see, the first regression will regress only one

Hi all, I´m trying to use write function to save the output of a program (my constructed "H" matrix) randz<-matrix(rnorm(1000000),500,2000) H<-matrix(0,500,2000) H[1,]<-randz[1,] for (j in 1:2000){ for (i in 2:500){ if(i<251) H[i,j]<-0.6*H[i-1,j]+randz[i,j] else H[i,j]<-H[i-1,j]+randz[i,j] }} write(H, file = "datad.txt",2000) If I ommit the 2000 on

If I do > example(lm) ... > mycoef <- function(object, ...) UseMethod("coef", object) > mycoef(lm.D9) Error in mycoef(lm.D9) : no applicable method for "coef" which is pretty surprising, as coef has a default method. After a bit of digging, this comes from do_usemethod having defenv = environment where the generic was defined */ defenv =

Hi I want to simulate random numbers normal distributed with this size (2000,10000). I tried this but my computer exhaust there is a fast way to make it? randz<-matrix(rnorm(2000000),2000,10000) [[alternative HTML version deleted]]

Hi I wrotte this function but when I get the "tmp.xls" file it shows data in a rare way. I mean not appears a matrix with 2000 rows and 100 columns. Can anyone help me, guide me? kim<-function(){ tmp<-randz<-matrix(rnorm(200000, mean=0,sd=0.01 ),2000,100); dim(tmp) write(tmp,file="tmp.xls") Thanks in advance. On the other hand, everytime I execute a function

Hi, I'm using R to do a time series analysis. In the model, I use the lags of some variables. such the lags of the variables have different length, I just can't use them directly in the lm function. Intuitively, I feel that "subset" might be useful, but I do not know how to use it. Anyone can give me an example syntax? Thanks. Harry [[alternative HTML version deleted]]

Hi all, I have this information on a file ht.txt, imagine it is a data frame without labels: 1 1 1 8 1 1 6 4 1 3 1 3 3 And on other table called "pru.txt" I have sequences similar this 4 1 1 8 1 1 6 4 1 3 1 3 3 1 6 1 8 1 1 6 4 1 3 1 3 3 1 1 1 8 1 1 6 4 1 3 1 3 3 6 6 6 8 1 1 6 4 1 3 1 3 3 I want to now how many positions are identical between each row in pru compared with ht. n and m

Hi all, I suppose this is a very simple question, but as I've lost already a bit of time with it, without being able to get what I wanted, I'm addressing the question to the group in the hope someone can help me. I pretend to plot the richness of herbaceous species (RichHN) as a function of time since remnant isolation (Isol) conditioned to the area of the remnant (fArea - this is a

Een ingesloten tekst met niet-gespecificeerde tekenset is van het bericht gescrubt ... Naam: niet beschikbaar Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070820/920567a6/attachment.pl

Hello! I'd like to create an assets-based economic indicator using data from a national household survey. The economic indicator is to be the first principal component from a principal components analysis, which (given the source of the data) I believe should take in consideration the sampling weights of the observations. After running the PCA with svyprcomp(), from the survey package, I

Any better (more efficient, built-in) ideas for computing coef[1]+coef[2]*x+coef[3]*x^2+ ... than polynom <- function(coef,x) { n <- length(coef) sum(coef*apply(matrix(c(rep(x,n),seq(0,n-1)),ncol=2),1,function(z)z[1]^z[2])) } ? Ben -- 318 Carr Hall bolker at zoo.ufl.edu Zoology Department, University of Florida http://www.zoo.ufl.edu/bolker

Hola, tenemos como equivalencia no paramétrica al Anova de 1 vía el test Kruskal Wallis, pero ¿existe un test no paramétrico equivalente al Anova de 2 vías? Gracias. Saludos. -- Sec. Apoyo Estadístico. Servicio de Apoyo a la Investigación (SAI) Vicerrectorado de Investigación e Internacionalización. Universidad de Murcia Edif. SACE. Campus de Espinardo. 30100 Murcia @. mariaelvira.ferre en

Hi all, I want to make multiple least squares estimation on two matrix (without intercept) imagine matrix "a" and matrix "b", I want to "regress" each colum on matrix "a" (dependent variable) on each column of matrix b, I mean, regress first colum on a to first column on b. Second column on a to second column on b. Imagine "a" and "b"

Hi, I''m having a problem with dates (laugh if you must :). $ rails --version Rails 2.3.5 $ ruby --version ruby 1.8.7 (2010-01-10 patchlevel 249) [x86_64-linux] I have a bunch of data in a text file that I''d like to place in my DB but am ending up with odd results when I save the data. Here''s my model: class CreateWaits < ActiveRecord::Migration def

excuse for to find a desperate solution... but i'm boried to spent hours... ;) not in asterisk !!! ;) i use the festial-script.pl of Donny Kavanagh... but i want to change the language that festival uses, depending on a variable for the callerid. English/Spanish How i can tell the script the correct voice that festival needs to use?? like --language spanish --language english ...in normal

I'm having trouble understanding how functions in the subset argument for lmList search for the objects they need. This trivial example produces "Error in rownames(fakedf) : Object "fakedf" not found": library(nlme) fitbyID <- function() { fakedf <- data.frame(ID = gl(5, 10, 50), A = sample(1:100, 50), B = rnorm(50))

Hi all, When somebody leaves a message in the voicemailbox, is there a way to know the file name of it? I need to return the voicemail file name in the deadagi command. Thanks, Anahi _________________________________________________________________ Convierte las fotos que m?s te gustan en tu nuevo fondo de escritorio para el ordenador. Es f?cil y adem?s gratis

Dear List Has anyone used R to distnguish between alternative forecasting models? In particular is the Diebold Mariano test available for use within R. Any assistance would be greatly appreciated. Graham Leask Lecturer in Strategy Economics & Strategy Group Aston University Aston Triangle Birmingham B4 7ET Tel: 0121 204 3150 E Mail: g.leask at aston.ac.uk

Hi people, when I try to get the billsec in the dialplan, it is 0... but if after that I check the database, it is right (not 0). I'm trying to get it in the h extension, like: exten => h,1,Noop(End) exten => h,n,Noop(Time is ${CDR(billsec)}) .... Is it updated after the extension h is executed? In that case, how can I get the call duration in the h extension? Thanks, Anahi

hello theora-list, i'd like to edit (especially crop) theora files I downloaded here and there. But as far as I know, there is no theora editor. Is there a cli theora files editing ? How could I export .ogg files into whatever format that i can then edit with software such as avidemux ? thanks very much thomas