similar to: Re : Frequency count

Hiya, I have been struggling to save the result from the FOR loop. What is the best way to do it, as I need the result to merge with another dataset for further analysis ? for (dd in ((M-10):M)){ + dist<-(32-dd) + r<-1/2*(1-exp(-2*dist/100)) + map<-c(dd,round(r,4)) + print(map) + next + } Thanks. Stella ___________________________________________________________________________ This

Dear members, I want to create a sequence of numbers for the multiple records of individual animal in my dataset. The SAS code below will do the trick, but I want to learn to do it in R. Can anyone help ? data ht&ssn; set ht&ssn; by anml_key; if first.anml_key then do; seq_ht_rslt=0; end; seq_ht_rslt+1; Thanks in advance. Stella

Dear R list, I am reading in text file data prepared in Access database by someone. One of the field contains timestamp data, how can I separate the timestamp data into two varaibles: date and time. Can I specify the field is in timestamp format when I first reading in ? My reading in data are as below: 449 LWT 22/10/2003 15:43:00 441 143 449 LWT 17/11/2003 15:25:00 421 169 449 LWT

Hi, I am trying to achieve 3-way conferencing taking hint from wiki link http://www.voip-info.org/wiki/view/Asterisk+n-way+call+HOWTO Here is the scenario: 1. user "ua1" calls user "ca1" 2. "ua1" then presses the feature code "*0" to redirect "ca1" to conference room 300 3. "ua1" then dials the user "33" 4. user

Dear lists, I want to construct a loop in R, but don't know how to do it. I can do it in SAS, but I prefer in R (which I am hoping I will off SAS for good soon). Could anyone help me to convert the SAS codes to equivalent R codes. Basically, the following codes were written to establish the sire gametes or phases for daughter design for one markers two alleles. Here are the SAS code: do

Hi, i am using the n-way-call dialplan solution found on voip-info. i have added its entry in applicationmap of features.conf file. the problem is......its not working. to activate the n-way call i dial *0 but nothing happens. i have played around with dtmf and codec settings but no success. the extensions and sip configuration is below if you want to have a look. I dont have any clue why its not

Hello! I need to use n-way call as it described here: http://habrahabr.ru/sandbox/52259/ It is in russian, but dial plan is quite clear. It works in asterisk 11: -- Blind transferring OOH323/7272-6385 to '0' (context fromtransfer) priority 1 -- Executing [0 at fromtransfer:1] NoOp("OOH323/7272-6385", "") in new stack -- Executing [0 at fromtransfer:1]

Dear R community, In the excellent nlme package the default diagnostic plot graphs the innermost residuals against innermost fitted values. I recently fit a mixed-effects model in which there was a very clear positive linear trend in this plot. I inferred that this trend occurred because my fixed effect was a two-level factor, and my random effect was a 12-level factor. The negative residuals

I don't get the same result, do you have a package loaded that would change the default behavior (such as Hmisc)? > list(grp.1) [[1]] [1] 1 2 Levels: 1 2 > list(grp.1[mask]) [[1]] [1] 1 Levels: 1 2 > library(Hmisc) <<snip>> > list(grp.1[mask]) [[1]] [1] 1 Levels: 1 --Matt > version _ platform i386-pc-mingw32 arch i386 os

Hi all, I'm using R 1.5.1 on Windows 2000. The following snippet of code doesn't seem to do anything - no error is reported, and there is no name change. names(myFrame[,c(1:3)]) <- c("name1", "name2", "name3") This code however works nicely: names(myFrame)[c(1:3)] <- c("name1", "name2", "name3") Can anyone suggest why

Hello, I have been searching for ways to obtain these for combinations of fixed factors and levels other than the 'baseline' group (contrasts coded all 0's) from a mixed-effects model in lme. I've modelled the continuous variable y as a function of a continuous covariate x, and fixed factors A, B, and C. The fixed factors have two levels each and I'd like to know whether

Greetings all, I'm trying to figure out how to calculate the inverse CDF (i.e. a quantile) for a non-central F distribution. I could put together a quick numerical solver routine using the CDF, but I wonder if there's a function that I've missed that would be more efficient? Thank-you, Andrew Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885

I would appreciate some thoughts on using the bootstrap functions in the library "bootstrap" to estimate confidence intervals of ICC values calculated in lme. In lme, the ICC is calculated as tau/(tau+sigma-squared). So, for instance the ICC in the following example is 0.116: > tmod<-lme(CINISMO~1,random=~1|IDGRUP,data=TDAT) > VarCorr(tmod) IDGRUP = pdLogChol(1)

I have two lists: xa <- list( X=c(1,2,3), Y=c(4,5,6), Z=c(7,8,9) ) xb <- with( barley, tapply( X=seq(1:nrow(barley)), INDEX=site , FUN=function(z)yield[z])) I can convert xa to a dataframe easily with: as.data.frame(xa) But if i try the same with xb I get: as.data.frame(xb) Error in as.data.frame.default(xb) : can't coerce array into a data.frame What

Why does stop("we are done") print "Error in eval.with.vis(expr, envir, enclos) :" ? It would seem to me that a plain stop() is not an error, and that it would make more sense to have an error() function that is different from a stop(). Is there a rationale here that I am missing? sincerely, /iaw

Hi All, I download the app_asyncgoto.c, compile the app_asyncgoto.so. Then according to this page http://www.voip-info.org/wiki/view/Asterisk+n-way+call+HOWTO ; when I dial ,there have this warning: -- Executing AsyncGoto("SIP/111-086497c8", "SIP/113-08674628|dynamic-nway|111|1") in new stack Feb 2 16:53:10 DEBUG[4218]: app_asyncgoto.c:95 asyncgoto_exec: Attempting

Hello, is there any function to print (using "cat" or "print" etc.) a table of values row per row in a exact order, even if the value has 3 digit and in the next row eg. 5 digits? This means I want to print a number with 3 digits with 2 spaces in front and a number with 4 digits only with one space in front. regards Andreas

Dear R Users What's the easiest way to load an excel data file to R? Thanks Rick --- Outgoing mail is certified Virus Free. Checked by AVG anti-virus system (http://www.grisoft.com). Version: 6.0.351 / Virus Database: 197 - Release Date: 19/4/2002 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

I have a 12" G4 iBook with the latest os (Panther) with the recent updates. When I start R by double clicking, it expands for a fraction of a second and then, nothing happens. There is mention in the readme about using something called "i-install" to uninstall a mistakenly installed library, but I do not know how to do this. Any help would be appreciated. Thanks, Martin

Dear R people, I have a naive question: after fitting "lm" to a data, I can't extract the pvalue corresponding to a specific covariate in a direct way. Could anyone give me a hint? Thank you very much. Frank