similar to: R puts '+' within my numbers

Does anyone know of an R-function or method to compare two C-statistics (Harrells's C - rcorr.cens) obtained from 2 different models in partially paired datasets (i.e. some similar and some different cases), with one continuous independent variable in each separate model? (in a survival analysis context)? I have noticed that the rcorrp.cens function can be used for paired data. Thanks

Dear R users, Does anyone know if there's a median filter available in R? I have a considerable amount of images as 256 by 256 matrices and want to smooth them with a 5x5 median filter. Until now I'm not lucky in searching the R-help list/files. Thanks, Hanneke ir. J.M. (Hanneke) Schuurmans (PhD student) Department of Physical Geography Faculty of Geosciences, Utrecht University, The

Dear R-users, I use unlist of POSIX dates to extract the year, hour etc. With that I can search for files in my database which are in the form 'yyyymmddhh_synops.txt' However, I get stucked during midnight where unlist just gives NA's. The script is given below, the problem accurs at acc.period[16] (midnight). However when I write out the character, unlist works well. But

Dear R-users, For a data frame (say in this example X) I want to look up the corresponding value in a 'look-up data frame' (in this example Y). The for-loop works but is very time-consuming because 'X' in reality is very big. Therefore I would like to have a solution with apply. However, I do not succeed. Any suggestions? Thanks in advance, Hanneke

Hello: I would like to be able to do a 3D scatter plot from 3 variables, 2 independent and 1 dependent. The closest R function I could find for this is "cloud". However cloud uses, as input, a matrix where the value of each matrix element is the dependent variable value at that matrix coordinate. My problem is that the independent variable values are floating point and can be of

Hello all! I have a problem to extract values greater that for example 1820. I try this code: x[x[,1]>1820,]->x1 Please help me! Thank you! The data structure is: structure(c(2.576, 1.728, 3.434, 2.187, 1.928, 1.886, 1.2425, 1.23, 1.075, 1.1785, 1.186, 1.165, 1.732, 1.517, 1.4095, 1.074, 1.618, 1.677, 1.845, 1.594, 1.6655, 1.1605, 1.425, 1.099, 1.007, 1.1795, 1.3855, 1.4065, 1.138, 1.514,

Dear R-users, I'm experiencing a segmentation fault when using hdf5load(file,load=FALSE). Library(hdf5) loads without problems but when loading a file, R crashes. I compiled R under Unix (Solaris for Sun). There is nothing wrong with the files, as I can run the same script at another place where R runs under Linux. Is it possible it has something to do with the hdf5 libraries where

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

Hello all! I have a problem to arrange data in another form. My initial data is like this: 'data.frame': 421 obs. of 58 variables: $ 01A: num NA NA NA NA NA NA NA NA NA NA ... $ 01B: num NA NA NA NA NA NA NA NA NA NA ... $ 03A: num NA NA NA NA NA NA NA NA NA NA ... $ 03B: num NA NA NA NA NA NA NA NA NA NA ... $ 05A: num NA NA NA NA NA NA NA NA NA NA ... $ 05B: num NA NA NA NA

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

Dear R users, I am using lmer combined with AIC model selection and averaging (in the MuMIn package) to try and assess how isotope values (which indicate diet) vary within a population of animals. I have multiple measures from individuals (variable 'Tattoo') and multiple individuals within social groups within 4 locations (A, B, C ,D) crucially I am interested if there are

Hi, I have a problem with the df. I read in a big csv file. Tabelle <- read.csv("C:\\Users\\Public\\Documents\\Bachelorarbeit\\eingabe8_durchnummeriert.csv" , header = T , sep=";") then I try this: > ygamma <- glm(Tabelle$sb_ek_ber ~1+ Tabelle$FAHRL_C + Tabelle$NUTZKREIS + Tabelle$schw_drittel_c   , family = Gamma) >  anova(ygamma, test="Chisq")

This is on Solaris 10 x86, do not see this behavior on Solaris 10 sparc. Seen on multiple machines. Sshd debug: debug1: server_input_channel_open: ctype session rchan 256 win 16384 max 16384 debug1: input_session_request debug1: channel 0: new [server-session] debug2: session_new: allocate (allocated 0 max 10) debug3: session_unused: session id 0 unused debug1: session_new: session 0 debug1:

Dear all, I?m having some problems getting optim with method="Nelder-Mead" to work properly. It seems like there is no way of controlling the step size, and the step size seems to depend on the *difference* between the initial values, which makes no sense. Example: f=function(xy, mu1, mu2) { print(xy) dnorm(xy[1]-mu1)*dnorm(xy[2]-mu2) } f1=function(xy) -f(xy, 0,

Hi all Assume I have a data set xx; Group: 1=group1 ?, 2=group2 IQ: ?1= High, 0 =low fit <- glm(IQ ~group, data = xx, family = binomial()) summary(fit) Results ?????? ????????????Estimate Std. Error z value Pr(>|z|) (Intercept) -2.55456??? 0.210 -12.273? < 5e-16 *** group????????? 0.36180 ?????0.076?? 3.952 ????5.24e-05 *** the odd ratio = exp(0.36180 )= 1.435912 My question

Dear forum members, How can I force a negative slope in a linear regression even though the slope might be positive? I will need it for the purpose of determining the trend due reasons other than biological because the biological (genetic) trend is not positive for these data. Thanks. Julia Example of the data: [1] 1.254 1.235 1.261 0.952 1.202 1.152 0.801 0.424 0.330 0.251 0.229

I have a dataframe called "data" with 5 records (in rows) each of which has been scored on each of many variables (in columns). Five of the variables are named var1, var2, var3, var4, var5 using headers. The other variables are named using other conventions. I can create a new variable called var6 with the value 15 for each record with this code: > var6=var1+var2+var3+var4+var5

Been trying to use the short format for a collection in rhtml, i.e. <% collection.each {|x| puts x} %> Doesn''t seem to work though [outputs nothing] ? -- Posted via http://www.ruby-forum.com/.

In simple RHTML template a want to show values based on examples from rails docs. In first example ''puts "Frodo"'' sends the value to console but not to browser. In second example I want avoid a huge of <%%> so I decided to put entire cycle into one <%%> and print values to browser with ''puts''. Still the same results. Can someone clarify,

Bugs item #17697, was opened at 01-02-2008 09:07 You can respond by visiting: http://rubyforge.org/tracker/?func=detail&atid=218&aid=17697&group_id=35 Category: Incorrect behavior Group: current Status: Open Resolution: None Priority: 3 Submitted By: Nobody (None) Assigned to: Nobody (None) Summary: "puts Frame.show" crashes application Initial Comment: I know that nobody