similar to: Sampling with unequal probabilities

Hi ,everyone. I want to draw a boxplot figure, main script shown as below. when I run the command "boxplot(bjavee~yearname)",it is wrong with the error information " model.frame.default(formula = bjavee ~ yearname) : variation 'bjavee' list is wrong". I want to add the yearname as the x-axis tick. And yearname is a array with the length 12,while bjavee is a 5*12

Hi ,everyone. I want to draw a boxplot figure, main script shown as below. when I run the command "boxplot(bjavee~yearname)",it is wrong with the error information " model.frame.default(formula = bjavee ~ yearname) : variation 'bjavee' list is wrong". I want to add the yearname as the x-axis tick. And yearname is a array with the length 12,while bjavee is a 5*12

Dear all, I have a problem with the command 'getURL' from the RCurl package, which I have been using to obtain a ftp directory list from the MOD16 (ET, DSI) products, and then to download them. (part of the script by Tomislav Hengl, spatial-analyst). Instead of the list of files (from ftp), I am getting the complete html code. Anyone knows why this might happen? This are the steps i

CentOS Errata and Security Advisory CESA-2009:0256 firefox security update for CentOS 5 x86_64: https://rhn.redhat.com/errata/RHSA-2009-0256.html The following updated file has been uploaded and is currently syncing to the mirrors: x86_64: updates/x86_64/RPMS/firefox-3.0.6-1.el5.centos.i386.rpm updates/x86_64/RPMS/firefox-3.0.6-1.el5.centos.x86_64.rpm

Can someone help me with the if loop below? In the subroutine, I initialize all of (joint12,marg1,marg2,cond12,cond21) as FALSE, and call with only one of them being TRUE: ,...,joint12=FALSE,marg1=FALSE,marg2=FALSE,cond12=FALSE,cond21=FALSE,,,, joint12 seems to always kick in, even though I call with, e.g., marg1 being TRUE and everything else being FALSE. My attempts with if... else if were

Is something wrong in the initialisation part that we don't see? joint12 <- marg1 <-F marg1 <-T if (joint12) { print ("joint 12") cat (joint12) } if (marg1) { print("marg 1") cat(marg1) } Would probably be my diagnostic approach On Fri, 9 Aug 2024, 04:45 Steven Yen, <styen at ntu.edu.tw> wrote: > Can someone help me with the if loop below?

The following (using if else) did not help. Seemed like joint12 always kicked in. ??? me1<-me0<-NULL. ??? if(joint12){ ????? {me1<-cbind(me1,v1$p12);? me0<-cbind(me0,v0$p12)} ??? } else if(marg1) { ????? {me1<-cbind(me1,v1$p1);?? me0<-cbind(me0,v0$p1)} ??? } else if(marg2) { ????? {me1<-cbind(me1,v1$p2);?? me0<-cbind(me0,v0$p2)} ??? } else if(cond12){ ?????

Thanks. Hmm. The loop is doing what it is supposed to do. > try1<-function(joint12=FALSE,marg1=FALSE,marg2=FALSE, +??????????????? cond12=FALSE,cond21=FALSE){ + # *************************************************** + # Testing if loop + # *************************************************** + if(joint12){ +?? {print ("joint12"); cat(joint12,"\n")} +?? {print

OK. The fact it's in a function is making things clearer. Are you trying to update the values of an object from within the function, and have them available outside the function. I don't speak functional programming articulately enough but basically v <- 1 funA <- function() { v <- v+1 } funA() cat (v) # 1 You either return the v from the function so funB <- function() {

Greetings, I''d like to remove the cruft (no longer maintained) areas of yum.puppetlabs.com. I fear they only cause confusion and offer less-than-desirable experience for our users. I''d like to remove /base /prosvc /porsvc.unsigned /SRMS /sources I''d like to do this some time next week. The stuff mostly found in base is available elsewhere and signed properly. The

[PLEASE REPLY _OFF_ THE LIST, i.e. DON'T CC to r-help at ...] Hi, I don't see this sort of thing very often on the mailing lists, so list moderators and others should feel free to tell me if it breaches list etiquette and/or delete my post if necessary. But I can't see what harm it could do... I am just wondering approximately how many people use / have used some of the R stuff

"Or use <<- assignment I think. (I usually return, but return can only return one object and I think you want two or more" You can return any number of objects by putting them in a list and returning the list. Use of "<<-" is rarely a good idea in R. -- Bert On Fri, Aug 9, 2024 at 1:53?AM CALUM POLWART <polc1410 at gmail.com> wrote: > > OK. The fact

What is the correct URL for the SRMS for packages found in: http://mirror.centos.org/centos-6/6/SCL/ ? -- *** E-Mail is NOT a SECURE channel *** James B. Byrne mailto:ByrneJB at Harte-Lyne.ca Harte & Lyne Limited http://www.harte-lyne.ca 9 Brockley Drive vox: +1 905 561 1241 Hamilton, Ontario fax: +1 905 561 0757 Canada

Hi: I am practicing with the attached shapefile and was wondering if I can get some help. Haven't used 'rgdal' and 'maptools' much but it appears to be a great way bring map data into R. Please take a look at the comments and let me know if I need to explain better what I am trying to accomplish. library(rgdal) library(maptools) library(ggplot2) dsn="C:/Documents and

Stop advertising. Le 26 sept. 2010 09:46, "Faisal Hanif" <faisal at vopium.com> a ?crit : Hi Abdul-Basit, If you need only different intervals of billing you can easily do it using any AGI as we are doing it in Perl AGIs using post call billing. But if you need realtime billing then the most stable and flexible option is to use FastAGI+ AMI. I have tested it in JAVA and it

POTS line already did. It's like a zombie channel. When this happens we issue "hangup request DAHDI/x-x". Anyone faced similar issue? Current configuration: Ubuntu 10.04 64 bits, Asterisk 1.6.2.7, DAHDI 2.3.0. Thank you, Alex Saavedra --0016364d1dc91459020496ed0bbe Content-Type: text/html; charset=UTF-8 Content-Transfer-Encoding: quoted-printable <div

I am using levelplot() from lattice with grids that have unequal cell sizes. This means that the boundary between two cells is not always half-way between nodes, as levelplot() assumes. The result is that some cell sizes are rendered incorrectly, which can be painfully obvious if using relatively large cells. Is there any work-around? I am using the conditioning capability of lattice and

I ask for a histogram of counts with slightly uneven breaks but histogram() ignores the 'type' argument and provides density scaling instead. x = sample(1:3, 100, replace=TRUE) histogram( ~ x, breaks=c(0,1.5,2.5,3.5), type='count') My real application has time data with month boundaries for the breaks and I DO want counts on the y-axis. A work-around is not difficult but

Dear R folks - Sorry to be so dense, but can you help me with two programming problems? 1. Randomizing within factors. I have factored the values in an array according to a series of ranked breakpoints. Now I want to randomize the values that fall within each of the breakpoint boundaries, without replacement. I want to end up with an array in which the values within each factored category

Hi I am trying to run an ANOVA on data with unequal variance. I am new to nlme, but to my understanding I need to use the gls function. I have single response variable (distance which is continuous) and the explanatory variable is individual ID (class variable: individuals differ in the variance in their distance values hence the need to using nlme). So I would create a model