similar to: scale,centre,and get more interactions

Hi, I am playing with clogit(case~spontaneous+induced+strata(stratum),data=infert) from clogit help file. This line works. 1. But, why strata(stratum) doesn't have a coefficient like spontaneous and induced? 2. When I remove strata(stratum) from the command, this function seems to keep running forever. Why? 3. I think the equation for clogit looks like P=1/(1+

Dear list, I'm using lattice::levelplot to plot a coloured image of 3D data. The range of the z values goes from negative to positive, but is not exactly centred around 0. I would however like to map a diverging colour scale with white falling exactly at 0, and both extremes being symmetrical in the legend to better contrast the opposite change in colour saturation. The following dummy

I am trying to create a plot title in R with substitution by a numeric variable (Figure number N) within the text which is bold and has a subcripted part as well. Here is what I have: title <- expression(bold(paste("Figure ", N, ": Plot ", C[max], " versus CrCL"))) plot(1, main="") # Simple plot for testing N <- 5 mtext(title, line=3, font=2,

Hello List I'm moving this over from the bioC list as, although the problem I'm working on is biological, the current bottle neck is my poor understanding of R. I wonder if someone would help me with the following function. cumulMetric <- function(deMirPresGenes, deMirs){ ??? #need to match position of each miR in deMirPresGenes with its FC to form a vector of FC in correct order ?

Dear R-help, I am using R-3.3.2 on Windows 10. I teach on a course which has 4 computer practical sessions related to the development and validation of clinical prediction models. These are currently written for Stata and I am in the process of writing them for use in R too (as I far prefer R to Stata!) I notice that predictions made from a Cox model in Stata are based on un-centred variables,

Hi, I am trying to run a conditional logistic model on a nested case-control study using cph() and then estimate survival based on the model. The data came from Prof Bryan Langholz website where he also has the SAS code to this, so I am trying to replicate the SAS results. The data attached. Basically, the variables are: rstime: risk set age rsentry: fake entry time, just before rstime setno:

Hello, I want to ask how to get call duration. -- Necati DEM?R http://demir.web.tr http://friendfeed.com/ndemir ndemir ~ demir.web.tr --------------------------------------- -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20100603/d4085dc6/attachment.htm

I would like to centre titles for pairs of plots in a 3-x-2 array. Each row of the array corresponds to a calendar year and I would like to have the year value centred between the two plots in the row, and just above their upper edges. I have attached an example in "demo.pdf" showing roughly what I want. I managed to produce the example using

Dear List members I am trying to find a way to generate a random point in a circle centred in a geographical location. So far I have used the following formula (see code below): random_x = original_x + radius*cos(angle) random_y = original_y + radius*sin(angle) where radius is a random number between 0 and the radius of the circle and angle is between 0 and 360 degrees The code bellow works

Hi, Is the following function built in somewhere? concat = function(v) { res = "" for (i in 1:length(v)) res = paste(res,v[i],sep="") res } e.g. > concat(c("12","3","45")) [1] "12345" Cheers, Andy -- Andy Fugard, Postgraduate Research Student Psychology (Room F3), The University of Edinburgh, 7 George

Hi there, Am I correct to believe that tetrachoric correlation is a special case of polychoric correlation when there are only two levels to the ordered factor? Thus it should be okay to use hetcor from the polycor package to build a matrix of correlations for binary variables? If this is true, how can one estimate 95% confidence intervals for the correlations? My guess would be mat =

Hello Sirs, > > Have you ran them through mp3check and mp3_check? > >Yep. Found a couple with a few bad frames and deleted them. Other than >that, nothing worse than a missing ID3 tag or two. Could you please give me the address of these tool(mp3check, mp3_check). I have many bad coded mp3, When Icecast streaming these bad! mp3s it kicks the connected clients. I want to delete

Hi Christian, Chuck (and lists) It seems that the problem may be the strange behaviour of 'unstack' inside a function. See this thread in the R mailing list: http://tolstoy.newcastle.edu.au/R/help/04/03/1160.html Anyway, I got round the problem by using 'aggregate' instead of converting to a list and then tapply to sum values of metric. Probably more efficient as well. Thanks

Dear R Users, I am trying to plot a barchart with a line graph superimposed (using par(new=TRUE)). There are 12 bars and 12 corresponding points for the line graph. This is fine, except that I'm encountering two problems: 1) The position of the points (of the line graph) are not centred on the middle (horizontally) of each corresponding bar. In fact, whilst the first point is located on the

I'm interested in computing densities of a number of small ego-centred or local networks. I've installed both the R program and the sna package. My next task is importing the data into R. Our data was entered into Spss in an edgelist format in which each case specifies the link between an alter and an ego. For example: 10101 10102 1 10101 10103 3 10102 10104 2 10201 10203

This is a really rookie question: when should i use TE110P ISDN PRI Card? -- Necati DEM?R --------------------------------------- -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20100621/07d2917f/attachment.htm

Hello all, There's a new mailing list for researchers in psychology who are learning and using R. New users of R are especially welcome. To join, venture to http://www.jiscmail.ac.uk/archives/psych-r.html and click "Join or leave the list (or change settings)" You can also peruse the archives from that page. Best wishes, Andy -- Andy Fugard, Postgraduate Research

I wanted to produce a layout like this: |----------| | | | 1 | | | |-----------------------| | | | | 2 | 3 | | | | |-----------------------| Three equal sized plots with the top one centred about the other two. So far I am using nf<-layout(matrix(c(1,1,2,3), 2, 2,

I am a little bit confusing about the following help message on how to fit a GAM model with interaction between factor and smooth terms from http://rss.acs.unt.edu/Rdoc/library/mgcv/html/gam.models.html: ?Sometimes models of the form: E(y)=b0+f(x)z need to be estimated (where f is a smooth function, as usual.) The appropriate formula is: y~z+s(x,by=z) - the by argument ensures that the smooth

Dear all, I'm struggling to find the best (set of?) function(s) to do repeated measures logistic regression on some data from a psychology experiment. An artificial version of the data I've got is as follows. Firstly, each participant filled in a questionnaire, the result of which is a score. > questionnaire ID Score 1 1 6 2 2 5 3 3 6 4 4 2 ...