search for: stu

...to be able to pass different fitness functions to this procedure. I am trying to find a good way to do this. I was thinking of passing in the name of the function and then using eval. However, I haven't gotten this to work the way I'd like it to. Consider the following example: > stu <- function(x) {return( 1 + (2*x*x) - (3*x) )}> (x=0:3)[1] 0 1 2 3> stu(x)[1] 1 0 3 10> (fun="stu")[1] "stu"> eval( parse( text = paste( fun, "(", x, ")", sep = "" ) ) )[1] 10> Notice that the function I defined called &quot...

Dear List: I am generating a series of plots iteratively using Sweave. In short, a dataframe is subsetted row by row and variable graphics are created conditional on the data in each row. In this particular case, this code ends up generating 17,000 individual plots. In some cases, all student data (this is working with student achievement data) are available and my code below works very well in the sense that a line connects all points. However, in some cases there are missing data and I need to modify my code so that lines are connected through all points even when data are missing...

...2006), and when I try to calculate the variances for several variables, I get negative numbers. I thought it may be my data, so I ran the example on the help page: data(api) ## one-stage cluster sample dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc) svyvar(~api00+enroll+api.stu+api99, dclus1) variance SE api00 11182.8 1386.4 api00 11516.3 1412.9 api.stu -4547.1 3164.9 api99 12735.2 1450.1 If I look at the full matrix for the variances (and covariances): test<-svyvar(~api00+enroll+api.stu+api99, dclus1) print(test, covariance=T) v...

Hi all! I am working with GLMM using the binomial family I use the following codes I dropped no significant terms, refitting the model and comparing the changes with likelihood: G.1<-lmer(data$Ymat~stu+spi+stu*sp1+(1|ber),data=data,family="binomial") G.1b<-lmer(data$Ymat~stu+spi+(1|ber),data=data,family="binomial") anova (G.1,G.2) But, when I want to evaluate the significance of random effect (1|ber) I cannot use a likelihood-ratio test, probably because the link functio...

...ains an infected object. Sender = samba@lists.samba.org (may be forged) Recipients = marketing@microstar.ru Subject = hi Message-ID = [unknown-id] Antivirus filter report: --- Dr.Web report --- Following virus(es) has been found: infected with Win32.HLLM.MyDoom.32768 Dr.Web detailed report: ns.stu.neva.ru:1332/[text/plain] - Ok ns.stu.neva.ru:1332/message.zip infected with Win32.HLLM.MyDoom.32768 Dr.Web scanning statistic: Infected : 1 --- Dr.Web report --- The original message was stored in archive record named: drweb.quarantine.T0fO52 In order to receive the original message, please s...

...sum( (y - (par[1] * (x1 + 1) * x2^(-par[2]) * x3^par[3])^2 ) } y, x1, x2, x3 are all vectors. I am a bit new to R and I have not been able to find a good description of how to iterate over rows in a data.frame. What is a straightforward way to do this? What am I missing? Thanks Stu [[alternative HTML version deleted]]

..."CG") Error in optim(c(0.66, 0.999, 0.064), pe, NULL, method = "CG") : objective function in optim evaluates to length 6 not 1 I don't know what the error message should mean to me. Particulars for the problem are listed below. Thanks in advance for any help!! Stu > y <- c(0.111111111111111,0.0921052631578947,0.0564516129032258,0.0608108108108108 ,0.0128205128205128,0.0136078431372549) > x3 <- c(3600,169200,185400,255600,271800,342000) > x2 <- c(1,3,4,5,6,7) > x1 <- c(0,0.978723404255319,0.087378640776699,0.549295774647887,0...

