similar to: Need some hint on faster data manipulation.

Hi Can anybody please suggest me about the documentation of "factanal.fit.mle()" (Not factanal()------ searching factanal.fit.mle() in R always leads to factanal()). Is there any function for doing principal component factor analysis in R. Regards Souvik Bandyopadhyay JRF, Dept Of Statistics Calcutta University [[alternative HTML version deleted]]

R Gurus,I have a vector of nearly 90,000 characters from which I have to extract the index of the characters which are repeated. So suppose if x<-c("a","a","b","a","b","c","d") then my output would be a vector having the index where the values are repeated i.e (1,2,3,4,5). I have been able to isolate out the values that are

Hi Is there any way to fit a bivaraite longitudinal mixed model using R. I have a data set with col names resp1 (Y_ij1), resp2 (Y_ij2), timepts (t_ij), unit(i) j=1,2,..,m and i=1,2,..n. I want to fit the following two models Model 1 Y_ij1, Y_ij2 | U_i = u_i ~ N(alpha + u_i + beta1*t_ij, Sigma) U_i ~ iid N(0, sigu^2) Sigma = bivariate AR structure alpha and beta are vectors of order 2.

Hi, Using indicator variables I have been able to fit and run the code for fitting a bivariate mixed model using unstructured covariance matrix The code is lme.fit1<- lme(one.var~-1+indic1+indic2+I(indic1*d.time)+I(indic2*d.time), random =~ -1+indic1+indic2|m.unit, weights = varIdent(~1|indic1) ,data = new.data) My variables are one.var :- the two response variables stacked one after

Hi, I have a bivariate longitudinal dataset. As an example say, i have the data frame with column names var1 var2 Unit time trt (trt represents the treatment) Now suppose I want to fit a joint model of the form for the *i* th unit var1jk = alpha1 + beta1*timejk + gamma1* trtjk + delta1* timejk:trtjk + error1jk var2 = alpha2 + beta2*timejk + gamma2* trtjk + delta2* timejk:trtjk +

dear list, as a part my problem. I have to estimate some parameters using ML estimation. The form of the likelihood function is not straight forward and I had to use a for loop to define the function. I used "optim" to maximise the result but was not sure of the programme. To validate my results, I tried to write a function to obtain the MLE of a bivariate normal in the same manner. On

Hi, all. I'm trying to cluster 12,500 objects using hclust from package mva. The distance matrix takes up nearly 600 MB. The distance matrix also needs to be copied when being passed to the fortran routine that actually does the clustering (it's modified during the clustering), so that's 1200 MB. I'm actually on a machine with 2.5 GB of memory (and nothing else running), so I

I have a data matrix with genes as columns and samples as rows. I want to create all possible gene ratios.Is there an elegant and fast way to do it in R and write it to a dataframe? Thanks for any help. Som. -- View this message in context: http://r.789695.n4.nabble.com/Calculating-all-possible-ratios-tp4627405.html Sent from the R help mailing list archive at Nabble.com. [[alternative HTML

Hi, I am trying to run libvchan on ARM Fast Model for VExpress platform for using inter domain communication. But there are some errors. I am attaching the strace outputs for vchan-node1 run in dom0 and domU. Has anybody faced these issues before. Xl console 1The command is Usage: vchan-node1 [client|server] [write|read] domid nodepath my nodepath is - /local/domain/1/data/vchan-2-1 Regards

Hi Shakeel, Assuming that you are starting with a bunch of dates: # make a vector of character strings that can be converted to dates rep_dates<-paste(sample(1:30,500,TRUE),sample(1:12,500,TRUE), sample(2013:2017,500,TRUE),sep="/") # if this isn't your format, change it date_format<-"%d/%m/%Y" # create a data frame with a column of dates

Once last night and the night before a XenU go into a Zombie state, requiring a reboot. I''m not quite sure what is happening, and am looking for advise on how to diagnose the problem. .. While writing this email, it crashed again. This time I had an ''xm console calcutta'' capturing the output. I read a suggestion in these archives suggesting that I can just

Folks, Several of you have received rather blunt emails from me when you sent me an off-list question or have continued an on-list thread off-list. It''s not that I''m trying to be a jerk but please consider: - There have now been over 50,000 downloads of Shorewall. - There are over 1,000 subscribers to this list. Assuming that each person downloading Shorewall has downloaded

Dear List: I have a data manipulation problem that I was unable to solve in R. I did it in SQL, and it may be that the solution in R is to do it in SQL, but I wondered if people could imagine a vector-based solution. Imagine a list A[i] of observers who observe some set of events B[j]. Each observer i may observe one or more events, and each event j may have been observed by one or more

Has anyone figured out the way of getting the caller id for BSNL on Asterisk 1.4.2 I have tried following link http://bugs.digium.com/view.php?id=6683&nbn=24 but was not able to get it, although did not ge any error too. I always get the caller id as asterisk. Can someone please help. Regards, Sanjay Rajdev

Hello, is it possible to rout ALL RTP Data over Asterisk, like SIP1 <---RTP---> Asterisk <---RTP---> SIP2 I know it seems quite useless. But I want to simulate a IAX -> SIP connection and have no Phonecard installed on my computer ;) Thanx, Kalle

Hello all, I have a number of Polycom phones 601's and 430's and I'm seeing: Got SUBSCRIBE for extensions without hint. Please add hint to s to context local-hints in the CLI over and over. I have: [local-hints] exten => 110,hint,SIP/110 exten => 111,hint,SIP/111 exten => 112,hint,SIP/112 exten => 113,hint,SIP/113 exten => 114,hint,SIP/114 The hints seem to be

Hi Centos Team, I was here the hard days in which CentOS 6 appeared would never appear, and then was CentOS 5.6 (maybe 5.7...do not remember well) that struggled to see the light and I was very afraid CentOS would not survive... but as I was worried in the past actually I'm very impressed of the speed of latest releases: what has changed? Just out of curiosity... Thanks for your hard work

Respected Sir, With reference to my mail to you dated 8th January,2013, and the reply by you dated 9th January, 2013, I am sending this mail to you. I had a problem regarding running a program in the latest version of the "BRugs" package in R 2.15.1 and 2.15.2. I want to mention here that this program runs well to others, who are running it using the earlier version

Hi List, I am confuse about this feature. When we use Wait(20) in active call session then it's work. But when we use after hangup the call then Asterisk don't wait from define time. Ex:- [call_log] exten => 4368,1,Answer() exten => 4368,n,Flite("Welcome") exten => 4368,n,Set(__StartTime=${STRFTIME(${EPOCH},Asia/Calcutta,%Y-%m-%d %H:%M:%S)}) exten =>

I have two server Server_A(outbound call) for agent login and agent make a outbound call from here and pass into server Server_B call extension.conf exten => _91XX.,1,AGI(agi://127.0.0.1:4577/call_log) exten => _91XX.,n,Dial(SIP/${EXTEN}@192.168.53.197,,tToR) exten => _91XX.,n,hangup() Server_B[192.168.53.197] for call forwarding extension.conf exten =>