similar to: My First Function: cryptic error message

I wonder if someone could help me adjusting the following code to parallelized snow code: #Creating a data set (not needed to be parallel) n<-100 p<-100 x<-matrix(rnorm(n*p),p) y<-rnorm(n) # Bootstrapping nboot<-1000 alpha<-0.05 rhoboot <- array(0, dim=c(p,nboot)) bootranks <- array(0, dim=c(p,nboot)) bootsamples <- array( floor(runif(n*nboot)*n+1), dim=c(n,nboot)) for

Dear R listers, I have a data file looks like the following: Testing marker: s_1 --------------------------------------------- Allele df(0) -LnLk(0) df(T) -LnLk(T) ChiSq p 3 7995 29320.30 7994 29311.85 16.90 4e-05 (2229/8000 probands) Testing marker: s_2 --------------------------------------------- Allele df(0)

I am involved in a study where, as in most of life, men demonstrate themselves to be recalcitrant. So while we have many probands and most of their mothers we only have about 50% of the trios being complete. I have been running tdt and trio.types. It appears as if it is ignoring the duos. Sometimes a duo can be informative. For instance Father ..missing Mother 1/2 Proband 1/1 This duo shows that

Hi, can anyone help me fit betabinomial model to the following dataset where each iD is a cluster in itself , if i use package aod 's betabinom model it gives an estimate of zero to phi(the correlation coeficient ) and if i fix it to the anova type estimate obtained from icc( in package aod) then it says system is exactly singular. And when i try to fit my loglikelihood by

In October of last year Marc Blanchet of the Canadian company Viag?nie made a presentation on how he and others had build IPv6 support into Asterisk and furthermore demonstrated that it worked. Marc Blanchet went into some details on how it was done and the amount of work that had gone into it. A question is this connection: When will these additions make their way into the Asterisk mainstream

Greetings all The help file for GAMM in mgcv indicates that the log likelihood for a GAMM reported using summary(my.gamm$lme) (as an example) is not correct. However, in a past R-help post (included below), there is some indication that the likelihood ratio test in anova.lme(mygamm$lme, mygamm1$lme) is valid. How can I tell if anova.lme results are meaningful (are AIC, BIC, and logLik

I tried to reproduce a result from a former colleague which he got with S-plus bootstrap method. I don't have S-plus at hand. In R, there are 2 packages related to bootstrap method, bootstrap and boot. The former has a function called 'bootstrap' but this does not seem to conform either to the function used in S-plus nor to that described in MASS, 3d ed., p.144. The latter seems to be

Bottom Line Up Front: How does one reshape genetic data from long to wide? I currently have a lot of data. About 180 individuals (some probands/patients, some parents, rare siblings) and SNP data from 6000 loci on each. The standard formats seem to be something along the lines of Famid, pid, fatid, motid, affected, sex, locus1Allele1, locus1Allele2, locus2Allele1, locus2Allele2, etc In other

Hello, i only got a small problem. i try to create automatic new dataframes, or png?s. the main problem i got is: how can i create automatic a new name for a file (read out by simply "for") - i tried to use "(paste...) but theres an errormessage, about a wrong declination. R told it is as.character, but need as.Real. Should i use another method than "paste"? i tried as

Suppose I have a Master.Rnw file that looks something like this: \documentclass[12pt]{mypaper} \usepackage[margin=1in]{geometry} \usepackage{setspace} \usepackage{url} \usepackage{indentfirst} \usepackage{fancyhdr} \usepackage{Sweave} \pagestyle{fancy} \lhead{sonographic rectal diameter and ADHD} \rhead{ } \usepackage{abbrevs} %\usepackage{natbib} %\usepackage{apacite}

Thanks for the advice. My question is more on how to do this? Let me use a biology gene analysis example to illustrate: In biology, there are always some house keeping genes which differ little even at pathological conditions. We know that at different batches, there are external factors affect the measurements. For example, overall signal intensity might be different due to lab reagents. A

Subject says it all... A Simple Asterisk Based Toll Fraud Prevention Script http://www.infiltrated.net/asterisk-ips.html Ramblings for admins/engineers to think about. Doesn't have to cost you umteen thousand dollars for stuff like IPS/IDS. Although a little on the crude side, quite effective. If you care to dabble with MySQL you can create quite an impressive hosts based IPS that is custom

Have you seen how these sub-one-cent companies take off on Good News? For the past two months, every one we have brought you has shown Amazing Appreciation. At such a low price even the smallest Gain means a Significant Percentage Return. We called (P)(P)(T)(L) as one to watch on Friday because of a Highly Anticipated Report from the Field. It moved up 13% on Friday and the news

Have you seen how these sub-one-cent companies take off on Good News? For the past two months, every one we have brought you has shown Amazing Appreciation. At such a low price even the smallest Gain means a Significant Percentage Return. We called (P)(P)(T)(L) as one to watch on Friday because of a Highly Anticipated Report from the Field. It moved up 13% on Friday and the news

I'm interested in further understanding the differences in using many classification trees to improve classification rates. I'm also interested in finding out what I can do in R and which methods will allow prediction. Can anybody point me to a citation or discussion? Specifically, I want to classify remotely sensed imagery where training data is extracted on class membership by the user.

I'm creating my own package for personal and I'm having trouble getting it to a point where R (v 2.3.1) will recognise it. I've followed two different tutorials for how to create the package structure and the DESCRIPTION file ( http://web.maths.unsw.edu.au/~wand/webcpdg/rpack.html , http://www.maths.bris.ac.uk/~maman/computerstuff/Rhelp/Rpackages.html#Lin-Lin ). I'm still getting

I am getting an error that I don't understand, and wonder if anyone could explain what's going on. I call a function defined thus: clogit.rds<-function(formula,data,extra.data,response.prob, na.action=getOption("na.action"),subset=NULL, control=coxph.control()){ method="exact" # only option for now mf<-match.call()

Hi, This problem may be related to a configuration problem but I believe it is a bug in the FOP since restarting the FOP server clears the problem. Here is the scenario: Using AgentCallBackLogin and have four agents logged in a call is made to one of the agents directly from an internal phone. Okay so far. Call is hung up and the same extension is used to call another agent okay again, no

Hi I have a 3 x 2 contingency table: 10 20 30 40 50 60 I want to update the frequencies to new marginal totals: 100 130 40 80 110 I want to use the ipf (iterative proportional fitting) function which is apparently in the cat package. Can somebody please advice me how to input this data and invoke ipf in R to obtain an updated contingency table? Thanks. By the way I am quite new to R. -- Dr

Hi, Short version: Is it possible to enable *all* builtins architectures? How? Longer explanation: I don't quite understand how compiler-rt targets are decided to be built while invoking cmake from the LLVM tree. I tend to believe that builtins and sanitizers (the two big parts in compiler-rt) are quite different monsters, but however they share the same build decisions and (if I'm not