similar to: PCM() non-convergence problem: how to get diagnostics?

Hello,everybody! I am a beginner of R. And I want to ask a question. If anybody would help me, thank you very much ahead. I want to plot something like a response surface, and I find the "rsm" package. Some commands are like this: #code head library(rsm) CR = coded.data(ChemReact, x1 ~ Time, x2 ~ Temp) CR.rsm = rsm(Yield ~ SO(x1,x2), data = CR) summary(CR.rsm) contour(CR.rsm,x1~x2)

Hola compañeros de la lista. Los molesto con la siguiente duda. En un diseño central compuesto (CCD) con dos factores (V1 y V2) y una variable de respuesta (R), utilizando valores codificados (-1.4142, -1, 0, 1, 1.4182), al aplicar la orden: rsm.segundo.orden <- rsm(R ~ Bloque + SO(V1, V2), data = DATOS.Codificados) Obtengo el siguiente modelo: R = 103.92 -2.16

This patch implements a checker to detect instructions which should be paravirtualized instead of direct writing raw instruction. This patch does rough check so that it doesn't fully cover all cases, but it can detects most cases of paravirtualization breakage of hand written assembly codes. Signed-off-by: Isaku Yamahata <yamahata at valinux.co.jp> Cc: "Luck, Tony"

This patch implements a checker to detect instructions which should be paravirtualized instead of direct writing raw instruction. This patch does rough check so that it doesn't fully cover all cases, but it can detects most cases of paravirtualization breakage of hand written assembly codes. Signed-off-by: Isaku Yamahata <yamahata at valinux.co.jp> Cc: "Luck, Tony"

Dear R-Listers, to measure the psychometric curve of pitch discrimination, one sequentially presents two tones of slightly different pitch to an observer (animal will do), and asks "which is higher". The pschometric curve is the fraction of correct responses plotted against the pitch difference. It passes through 50% (pure guessing) at zero and normally approaches 100% at large

Dear useRs, new and updated features in eRm 0.9-6 for extended Rasch modeling: - infit and outfit mean-square statistics added to functions personfit() and itemfit(). - new method: plotPImap() plots a person-item map (cf. Bond & Fox, 2007), i.e., a map of locations of item (and threshold) parameters, and the distribution of person parameters. - new options in plotGOF(): 'conf':

Dear useRs, new and updated features in eRm 0.9-6 for extended Rasch modeling: - infit and outfit mean-square statistics added to functions personfit() and itemfit(). - new method: plotPImap() plots a person-item map (cf. Bond & Fox, 2007), i.e., a map of locations of item (and threshold) parameters, and the distribution of person parameters. - new options in plotGOF(): 'conf':

Dear, I'm trying to run diagnostics on MCMC analysis (fitting a log-linear model to rates data). I'm getting an error message when trying Gelman-Rubin shrink factor plot: >gelman.plot(out) Error in chol.default(W) : the leading minor of order 2 is not positive definite I take it that somewhere, somehow a matrix is singular, but how can that be remedied? My code: library(rjags)

I'm using the fwdmsa package to identify deviant cases in a Mokken scale analysis. I've run into a problem. When I use scales comprising a few items, iI tend to get an error: Error in y[order(res[-msamp])][1:(length(samp) + 1 - length(msamp))] : only 0's may be mixed with negative subscripts I understand that the error is triggered when the algorithm is fetching cases to enter

Hi! I installed XEN 4.3.2 on my Gentoo Linux machine with EFI-Support. DETAILS: Motherboard: ASUS M5A97 Graphiccard: Radeon HD 5000/6000/7350/8350 Series ISSUE: When booting from xen.efi i get instantly a blank screen, but the XEN Hypervisor boots and i am able to SSH into it. Also i am able to issue commands in "blind mode". When i boot from grubx64.efi grub tells me "no

I have an analogue trunk to an AT&T Definity. It has a DISA context defined. From a Definity handset call the analogue port extension 1008 and wait for dial tone from asterisk. It takes between 3&4 rings. Likewise from Asterisk SIP handset <PBX Access No><PBX Extn> takes nearly 10 secs to ring. Is this configurable? Ian Cowley -----Original Message----- From:

Hola compañeros de la lista. Estoy aprendiendo a usar el paquete "rsm" para superficies de respuesta. Siguiendo este ejemplo todo va bien, hasta que trato de obtener el valor máximo predicho para la variable de respuesta en los valores de x1 y x2 estimados. ------------------------------------------------------------- library("rsm") ChemReact CR <- coded.data(ChemReact,

I've been trying several times now: # mount -t smbfs -o credentials=/etc/samba/private/credentials,codepage=cp850,iocharset=iso8859-1,noauto //machine/share\$ /mnt/share Causes: # ls /mnt/share/directory Sp?rsm?l.txt While this: # mount -t smbfs -o credentials=/etc/samba/private/credentials,codepage=cp850,iocharset=iso8859-1 //machine/share\$ /mnt/share Causes: # ls /mnt/share/directory

I would appreciate confirmation that the function vcov(model.nlme) gives the var-cov matrix of the fixed effects in an nlme model. Presumably the random-effects var-cov matrix is given by cov(ranef (model.nlme)? Rob Forsyth

In some trial simulation work I need to create batch files that will repeatedly generate pseudoreplicate datasets and then create non- linear mixed effects models using nlme. Inevitably these models sometimes fail to converge but I need the batch file to simply move on to another simulation rather than abort. I am using the try() function as in model<-try((nlme(...))) which handles

estimados, ajusté una superficie de respuesta con el paquete "rsm" de R y hago un grafico 2D de la superficie de respuesta con contour(), pero como el argumento de esta función es la superficie de respuesta ajustada y ésta se ha ajustado codificando las variables, de acuerdo a la recomendación de la documentación del paquete rsm, obtengo un gráfico con los ejes en unidades codificadas

Over the last couple of years I have written quite a few R programs for various "psychometric" techniques, and I am regularly updating and expanding what is there. I now have (wholly or partially), or have planned -- gifi package (update to homals on CRAN). Code for multiple correspondence analysis, nonlinear principal component analysis, nonlinear multiset canonical

On ia64 smp_rmb, smp_wmb, read_barrier_depends, smp_read_barrier_depends and smp_store_mb() match the asm-generic variants exactly. Drop the local definitions and pull in asm-generic/barrier.h instead. This is in preparation to refactoring this code area. Signed-off-by: Michael S. Tsirkin <mst at redhat.com> --- arch/ia64/include/asm/barrier.h | 10 ++-------- 1 file changed, 2

Hi all (and especially Frank), I'm trying to use x=T, y=T in order to run a validated stepwise cox regression in rsm, having multiply imputed using mice. I'm coding model.max<-fit.mult.impute(baseform,cph,miced2,dated.sexrisk2,x=T,y=T) baseform is baseform<-Surv(si.age,si=="Yes")~ peer.press + copy.press + excited + worried + intimate.friend + am.pill.times +

I have setup a Zentyal 3.2 Samba4 AD. I am using RSAT to administer the AD server. Windows7 and Linux (CentOS) clients are able to join the AD; uid/gid of directories/files show up as per the definitions in Samba4. With Mac OS X clients, there is a problem. The OS X nodes are able to join the domain (quite a few blogs on this subject). The AD connection shows green button. The AD user is