search for: pmaclean2011

Is there a way of executing SAS codes in R environment?   Peter Maclean Department of Economics UDSM [[alternative HTML version deleted]]

...ert tool (part of gnumeric) is very good at converting spreadsheets to csv-files. There is a wrapper in the "gnumeric" package on cran. Cheers, Thomas > Date: Fri, 4 Oct 2013 09:08:50 +0100 > From: Barry Rowlingson <b.rowlingson at lancaster.ac.uk> > To: Peter Maclean <pmaclean2011 at yahoo.com> > Cc: "r-help at r-project.org" <r-help at r-project.org> > Subject: Re: [R] Importing odf file into R > Message-ID: > <CANVKczMtrLD2rw_UWhRoJo5CUP- > 7taUcmiZO+swRRRvKONMUVw at mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1...

...You > > > > Sent from Windows Mail > [[alternative HTML version deleted]] > > > ------------------------------ > > Message: 3 > Date: Sun, 27 Oct 2013 08:32:55 +0000 > From: Barry Rowlingson <b.rowlingson@lancaster.ac.uk> > To: Peter Maclean <pmaclean2011@yahoo.com> > Cc: "r-help@r-project.org" <r-help@r-project.org> > Subject: Re: [R] Code Book from SPSS Data > Message-ID: > < > CANVKczMM9g6VT8o03myenSKwuYUdEBomQfOxMrk8nhdtk7fXiQ@mail.gmail.com> > Content-Type: text/plain; charset=ISO-8859-1 > &...

I am trying to estimate a gamma function using real data and I am getting the following error messages. When I set a lower limit; the error message is "L-BFGS-B needs finite values of fn" ? For other method the error message is: Error in optim(x = c(0.105286666666667, 0.3472275, 2.057625, 0.329675,? : ? non-finite finite-difference value [1] The codes works fine for simulated data

I am estimating a large model by groups. How do you save the results and?returns the associated quantiles? For this example I need a data frame n?? ?xi??????? mu????????beta 1?? 0.1033614? 2.5389580 0.9092611 2? ?0.3401922? 0.5192882 1.5290615 3?? 0.5130798? 0.5668308 1.2105666 I also want to apply gevrlevelPlot() for each "n" or group. ? #Example n <- c(1, 1, 1, 1, 1, 1, 2, 2, 2,

Hi Mackay and anybody (a) Is it possible to select randomly (let say five grids) and plot? (b) Is it possible to plot five nearest grid in one figure? The original question and improved codes: #I am ploting gridded time series data. I would like the actual lat #and lon value appear on the graph-if possible inside the graph as #numbers. If there is also more elegant ways to plot the graphs I

How to substring the following vector (in data frame b) >b a 11 12 1234 1245 124567 126786 145769 such that: a1??? a2 1???? 1 1???? 2 12??? 23 12??? 34 124?? 567 126???787 145???769 ? I tried logical commands with substr() did not work out.

I am estimating Value at Risk using PerfomanceAnalytics package. The?variables are stored in a data frame. I formated the data variables using zoo() and as.xtx() but it is not working. The working example is below. ##########################################################? reguire(zoo) require(PerformanceAnalytics) reguire(xts) ? year<- c(1991-12-30, 1992-12-30, 1993-12-30, 1994-12-30) R1

I have csv data that extend beyond the area I want for an existing map. I want using the boundaries of the polygon shape file as a cookie cutter so that I can overlay the csv data on map without including anything outside the map boundaries and create a dbf file or shapefile of the clipped data .  The reproducible example: ############################################### library(RColorBrewer)

In glm() you can use the summary() function to recover the shape parameter (the reciprocal of the dispersion parameter). How do you recover the scale parameter? Also, in the given example, how I estimate and save the geometric mean of the predicted values? For a simple model you can use fitted() or predicted() functions. I will appreciate any help. ? ? ? #Call required R packages require(plyr)?

How do I move the legend from default position (right and within the plot)?to the "bottomleft" of the plot? ? interaction.plot(YEAR, ID?GROWTH, legend=TRUE, col = 2:7,xlab="Year", ?????? ylim=c(0,2), ylab="Growth",leg.bty = "o") Peter Maclean Department of Economics UDSM

I am looking for study materisl?on how to conduct 'small area estimation' using structural equations models in R for both longtudinal and repeated cross-section data.?Tried google did not work. There are other regression technique my interest is on structural equation models.?

Do anyone now R packages that run/test stochastic dominance with respect to a function (SDRF) and/or stochastic efficiency with respect to a function (SERF)? Peter Maclean Department of Economics UDSM [[alternative HTML version deleted]]

I have data set with binary responses. I would like to conduct polychoric principal component analysis (pPCA). I know there are several packages used in PCA but I could not find one that directly estimate pPCA and graph the individuals and variables maps. I will appreciate any help that expand these reproducible scripts. #How to conduct polychoric principal component analysis pPCA using #either

I need help in indexing grouped data. In this excample (df1 data), the first child had a first immunization at age 2. The second child had the first, second and third immunization at age 5,10, and 12, the third child had first and second immunization at age 4 and 6 and the fourth child had the first immunization at age 2. I have df1 and I need to create df2 with and "ind' variable that

I am ploting gridded time series data. I would like the actual lat and lon value appear on the graph-if possible inside the graph as numbers. If there is also more elegant ways to plot the graphs I will appreciate more suggestions. ################################# library(ggplot2) library(lattice) month <- c("Jan", "Feb", "Mar", "Apr", "May",