search for: datai

Displaying 13 results from an estimated 13 matches for "datai".

Did you mean: data
2012 Jul 03
2
EM algorithm to find MLE of coeff in mixed effects model
...i b.list=list() for (i in 1:n) { b.list[[i]]=t(t(halfb)+miu.m[i,]) } #generate a list, each page contains Xi,yi,Zi, data.list=list() for (i in 1:n) { data.list[[i]]=list(Xi=Xi,yi=y.m[i,],Zi=Zi,b=b.list[[i]]) } #update sigma^2 t1=proc.time() tempaa=c() tempbb=c() for (j in 1:n) { #tempaa[j]=Eh4new(datai=data.list[[j]],weights.m=weights.mat,beta=beta.old) tempbb[j]=Eh4newv2(datai=data.list[[j]],weights.m=weights.mat,beta=beta.old) } var.new=mean(tempbb) if (ECM==T){var.old=var.new} sumXiXi=matrix(rowSums(sapply(data.list,function(s){t(s$Xi)%*%(s$Xi)})),ncol=ncol(Xi)) tempbb=c() for (j in 1:n) { t...
2012 Jul 03
0
need help EM algorithm to find MLE of coeff in mixed effects model
...i b.list=list() for (i in 1:n) { b.list[[i]]=t(t(halfb)+miu.m[i,]) } #generate a list, each page contains Xi,yi,Zi, data.list=list() for (i in 1:n) { data.list[[i]]=list(Xi=Xi,yi=y.m[i,],Zi=Zi,b=b.list[[i]]) } #update sigma^2 t1=proc.time() tempaa=c() tempbb=c() for (j in 1:n) { #tempaa[j]=Eh4new(datai=data.list[[j]],weights.m=weights.mat,beta=beta.old) tempbb[j]=Eh4newv2(datai=data.list[[j]],weights.m=weights.mat,beta=beta.old) } var.new=mean(tempbb) if (ECM==T){var.old=var.new} sumXiXi=matrix(rowSums(sapply(data.list,function(s){t(s$Xi)%*%(s$Xi)})),ncol=ncol(Xi)) tempbb=c() for (j in 1:n) {...
2000 Oct 12
2
works in R-1.1.1 but not in R-development; why?
...1.1.1 gives errors with R-devel; my (mis?)understanding of scoping rules is that it should have worked in both. The problems seem related to using the name of the data frame for extracting weights or subsets within a function call. The problems can be reproduced as follows: ********************** datai <- data.frame( y = rnorm(10), x1 = rnorm(10), x2 = abs(rnorm(10)), x3 = rep(seq(1,5),2), counter = rep(c(1,2),c(5,5))) formula <- as.formula(y ~ x1) # the following fails in R-1.2.0 but not in R-1.1.1 # > Error in eval(expr, envir, enclos) : Object "datos"...
2005 May 25
3
Problem with systemfit 0.7-3 and transformed variables
The 'systemfit' function in systemfit 0.7-3 CRAN package seems to have a problem with formulas that contain transformed (eg. log) variables. If I have my data in a data frame, apparently systemfit doesn't "pass" the information of where the variables should be taken to the transforming function. I'm not entirely sure if this is a bug or just a limitation, I was just
2000 Jun 02
3
graphical parameters in plot
I'm using plot to make a simple plot but I want to control where the tick marks go on bot axes. The graphical parameters xaxp and yaxp seem to be ignored by the plot function so I tried setting them using op <- par(no.readonly=T) par(xaxp=c(-2.4,-2.2,5), yaxp=c(-2500,10000,6)) plot(...) par(op) but they are still ignored. Next I tried to use the axis function as follows op <-
2008 Mar 19
2
nil-error
...key, unless of course it''s likely to be non-unique and cause a bug. Anyone have an idea on what could be causing this? I''m running svn revision 314. Thanks, Ian --- backgroundrb_debug.log --- Client disconected 000000079^D^H{ : type:^Ldo_work:^Kworker:^Tagronomy_worker: datai^B<9E>^A:^Rworker_method:^Mtransmit typedo_workdata414workeragronomy_workerworker_methodtransmit undefined method `send_request'' for nil:NilClass /vendor/plugins/backgroundrb/framework/packet/packet_master.rb:43:in `ask_worker'' /vendor/plugins/backgroundrb/server/master_worke...
2009 Sep 20
1
Return a list from a .Call but segfaults
Hello, I call a function via .Call passing to it a raw vector(D) and an integer(I) The vector is a series K1,KData1, V1,VData1, K2, KData2, ... where the integer K1 is the length of Data1 and similarly for Ki (wrt Datai)(similarly for V*) There 2*I such pairs( (Ki,KDatai), (Vi,VDatai)) The numbers Ki(and Vi) are written in network order. I am returning a list of I elements each element a list of two elements corresponding to KData and VData (which are serialized R objects) . ll <- .Call("returnLis...
2005 Sep 23
2
Strange behaviour of as.Date function
Dear All, I'm happily extracting data of temperature from an oracle db under R via RODBC. After manipulating the extracted data I put them into a data.frame 'dati' which is as follows: > dati DATA tm. UDINE/RIVOLTO tm.TORINO/CASELLE 1 2005-07-01 22.35 23.80 2 2005-07-02 22.70 22.85 3 2005-07-03 23.80
2007 Apr 18
2
virtual device bus "xen bus"
Hi, in order for installers and hardware discovery programs to be able to load the right device drivers, we really want a virtual device bus in Xen, a'la power5 VIO. This "xen bus" is on the top of Jeremy's TODO list, but since I have heard rumors that some of our friends here might have (some of the) code already it would seem prudent to check whether that code could be
2004 Oct 08
1
reading partial file content
...6 20 > > 4 4 > > 20 18 > > > > As per the example I have two columns and 7 rows of data ineach. > > However is real life data situation I may not know how many columnsare > > present and how rows are present and also with the certain datais > > missing. Yes I am assuming the data is delimited myTab. > > > > > > My question or rather problem is I want read data from each colum saycol > >x > > (8,11,16,18??EUR??.20) and store it into a variable so that I couldperform > >some > > operation...
2007 Apr 18
2
virtual device bus "xen bus"
Hi, in order for installers and hardware discovery programs to be able to load the right device drivers, we really want a virtual device bus in Xen, a'la power5 VIO. This "xen bus" is on the top of Jeremy's TODO list, but since I have heard rumors that some of our friends here might have (some of the) code already it would seem prudent to check whether that code could be
2007 Dec 25
6
error: no marshal_dump is defined for class Thread
I get this error periodically in backgroundrb.log, but only on my production box (the one running ruby 1.8.5) - my dev box doesn''t get this error. As far as I can tell, everything is working ok, despite the error (I don''t really know how to determine if the jobs are getting threaded successfully, but they are running, that I am sure of). Here''s the dump: no
2004 Oct 07
3
Read.Table Reading a Text file
Dear R users and Helpers I am beginner with using R and interested in carrying out certain task for my statistical research. I am reading data for a text file, which could contain data in following pattern x y 8 10 11 14 16 16 18 15 6 20 4 4 20 18 As per the example I have two columns and 7 rows of data in each. However is real life data situation I may not know how many columns are present