search for: ph1

Hi, I have a data.frame with many variables for which I am performing the mean by subgroup, for a pair of variables at a time, where one of them for each pair defines the subgroup. The subgroups in the x$cm1 are 0, 1 and 2. x ph1 cm1 0.2345 2 1.2222 1 2.0033 0 0.0000 2 1.0033 1 0.2345 0 1.2222 2 2.0033 0 0.0000 1 1.0033 2 > meanbygroup <- as.vector(with(x, by(x$ph1, x$cm1, mean))) > meanbygroup if the ph1 has no missing values the above statements work fine: [1] 1.4137000 0.7418333 0.6150000 In the moment that...

Hi folks. If I have the following in my "data" event pH1 pH2 1 4.0 6.0 2 4.3 5.9 3 4.1 6.1 4 4.0 5.9 and on and on..... for about 400 events Is there a way I can get R to plot event vs. pH1 and event vs. pH2 and then do a loess or lowess line for each?? Thanks in advance David [[alternat...

...} if ( ( r = ao_close(device) ) != 1 ) { fprintf(stderr,"Error in ao_close()\n") ; exit(EXIT_FAILURE) ; } ao_shutdown() ; return 0 ; } char *Sine(int sr, int nchan, double dur, int *nBytes) { double f1, f2 ; double db, A ; double ph1, ph2, phInc1, phInc2 ; int frames, i ; short *sbuffer ; char *cbuffer ; f1 = 1000 ; // frequency in Hertz f2 = 1.5 * f1 ; // perfect fifth higher ph1 = ph2 = 0 ; phInc1 = ( 2 * M_PI * f1 ) / sr ; // phase increments phInc2 = ( 2 * M_PI *...

...stupid to create the correct subset for phase 2. Phase1: Sample 1000 parts with 80 defective. Phase2: Sample 100 parts out of these 1000 with 15 defective. 0:ok, 1:defunct. The table below gives the conditional sampling values. Please help me: library(survey) ss1 <- data.frame(id=1:1000, ph1.x=rep(c(1,0),c(10,990)), subset=rep(c(1,0),c(100,900)), ph2.y=rep(c(1,0,NA),c(15,85,900)), n1=rep(1000,1000), n2=rep(100,1000) ) table(ss1$ph1.y, ss1$ph2.x) > Phase1.x >Phase2.y 0 1 > 0 85 0 > 1 5 10 p2 <- twophase(id=list(~id,~id), strata=list(NULL,NULL)...

Hello, In the following I tried 3 versions of an example in R help. Only the two first predict command work. After : library(splines) require(stats) 1) fm1 <- lm(weight ~ bs(height, df = 5), data = women) ht1 <- seq(57, 73, len = 200) ph1 <- predict(fm1, data.frame(height=ht1)) # OK plot(women, xlab = "Height (in)", ylab = "Weight (lb)") lines(ht1, ph1) 2) height <- women$height # 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 weight <- women$weight # 115 117 120 123 126 129 132 135 139 14...

...ecurring graphics issue that I've not been able to resolve with R. If I make a series of regression estimates and then plot the estimated function for the regression lines over a scatter plot of the data, e.g., using a sequence of plot( ) and lines ( ) similar to those below plot(dep12spp13ph1$DAYSWETm2,dep12spp13ph1$AFTERDECOMP,pch=dep12spp13ph1$spp,cex=0.75,ylab="Proportion of biomass after leaching", xlab="Number of days in wetland",xlim=c(0,250),ylim=c(0,1)) xplot<- 0:243 xplotbsk64deg1<-bs(xplot,knot=62,degree=1,Boundary.knots=c(0,243)) lines(xplot,exp(xp...

...be gratefull. My macro: [macro-phonesystems] exten => s,1,NoOp(We are calling=${ARG1}) exten => s,2,GotoIf($[${GROUP_COUNT(ph0)}>=1]?100:3) exten => s,3,Set(GROUP()=ph0) exten => s,4,Dial(Sip/${ARG1:1}@ovh1,40,TwW) exten => s,5,NoOP(PH0) exten => s,100,GotoIf($[${GROUP_COUNT(ph1)}>=1]?200:101) exten => s,101,Set(GROUP()=ph1) exten => s,102,Dial(Sip/${ARG1:1}@ovh2,40,Tw) exten => s,103,NoOp(PH1) exten => s,200,GotoIf($[${GROUP_COUNT(ph2)}>=2]?300:201) exten => s,201,Set(GROUP()=ph2) exten => s,202,Dial(Sip/${ARG1:1}@phonesystems1,40,Tw) exten =>...

Dear friends: It's possible make a legend in the coplot output? My data are: "dados"<-structure (list ( Qu1 = c(1,2,3,4,5,6,7,8,3,2), Qu2 = c(3,4,5,6,7,8,9,9,6,3), Ku1 = structure (factor ( c(1,1,1,1,2,1,2,2,1,1),levels=1:2),.Label=c("Kmin","Kmax")), Ph1 = structure (factor ( c(1,1,3,2,3,3,2,3,2,1),levels=1:3),.Label=c("Pmin","Pmean","Pmax")), Th1 = structure( factor( c(1,1,2,4,3,3,2,4,2,1),levels=1:4),.Label=c("Tmin_in","Tmax_in","Tmin_out","Tmax_out"))), row.names = c("At1...

Hi I've 3 iax connections to my provider , each of them have own DID , PH1<----| | \/ PH2<-->|-----| <---------------------------> |----|<-- DID1 | A1 | <---------------------------> |ISP |<-- DID2 PH3<-->|-----| <---------------------------> |----|<-- DID3 I had iax phone on each of this connection , but now I wa...

...to this effect. # check model assumptions mcp.fnc(m3) # check significance of model terms pamer.fnc(m3) # Only the intercept is significant. For purposes of the # example, let’s perform a posthoc analysis on FreqBdc on # model m2. m2.ph <- mcposthoc.fnc(model = m2, var = list(ph1 = "FreqBdc")) if(try(require(languageR,quietly=TRUE))){ # Generate plot for Length X Freq with function plotLMER.fnc require(languageR) plotLMER.fnc(m2, pred = "LengthBc", intr = list("FreqBdc", levels(eeg$FreqBdc), "beg", list(1 : 3, 1 : 3)))...

...ave read all the documentation and samba passes all the diagnoses tests. In fact samba always functions correctly for the subnet it's on. Any suggestions? Thanks, Beau Atwater -- Beau Atwater beau@bwa.bellcore.com Sensors, Systems, and Telemetry ph1 732-758-3373 Bell Communications Research ph2 908-630-9620 fax 732-758-2804

...----------------------------------- # No tengo para nada claro lo que hace. No hace falta que nadie se lo mire exprofeso. Pero si ya conoce como funciona y le interesa me puede comentar # http://www.mitpressjournals.org/doi/pdfplus/10.1162/jocn_a_00143 m1.posthoc <- mcposthoc.fnc(m1, var = list(ph1 = "Affix", ph2 = "group"), mcmc = F, nsim = 100, mc.cores = 1) names(m1.posthoc) as.matrix(names(m1.posthoc$var)) as.matrix(names(m1.posthoc$summaries)) levels(eeg$LengthBd) summary(m1.posthoc) -- [[alt...