similar to: Simulate Correlated data from complex sample

I am attempting to use optim to solve a neural network problem. I would like to optimize coefficients that are currently stored in a matrix Y=270 x 1 X= 27- x 14 b1= 10x14 b2= 11x1 V= 10 x 14 set of prior variances. I have the following function: posterior.mode1=function(y,X,b_0,b2,V) { log.like=function(b1) { a_g=compute(b1) z_g=tanh(a_g); z_g=cbind(1,z_g)

Hi, I want to find a value of n1. I used the following code but I am getting the error - Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' n=10 a_g<-(1/(n*(n-1)))*((pi/3)*(n+1)+(2*sqrt(3)*(n-2))-4*n+6) a_s<-function(n1) { t1=(n1-1)/2; (t1*(gamma(t1)/gamma(n1/2))^2)/2-1-a_g } xm<-solve(a_s) Can anyone help me out. Thanks in

Dear List: We are generating data such that students are clustered in schools for some item response data for a simulation study. One component of our simulation is to generate measurement error from a logistic distribution with a mean of 0 and standard deviation of 1.7 to match the logistic curve of the Rasch model. We are generating an error term for each of the 40 hypothetical test items a

Hi all, im triyng to implement a bayesian model with R and c++. I have a strange problem. I can't reproduce the error with a small script and then i post the original one. The problem is after the line for(MCMC_iter2=0;MCMC_iter2<thin;MCMC_iter2++) For the first 34 iterations all work fine, after, all the simulations of mu_acc_P return an "nan". If i delete the line

Thanks for replying so quickly Guus. I will try increasing the re-keying time to 86400, and am now running the servers in debug 5 mode logging to a file, so will get all output. Its interesting what you say about UDP being given lower priority, the two end points of this VPN are actually within the same data center, but are being routed within the DC's core routers (because its going

Hi, I am using tinc 1.0.23 in router mode and am having trouble with intermittent TCP connect timeouts that do not occur when routing the TCP over the internet directly. I am using a pure IPv6 setup (i.e direct IPv6 connections, and tunneled IPv6 subnets). My configs are as follows: Site 1: ----------------------- tinc.conf: Name = rsukmhb AddressFamily = ipv6 BindToAddress =

----- "Akemi Yagi" <amyagi at gmail.com> escreveu: > De: "Akemi Yagi" <amyagi at gmail.com> > Para: "CentOS mailing list" <centos at centos.org> > Enviadas: Sexta-feira, 25 de Novembro de 2011 17:56:21 (GMT-0300) Auto-Detected > Assunto: Re: [CentOS] CentOS 6.1- RPS/RFS kernel patch > > On Fri, Nov 25, 2011 at 11:49 AM, Karanbir

On Mon, 2014-01-13 at 11:40 +0200, Michael S. Tsirkin wrote: > On Sun, Jan 12, 2014 at 03:32:28PM -0800, Michael Dalton wrote: [...] > > The last issue is how the rx_queue_attribute 'show' function > > implementation for mergeable_rx_buffer_size will access the appropriate > > per-receive queue EWMA data. The arguments to the show function will be > > the

On Mon, 2014-01-13 at 11:40 +0200, Michael S. Tsirkin wrote: > On Sun, Jan 12, 2014 at 03:32:28PM -0800, Michael Dalton wrote: [...] > > The last issue is how the rx_queue_attribute 'show' function > > implementation for mergeable_rx_buffer_size will access the appropriate > > per-receive queue EWMA data. The arguments to the show function will be > > the

On Mon, Jan 13, 2014 at 7:38 AM, Ben Hutchings <bhutchings at solarflare.com> wrote: > I don't think RPS should own this structure. It's just that there are > currently no per-RX-queue attributes other than those defined by RPS. Agreed, there is useful attribute-independent functionality already built around netdev_rx_queue - e.g., dynamically resizing the rx queue kobjs as

Dear List: I am working to understand some differences between the results of the svymean() function in the survey package and from code I have written myself. The results from svymean() also agree with results I get from SAS proc surveymeans, so, this suggests I am misunderstanding something. I am never comfortable with "I did what the software" does mentality, so I am working to

I have a share created on my linux server. I can connect to it with the "administrator" user just fine. I need to be able to create files and directories with the administrator user, even in directories that are not owned by administrator. administrator is in the users group. Everything under the share is in the users group, and the group has write permission. So why can't

Dear List members I am using the regsubsets function to select a few predictor variables using Mallow's Cp: > sel.proc.regsub.full <- regsubsets(CO2 ~ v + log(v) + v.max + sd.v + tad + no.stops.km + av.stop.T + a + sd.a + a.max + d + sd.d + d.max + RPA + P + perc.stop.T + perc.a.T + perc.d.T + RPS + RPSS + sd.P.acc + P.dec + da.acc.1 + RMSACC + RDI + RPSI + P.acc + cov.v + cov.a +

Hello everyone, I recently bought a Riello IPG 600 UPS ( http://www.riello-ups.com/products/1-ups/39-iplug) to prevent hardware damage to my server in case of power outages (already lost two hard drives because of this) and I managed to successfully configure it with NUT version 2.7.4 available for Ubuntu 17.04: nut/zesty-updates,zesty-updates,now 2.7.4-5ubuntu2.1 all [installed]

Dear all, I have used R.rps to produce a Google API chart (googleVis) with a data request in another server. But i don't understand how is possible to scheduling a request data to the server and after produce a update of the charts. Thanks in advance. Alessandro Oggioni

http://bugzilla.mindrot.org/show_bug.cgi?id=42 hampus.lind at rps.police.se changed: What |Removed |Added ---------------------------------------------------------------------------- Status|NEW |RESOLVED Resolution| |FIXED ------- Additional Comments From hampus.lind at rps.police.se

Hello, I'm setting up nut on a Riello IDG800 ups, using riello_usb driver. I backported nut 2.7.1-1 debian packages from sid to wheezy, and installed nut-client, nut-server and nut-cgi. All seems to work quite nicely, but I want to submit the output of upsc, because some values are not correct (battery.charge, battery.runtime, ups.temperature, for instance): battery.capacity: 7

Hello. Thanx for your answer. It work ! Riello Idialog present in list as IDL 600 But I have 2 bugs battery.charge: 255 ups.temperature: 255 battery.capacity: 7 battery.charge: 255 battery.runtime: 3932100 battery.voltage: 13.7 battery.voltage.nominal: 12 device.mfr: RPS S.p.a. device.model: UNV1 device.serial: device.type: ups driver.name: riello_usb driver.parameter.pollinterval: 2

Hello Guys, I'm using the package RankProduct to determine the differential expression of genes in 2 classes. The importance of the thing is to know: - How can I determine which class is the class of up-regulated genes and which to the down regulated gene - How can i interpret the results : $pfp, $RPs, $RPrank, Orirank,

Hello guys, Today we faced an unique problem. We have about 120 stations and we are reaching 40k simultaneous listeners monday/friday (~ 11 AM GMT -3). We also use url auth to log, audit and build audience reports. But today, at the end of the day we were about 10k simultaneous and listeners start get errors (403, busy, please try again later), and after digging, delete logs and restart