similar to: numeric(0) -> NA

[[this bounced first, because it has 'help' in the Subject line ... -- Martin Maechler ]] I am a complete novice R programmer. (Though I know C quite well) I am trying to write some R code to do the following simulation. There is a 2-frame "movie" of noise and signal dots. the noise dots have random positions in each frame. The signal dots are placed randomly in frame 1,

Here is a programming question. The code I am using is quite slow and I was wondering if there is a way to get rid of the for loop. I am dealing with "interaction" in 2x2 table, and am using Edwards's G_I (Likelihood, p. 194). I label the cells in the table as follows stim response "y" "n" total -------------------------------- y hit miss nsignal

Dear colleague, I not sure this R code is correctly ? I would to show the number of Sample Size at Sample Size Axis that line draw from Power Axis (80%) from R code. How I show this and select the most appropriate of this power (.79955687 - 80983575). Thank for your help and answer. Best Regards, Nikom Thanomsieng, Email: nikom at kku.ac.th .... #Power analysis: Sample size for

Hi there, Could somebody help me disect this legacy R script I inherited at work, I have two questions: 1. I've tried to upgrade our R version from 1.6.2 (yeah, I know), to R 2.0, but some of the lines in this script are not compatible with R 2.0, could someone help me figure out where the problem is? 2. the jpeg generated (attached) seems to be off on some of the data, is there a better way

Ok, I have a new, multipart problem that I need help figuring out. Part 1. I have a three dimensional array (species, sites, repeat counts within sites). Sampling effort per site varies so the array should be ragged. Maximum number of visits at any site = 22 Number of species = 161 Number of sites = 56 I generated the array first by;

  Hello All,   I am working with dlm for the purpose of estimating and forecasting with a Kalman filter model. I have succesfully set up the model and started generating results. Of course, I need to somehow be sure that the results make sense. Without any apparent target to compare with, my natural selection is the results by odinary least square. The idea being that if I choose a diffuse prior,

Hi all, I am using c++ functions in R by Rcpp and inline packages. The code is quite simple, but the R session always automatically crash after some running time. Does anyone here familiar with Rcpp and inline? What¡¯s the problem in the following code? I have checked the input values, no NA and other strange value exists. Thank you for your attention! > mkc <- cxxfunction(

Hi all, I am using c++ functions in R by Rcpp and inline packages. The code is quite simple, but the R session always automatically crash after some running time. Does anyone here familiar with Rcpp and inline? What¡¯s the problem in the following code? I have checked the input values, no NA and other strange value exists. Thank you for your attention! > mkc <- cxxfunction(

Hello everybody, I was wondering if anybody knows how to make channel IDs different if all call are coming from the same host: core show channels Channel Location State Application(Data) SIP/sip-ny1.stanapho 16462257222@in-sip:4 Up Playback() SIP/sip-ny1.stanapho 16462257222@in-sip:4 Up Playback() Thanks in advance. Serge -------------- next part

Hallo, I'm dynamically generating buttons depending on the number of rows of my dataframe. Every button is supposed to call a function which generates a plot with the values of one of my dataframe rows. My code looks like this: base <- tktoplevel() plotten <- function(mat, namen, titel) { midpts <- barplot(height=mat, names.arg = namen, main = titel, las=2)

I am trying to implement a simple r-svm example using the iris data (only two of the classes are taken and data is within the code). I am running into some errors. I am not an expert on svm's. If any one has used it, I would appreciate their help. I am appending the code below. Thanks../Murli ####################################################### ### R-code for R-SVM ### use leave-one-out

Dear R-helpers, I don't see a difference between the following two plots of effect objects, which I understand should be different. What am I missing? require(doBy) require(effects) data(budworm) m1 <- glm(ndead/20 ~ sex + log(dose), data=budworm, weight=ntotal, family=binomial) m1.eff <- all.effects(m1) plot(m1.eff, rescale.axis = FALSE, selection = 2, main = 'rescale =

I have been attempting some basic power calculations using R and I am not getting the results I expect. I had a homework assignment in SAS, but I want to learn R as well, so I was attempting to reproduce my result. (No one else in the class is doing R, so there's no need to obsfucate the answer, the SAS code is what I get my grade for.) The code I am using is: # You assume that the

Dear all, I searched the mail archives and the R site and found no guidance (tried "merge", "cbind" and terms like "coalesce" with no success). There surely is a way to coalesce (like in SQL) columns in a dataframe, right? For example, I would like to go from a dataframe with two columns to one with only one as follows: From Name.x Name.y nx1 ny1 nx2 NA

Thanks; I did not think of that. I have just run a few basic tests with both rsync and OpenSSH in their default settings, when it comes to compression. SSH compression seems to have a very slight edge. However, I surmise that, given the number of knobs available on both sides (OpenSSH, in particular) one can tinker with settings almost endlessly in either side, probably being able to end up with

I sent a message a couple days ago about doing calculations for power of the ANOVA. Several people got back to me very quickly which I really appreciated. I'm working now on a similar problem, but instead of a balanced ANOVA, I have an unbalanced one. The first part of the question was: You assume that the within-population standard deviations all equal 9. You set the Type 1 error rate at á

Hi In the following codes, I had to choose between printing (= TRUE) or deliver something for grab (ei, vi). Is there a way to get both--that is, to print and also have ei and vi for grab? Thanks. Steven ... out<-round(as.data.frame(cbind(ap,se,t,p)),digits) out<-cbind(out,sig) out<-out[!grepl(colnames(zx)[1],rownames(out)),] if(printing){ cat("\nAPPs of bivariate ordered

Dear list I try to run a bootstrap with lmer. I often, but not always, get the error message: Error in objective(.par, ...) : Leading minor of order 6 in downdated X'X is not positive definite (with the number (here 6) varying) In R-archives I came across some threads that treated this problem, nevertheless they refer to lmer when using it with family = "binomial", so the

?s 09:51 de 11/08/2024, Steven Yen escreveu: > Hi > > In the following codes, I had to choose between printing (= TRUE) or > deliver something for grab (ei, vi). Is there a way to get both--that > is, to print and also have ei and vi for grab? Thanks. > > Steven > > ... > > out<-round(as.data.frame(cbind(ap,se,t,p)),digits) > out<-cbind(out,sig) >

Greetings, I sent the message below to the developer of the contributed R package "sspir", but have yet to receive any response. I would be very grateful for any advice people have on the matter. Thanks, Mark -------- Original Message -------- Subject: possible bug in sspir? Date: Tue, 19 May 2009 16:08:41 -0700 From: Mark Scheuerell <mark.scheuerell at noaa.gov> To: