similar to: help for color parameter

Dear ReXperts, I have the below text file output. I need to extract the T, QC, QO, QO-QC and WT columns for the data between T = 10 and T=150. Any ideas? Thanks in advance. ======================================================================================== 1 D C ---CAT-- T THETA QC QO QO-QC QC/QO WT FSD 8 1 0 1.0000E+01

Dear all: I have been thinking to import below one data file (.txt)into R by read.table(..,skip=1, header=T). But How can I deal with the repeated rows of TABLE NO.1 and names of data variables in the middle of this data file. The similar block will be repeated 100 times, here only show 4 of them and within each block, data records also can vary, here only paste 4 rows for example. I appreciate

Dear all: I have been thinking to import below one data file (.txt)into R by read.table(..,skip=1, header=T). But How can I deal with the repeated rows of TABLE NO.1 and names of data variables in the middle of this data file. The similar block will be repeated 100 times, here only show 4 of them and within each block, data records also can vary, here only paste 4 rows for example. I appreciate

R 2.15.1 OS X and Windows Colleagues, I have a file that looks that this: TABLE NO. 1 PTID TIME AMT FORM PERIOD IPRED CWRES EVID CP PRED RES WRES 2.0010E+03 3.9375E-01 5.0000E+03 2.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 1.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00 2.0010E+03 8.9583E-01

Dear List, I will appreciate any advice regarding how to convert the following numbers [I got in return by taxondive()] in numeric integers without the e.g. 6.4836e+01 abbreviations. Thank you very much in advance, Gian > taxa_dive Species Delta Delta* Lambda+ Delta+ S Delta+ Nat1 5.0000e+00 6.4836e+01 9.5412e+01 6.7753e+02 8.7398e+01 436.99 Nat2

R 2.15.1 OS X Colleagues, I am reading a 1 GB file into R using read.table. The file consists of 100 tables, each of which is headed by two lines of characters. The first of these lines is: TABLE NO. 1 The second is a list of column headers. For example: TABLE NO. 1 COL1 COL2 COL3 COL4 COL5 COL6 COL7 COL8 COL9 COL10

Hello, Are there any R functions available for performing a serial correlation test for short time series (e.g, series having between 10-14 observations)? Many thanks! Isabella R. Ghement, Ph.D. Ghement Statistical Consulting Company 301-7031 Blundell Road, Richmond, B.C., Canada, V6Y 1J5 Tel: 604-767-1250 Fax: 604-270-3922 E-mail: isabella at ghement.ca Web: www.ghement.ca -----Original

Hello, I want to compute the P-value from the joint cumulative distribution of an n-dimensional order statistic in R, using the formula found on http://cmgm.stanford.edu/%7Ekimlab/multiplespecies/Supplement/methods_network.html My data consists of three different techniques (G2D, POCUS and RANDOM), and each has associated with it a number of rankings (integer between 0 and 1000), like for

Greetings, I am working in the sem package on a model with 3 exogenous variables (2 are nominal-categorical), and 4 endogenous, continuous variables. To use sem with the nominal variables, I created dummy variables. Now, in my sem output I have estimates for path coefficients for the relationship between each level of the nominal variables and the endogenous variables they are associated

System info: Linux slackware R Version 1.5.1 (2002-06-17) ESS 5.1.21 Emacs 21.2.1 ------------------- Colleagues I have some text files in a certain directory and want to process them sequentially. So I want to pass the filename identifier into a function and loop through the processing of each file in turn.I've called the file identifier f in this code segment: image.dvm <- function

Dear Rxperts, I have the following data: Study Study.Name C Category TC Time QC QO SD FSD Theta 1 NONE 0 P(22) 0 0.00 7.5596 0 0 8.0361e-03 0 1 NONE 6 G(50) 0 0.00 1.0000 0 0 0.0000e+00 0 1 NONE 2 F(02) 0 0.00 100.0000 0 0 0.0000e+00 0 1 NONE 3 F(03) 0 0.00 13.2280 0 0 1.6732e-02 0

(I tried to send this earlier, but it doesnt seem to have come through, due to problems on my system) Hola: Both cannot be correct: > m <- matrix(1:4, 2) > svd(m) $d [1] 5.4649857 0.3659662 $u [,1] [,2] [1,] -0.5760484 -0.8174156 [2,] -0.8174156 0.5760484 $v [,1] [,2] [1,] -0.4045536 0.9145143 [2,] -0.9145143 -0.4045536 > La.svd(m) $d [1]

Hi, I'm having a problem converting a Matlab program into R. The R code works almost all the time, but about 4% of the time R's optim function gets stuck on a local minimum whereas matlab's fminsearch function does not (or at least fminsearch finds a better minimum than optim). My understanding is that both functions default to Nelder-Mead optimization, but what's different about

Dear all, I am currently learning to run logistic regression models with R. Would someone tell me what this exactly means: Estimated scale (compare to 1 ) 1.746724 If the value is higher or lower than 1, what should I do? The complete results of the model were as follows: Generalized linear mixed model fit using Laplace Formula: success ~ diffca + (1 | partner1) + (1 | partner2) + (1 |

I was using the hist() function to create a frequency table of some network traffic data. The range in values is rather large, from 0 till just under 10e12. Calling hist(x, breaks=c(0,1000,1e6,1e9,1e12),plot=F,freq=T) causes hist() to return : $breaks [1] -1.0000e+05 1.0100e+05 1.1000e+06 1.0001e+09 1.0000e+12 Is this recalculation of the breaks by hist() intended? Maarten van Gelder.

Hi, I am having problems trying to assess the significance of random terms in a generalized linear mixed model using lme4 package. The model describes bird species richness R along roads (offset by log length of road log_length) as a function of fixed effects Shrub (%shrub cover) and Width (width of road), and random effect Site (nested within Site Cluster). >From reading answers to previous

Dear R-Gurus I am a PhD student from South Africa working on chimpanzee behaviour. I am looking at patterns of shade utilization and am using generalized linear mixed models to examine the effects of various factors on whether chimpanzees choose to spend time in the sun or shade. I realise that the lme4 package and the outputs of the lmer functions have been discussed ad nauseum but I have been

I am attempting to model data with the following variables: timepoint - n=48, monthly over 4 years hospital - n=3 opsn1 - no of outcomes total.patients skillmixpc - skill mix percentage nurse.hours.per.day Aims To determine if skillmix affects rate (i.e. no.of.outcomes/total.patients). To determine if nurse.hours.per.day affects rate. To determine if rates vary between

the code is list below , i can run this code for laplace approximation but not importance sampling. the code work properly if use impSamp=0 , but error for any other number >0 . could you guys help me out. thanks so much ! source("http://www.mi.uib.no/~skaug/cash/lesaffre_dat.s") library(glmmADMB) fit_glmmADMB_ip <-

Hi here, i recently notice the software built by otter company know as AD-MB tools . i explore a little bit this tools under the R Package for mixed logistics model. it is a very interesting tool. but i get several questions. first question: i tried the R code below, and have the error , i can run the code for switching the setting to impSamp=0 , it means i can only run the code under