similar to: glm specification where response is a 2col matrix

Dear R-List, I?m new in R, so I hope my question is not to primitive, but I haven?t found a solution in the R-help. I?ve got a datatframe with 3 factors, called xf (9 levels), yf (9 levels), zf (3 levels) and one variable (rt): > xf yf zf rt > 1 1 1 67 > 1 1 1 56 > 1 1 1 60 [...] > 1 1 2 58 > 1 1 2 61 [...] > 9 1

I'm having trouble creating a looping variable and i can't see wher ethe problem arises from any hep gratfully appreciated First create a table x<-table(SURVEY$n_0,exposed) > x exposed False True Under 16 24 1 16-19 68 9 20-24 190 37 25-34 555 204 35-44 330 87 45-54 198 65 55-64 67 35 65+

Hi Sherri, There are examples of topographic maps which you have been pointed to, however, I suspect that you want to know where you can obtain topographic data from rather than a canned example. There are quite a few intricacies to the process so I will go through them for you. (1) Topography files can be found in the geomapdata library. You will probably want to use the maps package too (if

Dear all, Further to my previous message I now made a one-line convencience function to export your currently active graphics window/plot to either Word or Powerpoint in Office-native vector-based DrawingML format using either export2ppt(file="plot.pptx") or export2doc(file="plot.docx") : see

Dear Researchers, I wrote two function to fit a circle using noisy data. 1- the fitCircle() is derived from MATLAB code of * zhak Bucher* from the link http://www.mathworks.com/matlabcentral/fileexchange/5557-circle-fit/content/circfit.m 2- the CircleFitByPratt() from MATLAB code of *Nikolai Chernov *from the link

Hi! Consider d <- data.frame(x=1:10,y=5+1:10, yf=rnorm(10,5+1:10)) x y yf 1 1 6 5.268621 2 2 7 8.623896 3 3 8 8.114830 4 4 9 10.125955 5 5 10 9.977261 ... I plot y and yf against x with xyplot(y+yf~x,data=d,col=c('red','green'),pch=c("a","b")) BUT - I would like that the plot of y against x is with type='l' and the

Hi, I am using the profiling feature of LLVM 3.1. Mostly, it works well, but seems to be wrong in a few cases. For example, I have found that the number of times a function executed is reported as less than the number of times that the basic block calling this function executed (the difference was about 25x!). I am copying out the exact functions below. The relevant parts are in bold font:

Hi! I would like to conduct a Cochran`s Q Test in R, but have not found any suitable function. My first idea was: J <- as.table(matrix(c(6,16,2,4),ncol=2, dimnames = list("C" = c("Favorable","Unfavorable"),"Drug A Favorable"=c("B Favorable","B Unfavorable")))) L <- as.table(matrix(c(2,4,6,6),ncol=2, dimnames =

I believe there is a memory leak in isoreg in the current version of R, as I believe the following shows > gc() used (Mb) gc trigger (Mb) max used (Mb) Ncells 120405 3.3 350000 9.4 350000 9.4 Vcells 78639 0.6 786432 6.0 392463 3.0 > for(k in 1:100) { + + y <- runif(10000) + isoreg(x,y) + } > rm(x) > rm(y) > gc() used (Mb) gc

I would like to use lattice graphics to plot multiple functions (or groups or subpopulations) on the same plot region, using different line types "lty" or colors "col" to distinguish the functions (or groups). In traditional graphics, this seems straightforward: First plot all the data using 'type="n"', and subsequently execute a series of "points"

I am running the following nls equation. I tried it with data that excel was fitting and got the error: singular gradient matrix at initial parameter estimates I thought it was due to a low number of points (6), but when I create a dataset, I get the same problem. If I remove the parameter "a," then it can find a solution. Does anyone know what I can do to fit this model?

Hi, I'm trying to generate a distance matrix between sample pairs (example below). I'm not very familiar with the loop command which I expect I will need for this. The example below demosntrates what I'd like to get out of the data - essentially, to calculate the proportion of positions where two samples differ. Any help much appreciated! Also, any notes on how the functions work

I'm trying to specify a model using lmer with a binary response and interaction term, but I get an error I can't parse (see below). Here is some sample data: Subject Concord Age Disc SVC999MX148SU-F yes u int TOU999JU030S1 yes u int TOU999JU030S1 yes u int TOU999JU030S1 yes u int TUT578MX037S2 yes g int COL140MX114S2 yes yf

Full_Name: Petr Klasterecky Version: 1.8.1 OS: Windows XP, Linux Submission from: (NULL) (195.113.27.212) The isoreg() function causes R to crash when called repeatedly. Consider the following simple script: { library(modreg) N <- 10 x <- rnorm(N) print("Original x values:") print(x) for(n in (1:N)){print(y <- isoreg(x[1:n])$yf)} } I am able to run (call) it several

Hi, I am segfaulting when using predict.loess() (checked with r62092). I've traced the source with the help of valgrind (output pasted below) and it appears that this is due to int overflow when allocating an int work array in loess_workspace(): liv = 50 + ((int)pow((double)2, (double)D) + 4) * nvmax + 2 * N; where liv is an (global) int. For D=1 (one x variable), this overflows at

Thank you for your hint. I really mean I am planning to store millions of files on the file system. Then may I ask that what is the maximum number of files which could be stored in one directory without affecting the performance of web server? At 2018-11-03 16:03:56, "Walter H." <Walter.H at mathemainzel.info> wrote: >On 03.11.2018 08:44, yf chu wrote: >> I have

Thank you for your advice. I know the issue depends on a lot of factors. Would you please give me some detail information about how to tune these parameters such as the size of cache,the type of cpu? I am not quite familiar with these detail. At 2018-11-03 22:39:55, "Stephen John Smoogen" <smooge at gmail.com> wrote: >On Sat, 3 Nov 2018 at 04:17, yf chu <cyflhn at

How can I pass an array to a define? It''s not documented in the puppet language guide. I''ve got: define lvm::create_vg ( $pvdisks ) { exec { ''pvcreate'': command => "/sbin/pvcreate -yf $pvdisks", unless => "/sbin/pvdisplay $pvdisks", ... } } class someclass { lvm::create_vg {

hi i got a strange error for openssh-4.0p1 and openssh-4.1p1 (didn't try other versions) with openssl-0.9.7g on AIX 5.1 openssl-0.9.7g and openssh build without errors, "make test" for openssl returns no errors, "make test" for openssh stops at the first connection test "make test" for openssh with openssl-0.9.6m returns no errors (i don't get errors for

Hello Jean-Marc, When you count the adaptation step size in case where the filter has already had minimal adaptation, you use the following expression: r = (0.7*r + 0.3*15*RER*e)/e*(power[i] + 10) = (0.7*leak_estimate*Yf[i] + 0.3*15*RER*(Rf[i] + 1))/((Rf[i] + 1)*(power[i] + 10)). Why do we need this weighted sum and the component 0.3*15*RER*(Rf[i] + 1)? Why use the correlation-based RER if