similar to: sampling from general multivariate pdf

(apologies if this looks like a re-post, I just sent a similar message to the r-help mail list. This version is via Nabble.) My intended application is error propagation using the ISO GUM Supplement 1 approach (propagation of distributions using Monte Carlo strategies). To automate uncertainty analysis I typically have the following data: (1) a measurement function y(x1,x2,...xn) (2) 'n'

Hello, I'm using dmnorm from the package {mnormt} and getting strange results. First, according to the documentation, dmnorm should return a vector of densities, and I'm only getting one value returned (which is what I would expect). I've been interpreting this as the joint density of all values in the x vector (which is what I want). Should a vector of densities be returned, and if

Dear R-help list members, I have created a package 'mnormt' with facilities for the multivariate normal and t distributions. The core part is simply an interface to Fortran routines by Alan Genz for computing the integral of two densities over rectangular regions, using an adaptive integration method. Other R functions compute densities and generate random numbers. The starting

Going through the examples and reading through the wiki's I still have not found exactly what I was looking for in matching uid's and gid's. using samba samba3-3.0.32-36 We currently have a domain controller Samba/LDAP PDC. samba-3.0.20b-1 Previous member servers samba-3.0.10-1.4 and I went to add a member server. Now I find that users and groups don't match. So from

Hi All. I'd like to generate a sample of n observations from a k dimensional multivariate normal distribution with a random correlation matrix. My solution: The lower (or upper) triangle of the correlation matrix has n.tri=(d/2)(d+1)-d entries. Take a uniform sample of n.tri possible correlations (runi(n.tr,-.99,.99) Populate a triangle of the matrix with the sampled correlations Mirror the

For me getting a member server to work I did not need winbind just ldap was sufficient. Did you made the trust account? Getent group and passwd must give you all users and groups. You must be able to chmod domainuser:domaingroup on your Member Server. What I recognized is that the member server with samba 3.028 is much to slow. It takes too long if you try to connect over My Network Places. This

I have been working the the pmt function in the {mnormt} package and which requires "S a positive definite matrix representing the scale matrix of the distribution, such that S*df/(df-2) is the variance-covariance matrix when df>2; a vector of length 1 is also allowed (in this case, d=1 is set)" is there a way that I can specify the scale covariance matrix instead? Or

[Hope this is the right list where to send...] An attempt to update package 'mnormt' involves the addition of a small new function called 'pd.solve'. When I come to the package checking stage, an error occurs in parsing pd.solve.Rd. The full transcript of the outcome is copied below (it includes details on my installation) but the critical point is where the \examples{} section

Hi all, I I need to manipulate the x axis in a stripchart. I will use one of the data sets included in R to explain what I need to do. attach(ToothGrowth) stripchart(len[supp=='VC']~dose[supp=='VC'], vertical=TRUE, group.names=c('A','A','A')) stripchart(len[supp=='OJ']~dose[supp=='OJ'], add=TRUE, vertical=TRUE, at=c(1:3)+.1,

I'm trying to plot a scatterplot of two variables using pch to plot different characters based on a third factor. Here is my example > data("ToothGrowth") > attach(ToothGrowth) > levels(supp) [1] "OJ" "VC" > plot(len ~ dose,pch=as.numeric(supp)) > legend(locator(1),pch=as.numeric(supp),legend=levels(supp)) The command as.numeric(supp) returns 2 2

Hello, I've problems installing several packages in my R on Fedora 27 64 bit. I found out that it has to do something with a missing compiler (libgfortran.so.3, see below). It works if I downgrade the current version of libgfortran to the specified version by downloading libgfortran-6.2.1-2.fc25.x86_64.rpm and manually installing it. However, I don't want to mess up my system,

This relates to a message from Christophe Pallier to r-help some time ago. Like myself, he finds aggregate very useful, but the interface a little cumbersome. I've implemented a more compact formula interface, found at the bottom of this message: data(ToothGrowth) # I used to aggregate like this: aggregate(list(len=ToothGrowth$len),

Hello, I'm using the eclat function of the arules package (1.0-6) for the identification of frequent itemsets. I need the tidLists, but if I set in the function tidLists=TRUE R crashes (Windows XP Professional SP3, 32 bit, R version 2.12.1 (2010-12-16), reproducible on two different computers) with two different error messages or non at all. Minimum examples are: library(arules)

Kaiya Liu wrote: > > Hi, > > I encountered a problem when I did interaction plots using the > interaction.plot() function. The labels for the trace factor appeared on > the plot when I did the plots in separate windows, but disappeared/cut off > when multiple plots were placed in one window. Can anybody help me? >From the example of ?interaction.plot: data(ToothGrowth)

I am new to Ruby and to Ruby on Rails and have spent the last few days trying to get my head around it. I think I am already hooked but am suffering badly from years of programming a non-OO RDBMS. I am trying to learn by attempting to produce a system for my work. I have a main table sdocs which has a field supplier_id which is an entry in another table (suppliers) to lookup supplier name etc.

I'm having trouble with a simple application with metRology. I need to estimate the uncertainty of the density thickness of seven sheets of film. This is calculated from measurements of mass, length and width of rectangular samples of film. It's not too hard to calculate the whole thing with a little Monte Carlo loop. I get about 0.07 with this: #sample area

Dear Sir I am getting errors in using following R code for L-moment ratio diagram Help in this regard > "plotlmrdia" <- + function(lmr, + nopoints=FALSE, + nolines=FALSE, + nolimits=FALSE, + nogev=FALSE, + noglo=FALSE, + nogpa=FALSE, + nope3=FALSE, + nogno=FALSE, + noexp=FALSE, + nonor=FALSE, +

Thank you for your guys reply for my previous question. But I got one more question about the boxplot. With the code in the R-help: boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset = supp == "VC", col = "yellow", main = "Guinea Pigs' Tooth Growth", xlab = "Vitamin C dose mg", ylab =

Dear List, ? I am creating a boxplot with two subsets, very similar to the example by Roger Bivand at ?boxplot (reproduced below).? I am trying to change the labels on the x-axis to have one number to cover both subsets.? I can do this in other plots by using axis=FALSE followed by a separate axis() command.? I have also tried variations in the names= argument but can't get it to work.?

I've developed a function for updating Flac metadata blocks inside Ogg files. Obviously it is desirable to use padding so that most edits don't require rewriting massive files. Although the flac encoder creates a padding metadata block by default in Ogg Flac files, this appears to be useless. When packets on an Ogg page change size, the page size may change by as different amount, hence