similar to: regression using GMM for mulltiple groups

Dear all, I know there have been various questions posted over the years about loops but I'm afraid that I'm still stuck. I am using Windows XP and R 2.9.2. I am generating some data using the multivariate normal distribution (within the 'mnormt' package). [The numerical values of sanad and covmat are not important.] > datamat <-

Dear R-helpers, I am using a vrtest on time series data. My commands are as follows; read.table("B.txt",sep="\t",fill=TRUE, na.strings = "NA") require(vrtest) rm(list=ls(all=TRUE)) datamat <- read.table("B.txt",sep="\t",fill=TRUE, na.strings = "NA") column <- 1 nob <- nrow(datamat) y <-

A multinomial logit model can be specified as a conditional logit model after restructuring the data. Doing so gives flexibility in imposing restrictions on the dependent variable. One application is to specify a loglinear model for square tables, e.g. quasi-symmetry or quasi-independence, as a multinomial logit model with covariates. Further details on this technique and examples with several

Hello, I am trying to simulate 10 relicates of 100-tables. Each table is a 2 x 3 and 80% pf the tables are true nulls and 20% are non-nulls. The nulls follow the Hardy Weinberg distribution (ratio) 1:2:1. I have the code below but the p-values are not what I am expecting. I want to use the Cochran Armitage trend test to get the p-values. num.reps=10 num.vars=1000 pi0 = 80 num.subjects = 100

Hello, I want to draw a histogram of the mean of sample observations drawn from multivariate t distribution. I am getting the following error corresponding to the code I used. Though I am getting the graph, but I am curious to know the warning message. Warning messages: 1: In if (freq) x$counts else { : the condition has length > 1 and only the first element will be used 2: In if (!freq)

Hello, I'm working with R version 1.5.0 in Windows. I've written a function (SummaryMat, segment below) which uses a loop to repeatedly call another function (PercentsMat, segment below). PercentsMat creates a matrix and adds rows to it each time it is called. I use deparse(substitute(...)) to get the names of the lists sent to PercentsMat to use them as row names in the generated

Dear all, I have written a function for calculating the volume of a tree (=trad) or snag (=h?gst). The included volume regreesion model includes ten parameter values, which are tree species specific. bj?rk.formh?jd.pars is an object which includes the parameter values (parameter set) for birch (=bj?rk). There is one row per tree in the data object. > relev.kols[1:5,

Hi, I got a computational problem with lme (nlme library R 1.9.1) using corSymm(). Here is the data: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 0.19639793 0.09127954 0.11733288 0.07598273 0.06545106 0.06211532 [2,] 0.22773467 0.10981912 0.16052847 0.38101187 0.18353474 0.24072918 [3,] 0.46743388 0.45733836 0.32191178 0.43356107 0.39159746 0.53984221 [4,]

Dear all, I am trying to make a barplot with clustered pairs of bars, using class=numeric data and the following command: barplot(c(bline_precip[10,9], bline_runoff[10,9], cccma_precip[10,9], cccma_runoff[10,9], csiro_precip[10,9], csiro_runoff[10,9], ipsl_precip[10,9], ipsl_runoff[10,9], mpi_precip[10,9], mpi_runoff[10,9], ncar_precip[10,9], ncar_runoff[10,9], ukmo_precip[10,9],

Folks, I have been using the VAR {vars} program to find a fit for the following bi-variate time series (subset): bivariateTS<-structure(c(0.950415958293559, 0.96077848972081, 0.964348957109053, 0.967852884998915, 0.967773510751625, 0.970342843688257, 0.97613937178359, 0.980118627997436, 0.987059493773907, 0.99536830931504, 1.00622672085718, 1.01198013845981, 1.01866618122606,

Dear All, I've successfully import my synteny data to R by using scan command. Below show my results. My major problem with my data is how am i going to combine the column names with the data( splt) where i have tried on cbind but a warning message occur. I have realized that the splt data only have 5 column instead of 6. Please help me with this!! I want my data to be a numerical

-- ZeeShan Afzal (Channel Development Manager) Cell No. (0321) 844 6631. (0300) 844 6631. Pakistan Office Products (POP) UAN : 111 767 767 URL : www.pop.com.pk S-15, 2nd Floor, Hafeez Centre, Main Boulevard, Gulberg III, Lahore. Pakistan. Tel : (+92-42) 571-545, 571-5154. Fax: (+92-42) 571-5265. Email : ZeeShan@POP.com.pk E-mail Disclaimer Notice: This communication may contain

When I run the VB6 aplication , this is the error log wine "c:\silverlab\silverlab.exe" fixme:ole:CoRegisterMessageFilter stub fixme:ole:OleLoadPictureEx (0x4143f104,774,0,{7bf80980-bf32-101a-8bbb-00aa00300c ab},x=0,y=0,f=0,0x406bfae8), partially implemented. fixme:ole:OleLoadPictureEx (0x4143f104,1086,1,{7bf80980-bf32-101a-8bbb-00aa00300 cab},x=0,y=0,f=0,0x406bfab8), partially

Hello Chris,  I 've found two other issues  with MAE and CRPS, giving warning  when running examples. I've the same issue on my data.  Hope that you could find some time to take a look here. Thank you Anna > library(ensembleBMA) Loading required package: chron > example(MAE) MAE>   data(ensBMAtest) MAE>   ensMemNames <-

Hi I''m very keen to learn RoR and have been following the Agile Web Development with Rails book (first edition). However, I''ve come across a problem which has had me stumped for days now. After following the instructions up to page 85 I have tried testing the code thus far, as suggested, by adding an item to the cart. But I am getting the following error: --

I need to replicate MySQL DB of mulltiple server on SiteA to my DR-Site Site_B... all DB are alocated on RHEL,SuSE,Centos,Debian, FreeBSD servers. I need a script to take Multiple MySql DataBase Backup and then import to SiteB, the replica can be done as cold or hotbackup and cron it Thanks -- Madunix_at_Gmail Sysadmin "Computers are useless. They can only give you answers" - Pablo

I am using read.csv to read a CSV file (produced by saving an Excel file as a CSV file). The columns containing dates are being read as factors. Because of this, I can not compute follow-up time, i.e. Followup<-postDate-preDate. I would appreciate any suggestion that would help me read the dates as dates and thus allow me to calculate follow-up time. Thanks John John Sorkin M.D., Ph.D. Chief,

Hello -- The question I have is about the gmm() function from the 'gmm' package (v. 1.4-5). The manual accompanying the package says that the gmm() function is programmed to use either of four numerical solvers -- optim, optimize, constrOptim, or nlminb -- for the minimization of the GMM objective function. I wonder whether there is a way to pass controls to a solver used while calling

I’m starting to use the GMM estimator with panel data in R. I´ve read the document «Panel Data Econometrics in R: The plm Package» (Croissant and Millo). In Stata before using the functions lag() or diff() we must sort the data by individual and by time. I would like to know if I have to do something like this in R. If you know any other interesting document about panel data in R please let me

hello I have covariance stacionary proces, and i want to estimate some parameter of this proces via gmm. My problem is with write "g" -function. 0 order autocovariance is not problem 1 and higher order autocavariance are problem, because add order from 0 mean that I "loose" one "observacion" if I have 100 observation and i am going to use mean, variance and first