similar to: extracting p-values from the Manova function (car library)

hi, i am getting an error when calling the impute.knn function (see the screenshot below). what is the problem here and how can it be solved? screenshot: ################## *** caught segfault *** address 0x513c7b84, cause 'memory not mapped' Traceback: 1: .Fortran("knnimp", x, ximp = x, p, n, imiss = imiss, irmiss, as.integer(k), double(p), double(n), integer(p),

Hi , is it possible to calculate ld-measures D, D', r and perhaps corresponding p-values with r IF THE PHASE IS KNOWN? The genetics - package provides the LD function only for ambigious phase. Thank you very much Bettina Kulle

Dear useRs, the new R package exams provides Sweave-based automatic generation of exams with multiple-choice questions and arithmetic problems. The package is available from CRAN: http://CRAN.R-project.org/package=exams It includes a vignette giving an overview of the main design aims and principles as well as strategies for adaptation and extension. Hands-on illustrations - based on example

Dear useRs, the new R package exams provides Sweave-based automatic generation of exams with multiple-choice questions and arithmetic problems. The package is available from CRAN: http://CRAN.R-project.org/package=exams It includes a vignette giving an overview of the main design aims and principles as well as strategies for adaptation and extension. Hands-on illustrations - based on example

Dear list: I would like to recreate how the artificial data set BregFix was generated in package flexmix (thanks Bettina and Friedrich). The data set is thoroughly described in Grun and Leisch's Computational Statistics & Data Analysis 51(11) :5247-5252 but references to the appropriate seed number(s) are missing (providing these details was certainly beyond the scope of the authors'

Hello, I have a dataset with proportions that vary around a fixed mean, is it possible to use betareg to look at variance in the dispersion parameter while keeping the mean fixed? I am very new to R but have tried the following: svec<-c(qlogis(mean(data1$scaled)),0,0,0) f<-betareg(scaled~-1 | expt_label + grouped_hpi, data=data1, link.phi="log",

Dear useRs, we are organizing the following session Topic: Innovative Tools in Data Analysis Organizers: Achim Zeileis and Bettina Gruen at the First Workshop of the ERCIM Working Group on Computing & Statistics June 19-21, 2008 Neuchatel, Switzerland URL: http://www.dcs.bbk.ac.uk/ercim08 To improve the quality of statistical data analysis the provision of innovative tools which make new

I am using 'tm' package for text mining and facing an issue with finding the frequently occuring terms. From the definition it appears that findFreqTerms and minDocFreq are equivalent commands and both tries to identify the documents with terms appearing more than a specified threshold. However, I am getting drastically different results with both. I have given the results from both the

Dear colleagues, I had a question wrt the car package. How do I evaluate whether a simpler multivariate regression model is adequate? For instance, I do the following: ami <- read.table(file = "http://www.public.iastate.edu/~maitra/stat501/datasets/amitriptyline.dat", col.names=c("TCAD", "drug", "gender", "antidepressant","PR",

Problem on flexmix when trying to apply signature developed in one model to a new sample. Dear R Users, R Core Team, I have a problem when trying to know the classification of the tested cases using two variables with the function of flexmix: After importing the database and creating a matrix: BM<-cbind(Data$var1,Data$var2) I see that the best model has 2 groups and use: ex2

Hello, I would like to annotate a graph with the expression 'alpha = 5%' (the alpha should be displayed as the greek letter). I tried > text(1,1,expression(alpha == 5%)) which gives a syntax error. escaping the percent sign (\%) or doubling (%%) does not help. What do I do? Thanks, Martin Keller-Ressel -- Martin Keller-Ressel Research Unit of Financial and Actuarial

Hallo! I want to use the package bayesmix. Trying the examples I had no success. The reason is, I think: > haveJAGS() [1] FALSE Have read of JAGS executable (I could not find JAGS.exe in my R directory). What is JAGS? Where can I find or download it? Karl

On Friday 05 October 2007 03:54, puppet wrote: > #829: Multiple source with recursion error > -----------------------------+--------------------------------------------- >- Reporter: paul | Owner: michael > Type: defect | Status: closed > Priority: normal | Milestone: unplanned > Component: library |

Hi, How do i convert from Maildir to dbox using the convert plugin? My users reside in /home/user/Maildir. In the dovecot.conf file i've set up convert_mail = dbox:%h but I am not sure about the syntax. Im getting the error: Couldn't create mailbox directory INBOX Mailbox conversion failed, exiting Thanks -- Kjartan

I've just seen this post: http://www.mail-archive.com/asterisk-users@lists.digium.com/msg41132.html and it took me back to play again with my dust collecting 2000W. Does anybody got DTMF to work? My sip.conf looks like this: [400] type=friend context=from-sip username=400 secret=verysecret disallow=all allow=g729 dtmfmode=rfc2833 host=dynamic nat=yes qualify=300 canreinvite=no My phone is

I am trying to cluster some binary data using k-means . As the regular "kmeans" available from stats package in R does'nt provide the option to change the distance method. I was wondering there is any package available to specify type of distance measure to be used in k means clustering in R. Especially distances like "Jaccard" which is good for binary data.

Dear R-helpers & Friedrich Leisch, I'm trying to fit a mixture of regression models on attached data set by doing the following: agl <- read.table("agl.txt") mod1 <- flexmix(resp~gng+csc|Subject,data=agl, model=FLXglm(family="binomial"),k=2) The result is a (varying) number of the following warnings: Warning messages: 1: non-integer #successes in a

En innebygd og tegnsett-uspesifisert tekst ble skilt ut... Navn: ikke tilgjengelig Nettadresse: https://stat.ethz.ch/pipermail/r-help/attachments/20080201/6cd3960c/attachment.pl

Hello all, I am trying to fit a truncated mixture model and I wrote a driver for flexmix following the example in the vignette, but it doesn't work for me: it assigns all data points to one component only, e.g.: > > source('bugged.R') > > Call: > flexmix(formula = x ~ 1, k = 2, model = truncatedmodel(lower = -4, > upper = 4)) > > prior size

Hi everyone I have got a quick question: I the "seqinr" package: *dist.alignment(x,"identity")* This is calculating the square root of pairwise distances. Does anyone know whether/how gaps are counted in this function? Thank you. Best wishes, Bettina [[alternative HTML version deleted]]