...d.J..Yoon = c(0L, 0L, 0L, 0L, 0L, 0L), Eugene.Koontz = c(0L, 0L, 0L, 0L, 0L, 0L), Jakob.Homan = c(0L, 0L, 0L, 0L, 0L, 0L), Kelvin.Kakugawa = c(0L, 0L, 0L, 0L, 0L, 0L), Kirk.True = c(0L, 0L, 0L, 0L, 0L, 0L ), Lars.George = c(0L, 0L, 1L, 0L, 0L, 0L), Soren.Macbeth = c(0L, 0L, 1L, 0L, 0L, 0L), Stu.Hood = c(0L, 0L, 0L, 0L, 0L, 0L), Tibor.Kiss = c(0L, 0L, 0L, 0L, 0L, 0L), Tom.White = c(0L, 0L, 1L, 0L, 0L, 0L), Unassigned = c(0L, 0L, 0L, 0L, 0L, 0L )), .Names = c("Adrian.Cole", "Alison.Wong", "Andrei.Savu", "Bruno.Dumon", "Edward.J..Yoon&quot...

Having issue updating a model with JSON. I simply *"model.attributes=json"*, which does update the attributes but not the id, or uuid in this case. I want to save or create the record (if it''s id exists). Seems like this would be documented, but I''m not finding anything. I''m using a MOM to pass JSON between apps. Thanks. -- You received this message because

...Linux groups for privileges. (I'm assuming Samba can still be used this way.) My question is: How does Samba calculate these group memberships? I'll look through the source code and see what I can see, but I am no programming wizard, so I doubt I'll find what I need there. Thanks! Stu... -- Stuart Reedy Working hard for a great university! stu at coe.uky.edu 859 257-7966 http://www.coe.uky.edu/~stu/

<!doctype html> <html> <head> <meta charset="UTF-8"> </head> <body> <div> <br> </div> <blockquote type="cite"> <div> On 19/05/2020 17:18 Stuart Henderson <<a href="mailto:stu@spacehopper.org">stu@spacehopper.org</a>> wrote: </div> <div> <br> </div> <div> <br> </div> <div> On 2020/05/19 17:04, Aki Tuomi wrote: </div> <...

Dear R-Help, I am using the studres function from the MASS package to compute studentized deleted residuals in a oneway anova. I'm having trouble interpreting the results in situations where a factor level has only one observation. Sometimes studres yields an NaN and sometimes it produces a numeric value for cases where a f...

<!doctype html> <html> <head> <meta charset="UTF-8"> </head> <body> <div> <br> </div> <blockquote type="cite"> <div> On 19/05/2020 16:48 Stuart Henderson <<a href="mailto:stu@spacehopper.org">stu@spacehopper.org</a>> wrote: </div> <div> <br> </div> <div> <br> </div> <div> On 2020-05-19, David Gessel <<a href="mailto:gesse...

On 2018-08-27, Stuart Henderson <stu at spacehopper.org> wrote: > On 2018-08-27, Zach Cheung <kuroro.zhang at gmail.com> wrote: >> After upgrading my VMware Fusion (10.1.3) Arch Guest to the latest with >> OpenSSH upgraded from 7.7p1 to 7.8p1, found that ssh from macOS Sierra >> (10.1...

...put all the related scripts I use in the same folder with a "master" script and then source() them in that "master" script. Problem is, the "master" script must first know where it is (without me having to open it and retype the path every time I move it). Any ideas? Stu ps. I found this on the internet but a) I couldn't understand it, and b) it didn't work: frame_files <- lapply(sys.frames(), function(x) x$ofile) frame_files <- Filter(Negate(is.null), frame_files) PATH <- dirname(frame_files[[length(frame_files)]]) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~...

...care about slowing down other traffic flows. H.323 packets should be serviced first no matter what. (2) Any traffic with source or destination public IP Address "abc.def.ghi.jkl" or private IP Address 192.168.16.4 is to/from from the videoconference device. After all the reading and studying, is it really this simple? Does anyone have ideas on how to test this? How do I watch packets to see which packets go into what PFIFO band? #!/bin/sh VTC1_PRIVATE="192.168.16.4" VTC1_PUBLIC="abc.def.ghi.jkl" IPTABLES="/usr/local/sbin/iptables" $IPTABLES -t m...

Hi all, My data genes [,1] [,2] [,3] [,4] [1,] 25 72 23 55 [2,] 34 53 41 33 [3,] 26 43 26 44 [4,] 36 64 64 22 [5,] 47 72 67 34 stu<-t.test(genes[,1:2],genes[,3:4]) > stu$p.value [1] 0.4198002 i get 1 pvalue for the entire col1:col2 Vs col3:col4. I am trying to get 5 p values for the 5 rows i have. I am trying to avoid a for loop coz my actual data has 56000 rows and its taking more than 4 minutes to compute. Thanks i...

...# vector of interest A = xyz(vec,-1) # creates a 6x6 matrix with vec on the subdiagonal where xyz is some function similar to diag, but with differing arguments. I can't believe there is not a simple way to do this... Thanks for your help, ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Stu Field, PhD Postdoctoral Fellow Department of Biology Colorado State University 1878 Campus Delivery Fort Collins, CO 80523-1878 Office: E208 Anatomy/Zoology Phone: (970) 491-5744 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ [[alternative HTML version deleted]]

...constraints: 0.0 < c2 < 1.0, and 0.0 < c3 < 1.0 y, x1, x2, and x3 are observed data. I have a total of 6 rows of data in a data set. Is "optim" in the stats package the right function to use? Also, I can't quite figure out how to specify the constraints. Thank you! -Stu

...i and ci to (hopefully) force par2 and par3 into the range 0.0 <= par2 <= 1.0 and 0.0 <= par3 <= 1.0. I am not at all sure that ui and ci are defined correctly. The call to constrOptim returns normally, but without a solution (apparently). Any suggestions appreciated. Thanks Stu # data y <- c(0.111111111111111,0.0921052631578947,0.0564516129032258,0.0608108108108108 ,0.0128205128205128,0.0136078431372549); x1 <- c(0,0.978723404255319,0.087378640776699,0.549295774647887,0.0596026490066225 ,0.61578947368421); x2 <- c(1,3,4,5,6,7); x3 <- c(3600,169200,1854...