similar to: looping through a data frame

Hi, I''ve set up puppet to get node definitions from LDAP as per the docs. It''s been working well, but I now want to use fully qualified domain names instead of simple domain-less hostnames for the node name. I replaced the ou=Hosts entries with equivalent ones using FQDNs, restarted the puppetmasterd, and tried a "puppetd --test" from one of the nodes. However, I

Hi, I am doing an analysis to see if these is tissue specific effects on the gene expression data . Our data were collected from 6 different labs (batch effects). lab 1 has tissue type 1 and tissue type 2, lab 2 has tissue 3, 4,5,6. The other labs has one tissue type each. The 'sample' data is as below:

Hi, First post to the list, so please feel free to set me straight if I'm not following some protocol or other :o) We need to use rsync to send files to a client, being a Windows user, we decided to try both the cwrsync implementation and also a straight cygwin/rsync install. I'm experiencing the following errors: using cygwin implementation: $ sh clientupload.sh

I realize batch mode is still experimental, but I was hoping there might be a workaround for a problem I am getting. I have been trying to run some tests and I get the below error when I use the --read-batch option to. I can successfully create an initial set of batch files, then a second set based upon a few modified test files from the first batch. When I first run the --read-batch option

Hi, I am new to R and am still trying to get a grip of it. I have data in this format: Run Lab Batch Y 1 1 1 1 608.781 2 2 1 2 569.670 3 3 1 1 689.556 4 4 1 2 747.541 5 5 1 1 618.134 6 6 1 2 612.182 7 7 1 1 680.203 8 8 1 2 607.766 9 9 1 1 726.232 and I want to make a boxplot of the Y values for each

Dear R People: I'm having some trouble with mpiexec and Rmpi. I would like to be able to pass in the number of "children" via the mpiexec command (from the command line). this is in SUSE10.1, with R-2.7.1 Here are my files: cat eb.R library(Rmpi) mpi.remote.exec(paste("i am",mpi.comm.rank(),"of",mpi.comm.size())) mpi.quit() hodgesse at

I am having trouble printing a table out to the GUI display when the table is created and printed within a loop. I get a "Error: syntax error message" If I comment out the print statement, the loop runs fine and I can print out the last iteration of the table. ...[multiple loops and calculations ending with.....]... +

Hi, I really feel I've looked everywhere, although I know this can't be a hard problem. I'd like to be able to call the graph below as a function, but I can't get the function to recognize variables beyond 'dframe'. I've read through many papers on writing functions in R, but I can't get this to work. data <- data.frame('date' = as.Date(rep(c(15101,

I've created a short program to print a table of learning curve factors. However, I cannot figure out how to format the table to: 1) Get rid of the [1]s in the first column and replace it with the values of N. 2) Line up the first row with the factors (decimal fractions). Thanks for any help. The complete program and output is as follows: > Lc<-seq(0.70,0.95,0.05) #Specify learning

Hello, I am trying to filter a data set like below so that the peaks in the Phase value are more obvious and can be identified by a peak finding function following the useful advise of Carl Witthoft. I have written the following for(i in length(data$Phase)){ newphase=if(abs(data$Phase[i+1]-data$Phase[i])>6){ data$Phase[i+1] }else{data$Phase[i] } } I get the following error which I have not

Hello all, I am working with a 1000x1000 matrix, and I would like to return a 1000x1000 matrix that tells me which value in the matrix is greater than a theshold value (1 or 0 indicator). i have tried mat2<-as.matrix(as.numeric(mat1>0.25)) but that returns a 1:100000 matrix. I have also tried for loops, but they are grossly inefficient. THanks for all your help in advance. Lanre

Dear all, Is there any function for calculating confidence limits for coefficients in an lm() object? I know of the confint() function in the MASS library working very well on my binomial GLMs and I have tried it (using glm () , family=gaussian) but it gives NAs according to below. Does the confint() function not accept gaussian GLMs? Could there be convergence problems in the GLM? Note the

Hi, When I run the survfit function, I want to get the restricted mean value and the standard error also. I found out using the "print" function to do so, as shown below, print(km.fit,print.rmean=TRUE) Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier") records n.max n.start events *rmean *se(rmean) median 200.000

Hi All, Please pardon me if I am missing something obvious here. How do I get the Kaplan-Meier estimate function that is created by survfit and plotted by the code. fit <- survfit(Surv(time, status) , data=aml) plot(fit) That is, I need a function that will give me the survival estimate at a given time: \hat{S}(t). Thanks in advance. Ritwik Sinha ritwik.sinha at gmail.com | +12033042111 |

Hello, Thank you to all those great folks that have helped me in the past (especially Dennis Murphy). I have a new challenge. I often generate time-series data sets that look like the one below, with a variable ("Phase") which has a series of flat-topped peaks (sample data below with 5 "peaks"). I would like to calculate the phase value for each peak. It would be great to

Hi I compouted a multiple linear regression with repeated measures on one explanatory variable: BOLD peak (blood oxygenation) as dependent variable, and as independent variables I have: -age.group (binaray:young(0)/old(1)) -and task-difficulty measured by means of the reaction-time 'rt'. For 'rt' I have repeated measurements, since each subject did 12 different tasks. -> so

Hello I've been looking at the tables generated by `SequenceToOffsetTable::emit`, and notice that when the generated data are strings, the data is basically un-grep-able, and very tricky to read, as they are emitted as an array of comma-separated char-literal: extern const char HexagonInstrNameData[] = { /* 0 */ 'G', '_', 'F', 'L', 'O',

Hi, I've got two suggestions how to speed up median() about 50%. For all iterative methods calling median() in the loops this has a major impact. The second suggestion will apply to other methods too. This is what the functions look like today: > median function (x, na.rm = FALSE) { if (is.factor(x) || mode(x) != "numeric") stop("need numeric data")

Dear Thomas, The head of my dataset > head(wsuv) parcel sp time censo treatment species 1 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1 2 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1 3 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1 4 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1

Hello, I am trying to run a linear regression analysis on my data set. For some reason most variables are removed due to singularities. My linear regression looks this way (I am using only partial data, which is selected by flags): fm<-lm(log(cplex6.time..sec..[flags]) ~ cplex6.cities[flags] + log(1/features.meanOver.frust[flags]) + log(1/features.meanOver.minDist[flags]) + [...]