similar to: betabinomial model

Buenos Días, en esta ocasión les escribo para tratar de solucionar una duda que puede resultar muy sencilla, pero de la que no encuentro solucion. necesito hallar la descomposición en valores singulares; tengo la siguiente matriz de datos: DPTO PERSONAS HACINA INTERNET COMPU TV_SUCRIPC Antioquia 4962749,1 1,61 0,39 0,47 0,77 Atlántico 2191360,7 1,97 0,3 0,31 0,42 Bogotá D.C. 7499198 1,67

Hello, I used a mac G5, R.2.1.1, and G77 3.4.4 and I would like to use and call a fortran subroutine. The trouble is that it seems I am not able to correctly load the compiled code. Here is what I have done: In the terminal this how I compiled my fortran code: R CMD SHLIB ~/Desktop/Fortan_kmeans/kmeans3.f There is the wrapper I have paste inside de kmeans3.f file: c

Nothing has been defined <;0;3;(0,23)> failure on __vt_9type_info:G(0,65)=ar(12,3);0;3;(0,23) at ;0;3;(0,23) Nothing has been defined <;0;7;(0,11)>

Okay, here's the debugger output: WineDbg starting... on pid 806e830 No debug information in 32bit DLL 'F:\Mapper6.exe' (0x00400000) No debug information in ELF 'wine' (0x00000000)

hi, a) i have something like: ecdfgrp1<-ecdf(subset(mydata,TMT_GRP==1)$Y); ecdfgrp2<-ecdf(subset(mydata,TMT_GRP==2)$Y); how can i plot the difference between these 2 step functions? i could begin with ecdfrefl<-function(x){ecdfgrp2(x)-ecdfgrp1(x);} ... what next? b) if i have a vector with repeated numeric values how can i get the subset without repeated values .e.g (0,4,0,2,2)

Hi, I need to write two inequalities depend on cumulative distribution (CDF) of beta binomial distribution where alpha and beta are unknown and need to find them. CDF of betabinomial(2,10,alpha,beta) <0.3<=CDF of betabinomial(3,10,alpha,beta) and CDF of betabinomial(5,10,alpha,beta) <0.8<=CDF of betabinomial(6,10,alpha,beta) How I can do that using r studio package? I tried to do

I have a problem, would really appreciate if anyone could help me. I am trying to use bootstrap to resample a dataframe n then use that sample of dataframe to calcuate parameter estimates using betabinomial model. I am kind of lost coz the bootstrap function requires us to give statistic (A function which when applied to data returns a vector containing the statistic(s) of

I find myself wanting to plot three graphs side by side 'as if' they were panels -- that is, with the same y-axis limits, no space between the graphs, and precise vertical alignment of the plot areas. However, I don't want strip titles; I want each graph to have its own x-axis label, on the bottom of the plot. The best way I have so far found to do this is to fake up a data frame that

Hello, i've noticed a new possible missed optimization while testing more trivial code. This time it's not a with a xor but with a multiplication instruction and the example is little bit more involved. C code: typedef short t; t foo(t a, t b) { t a4 = a*b; return a4; } argument "a" is passed in R15:R14, argument "b" in R13:R12, the return value is stored in

Hello, I've got an Asterisk server with a TDM400P with 2 FXO modules and 2 FXS modules. This server is connected to 2 PSTN lines and 2 analog phones. In my Zaptel configuration, I've defined 2 groups : one for the 2 FXO's and one for the 2 FXS. The asterisk server is just used to add a little IVR and Voicemail service. Eveything works fine, but sometimes the conversation is

Hi, all. I want to write some functions like glm() so i studied it. In glm.fit(), it calls a fortran subroutine named "dqrfit" to compute least squares solutions to the system x * b = y To learn how "dqrfit" works, I just follow how glm() calls "dqrfit" by my own example, my codes are given below: > qr <- >

I am creating habitat selection models for caribou and other species with data collected from GPS collars. In my current situation the radio-collars recorded the locations of 30 caribou every 6 hours. I am then comparing resources used at caribou locations to random locations using logistic regression (standard habitat analysis). The data is therefore highly autocorrelated and this causes Type

Hello! I have a curious problem, which I cannot solve. With my code I solve thousands of small linear programs with the package lpSolve automatically. But R crashes sometimes (~always, but always on different linear programs) in a strange way. For illustration, I tried to prepare a simple example, which shows the nature of the problem. The function aaa (see below) declares some constants (only in

Dear useRs, upon request, the VGAM package (currently version 0.7-1) has been officially released on CRAN (the package has been at my website http://www.stat.auckland.ac.nz/~yee/VGAM for a number of years now). VGAM implements a general framework for several classes of regression models using iteratively reweighted least squares (IRLS). The key ideas are Fisher scoring, generalized linear and

Hi all, I just installed centos 5.4 xen-kernel on intel core i5 machine as dom0. After some hours of syncing a raid10 array (8 sata disk) I noticed a steadily increasing loadavg. I think without reasonable i/o wait or cpu utilization the loadavg on this system should be very lower. If this loadavg is normal I would be greatful if somone could explain why. The screenshots below show that there is

On Sep 4, 2010, at 11:21 AM, Borja Ferrer wrote: > I've noticed this pattern happening with other operators aswell, but used xor in this example. As i said before, i tried with different register allocation orders, but it will produce always the same result. GCC is emitting longer code, but since LLVM is so nearer to the optimal code sequence i wanted to reach it. In LLVM, copies are

Dear all, I just installed R 2.5.0alpha and noticed new warnings related to the extractor "$" when using contributed packages. For instance: > library(RODBC) Warning message: $ operator is not valid for atomic vectors, returning NULL > library(aod) Package aod, version 1.1-18 > data(orob2) > m1 <- betabin(cbind(y, n-y) ~ 1, random = ~ 1, data = orob2) >

Hello, Thank you for your response. I need to install RankTail package since it contains the beta binomial distribution, CDF and inverse CDF in the usual form which I need to use. However rmutil package contain unusual forms for these functions. So it is easier for me to deal with the forms are contained in RankTail. I tried to install bioconductor package, using the following commands but I

Hi, I did a quick search for other packages that provide the beta binomial distribution and found "rmutil". > install.packages("rmutil") The package has the CDF (pbetabinom) and inverse CDF (qbetabinom) among other functions. HTH, Eric On Wed, Nov 1, 2017 at 7:50 AM, MCGUIRE, Rhydwyn < rmcgu at doh.health.nsw.gov.au> wrote: > Hi there, > > It looks like

Hi list, I’m trying to plot lines in the same graphic. There are lots of similar topic on the list history. I’ve read about xyplot (lattice) and other ways to do that but I only got complex error messages, and wasn’t able even to start a decent R scrip. Well I have some paired information (each H W ), were H is the y axis and W is the x axis, how could I plot all the four lines at the same image?