similar to: Running MCMCpack

Hello All, I'm trying to run a one-dimenional irt model using the packages MCMC and coda on a rather large set of roll-call voting data with many missing observations. Here's a sample of the code: Post10<- MCMCirt1d (Italy10, burnin = 1000, mcmc=50000, thin=100, verbose=TRUE, theta.constraints = list(V549=1, V443=-1)) The MCMCirt1d command seems to work fine, but when I try to

I am trying to install the game "Perfect World" (http://pwi.perfectworld.com/download) On my Mac using Wine. After going through the obnoxious process of installing wine (for me this took forever because I am hugely computer illiterate, but through this process I learned to find info in forums AND use the terminal! Go me!), PW did not work. It installed fine, it would even launch but

Dear all, I would like to compare two distance matrices to find the correlation between those and create a XY plot.could you please kindly help me and let me know how I can do that using R.(Mantel test?). Your help is very much appreciated!! Many thanks! Regards, Sima Fakheran ............................................................... Sima Fakheran Esfahani PhD Student Institute of

Hi everyone, Is it possible to use tapply(x,y,mean) if not all groups of x by y are of the same length (for example if you have one missing observation)? I tried tapply(x,y,mean,na.omit=T) but it doesn't work! Steffi -- --------------------------------- Stefanie von Felten Doktorandin ETH Z?rich Institut f?r Pflanzenwissenschaften ETH Zentrum, LFW A 2 Telefon: 044 632 85 97 Telefax: 044

Dear R-users Problem: Given the following system of ordinary differential euqations dM/dt = (-n)*M-h*M dS/dt = n*M-h*S+u*R dA/dt = h*S-q*A dI/dt = q*A-p*I dJ/dt = h*M-v*J dR/dt=p*I+v*J-u*R where M,S,A,I,J,R are state variables and n,h,u,q,p,v parameters. I'm able to calculate the likelihood value based on the solutions M,S,A,I,J,R of the ODE's given the data, but without an explicit

I'm using an SVM as I've seen a paper that reported extremely good results. I'm not having such luck. I'm also interested in ideas for other approaches to the problem that can also be applied to general problems (no assuming that we're looking for spirals). Here is my code: library(mlbench) library(e1071) raw <- mlbench.spirals(194, 2) spiral <-

I have been installed PWI and I can get into the game but the mouse icon does not show up. I can get in and build a character if I hunt around and click blindly. That is the farthest I have gotten so I don't know if the game is going to work or not. I am fairly new to the Linux system but I have been learning. Does anyone have any ideas as to how to fix this?

Hi Samba folks, I had a couple questions about password complexity checking. To preface, in smb.conf, we set: check password script = /usr/local/sbin/crackcheck -d /usr/share/cracklib/pw_dict Also, if I understand correctly: /usr/local/sbin/crackcheck comes from samba source rpm package. maybe we need to compile it ourselves. /usr/share/cracklib/pw_dict* comes from cracklib-dicts rpm

Hello, I installed PWI, then run and started to update, when finished updating, close it and open it again, and I get the following screen: [Image: http://s2.subirimagenes.com/imagen/previo/thump_6860587screen.png ][/img] see, I get the launcher in black. besides, I get a window with this error: err:wininet:CommitUrlCacheEntryInternal no free entries

Hi, there! I finally became a disciple of 'R', after having lost years of my life handling data with a popular, rather wide-spread spreadsheet-software. Now I want to plot the results of many runs of my simulation over time, so that the means +/- Standard error are on the y-axis, and time on the x-axis. I have tried 'boxplot', with timesteps as the grouping variable, but did

Dear all! After hours of trying around, I gave up: I have a 2-dimensional array, and I know how to split it into its rows and how to get the mean for every row using 'sapply'. But what I want is to calculate the mean over the first n rows, and then the second n rows, etc., so that I get a vector like: v == mean1(row 1:5), mean2(row6:10),... (trivial, you might say. I find it rather

I'm trying to get the game Perfect World International to run under Wine. The installer ran without an apparent problem, but when I try to run the patcher (what usually runs next under windows), I get the following error: almond at Skuld[18]% wine /Users/ralmond/.wine/drive_c/Program\ Files/Perfect\ World\ Entertainment/Perfect\ World\ International/patcher/patcher.exe > pwi.log

Morning, This weekend I attempted an upgrade of my primary samba server from 3.0.24 to 3.3.3. When testing this primary server after the upgrade I had a few issues, so rolled back the upgrade until I can find solutions. This server also has the OpenLDAP server local to and co-located with samba. The two things that initially didn't seem right are that each time I logged into a windows XP box

Hi all, I am trying to find estimates for 7 parameters of a model which should fit real data. I have a function for the negative log likelihood (NLL) of the data. With optim(method="L-BFGS-B",lower=0) I am now minimizing the NLL to find the best fitting parameters. My problem is that the algorithm does not converge for certain data sets. I have read that one should scale the fn

Dear all, is this a bug or do a miss something (this concerns version 0.63.1 on Linux)? > boxplot (c (1, 2, 3, 4, 5), ylim= c (-5, 5)) Warning: parameter "ylim" couldn't be set in high-level plot() function > boxplot (c (1, 2, 3, 4, 5), pars= list (ylim= c (-5, 5))) Warning: parameter "ylim" couldn't be set in high-level plot() function > and actually in

Dear all, this is a rather basic question; i am not sure how to structure my data well: I want to extraxt various measures from my raw-data. These measures are of different sizes, so I decided to store them in a list, like: run1 <- list(Dom = (my_vector), mean = (my_single_number)) I can do that in a for loop for 40 runs, ending up with 40 lists: run1, run2, ..., run40. To have all the

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Dear all, I would like to use the biological symbols for male and female as plotting symbols in a scatterplot (ideally filled and non-filled). R does not seem to have these symbols using pch= in plot() nor are they implemented via expression() or at least I did not find them. I found that the symbols are e.g. available in the wasysym and the marvosym package in LaTeX. I have coded two very rough

Dear all, well, today I just seem to run into these problems ... (I am still working on R 0.63.1 on a SuSE 5.2). I would like to have a data.frame that also contains character (and NOT factor) vectors. Two things that I tried (the second worked in R 0.62.2): > test.df <- data.frame (X= character (4)) > test.df$X [1] <- "a" Warning: invalid factor level, NAs generated or:

> Apart from ipa are there any other good tools out there for centralised user > auth? I am currently testing LDAP (openldap) combined with nss_ldap, configured with authconfig. I would start by testing IPA. Redhat is building out a set of enterprise management tools which include cobbler and spacewalk; I would think IPA will eventually be integrated into their mgt tools. That said, if