search for: ex1

Hi there fellow R-users, Im seeing some strange behaviour when I multiply a vector by a matrix Here is my script: > tr 1 2 3 4 5 6 0.2217903 0.1560525 0.1487908 0.1671354 0.1590643 0.1471667 > > ex1 a b c d e f 1 0.2309579 -3.279045 -0.6694697 -1.1024404 0.2303928 -1.5527404 2 -0.2865357 -2.519789 -0.1138712 -0.3571832 -2.6727913 -0.3296945 > is.vector(tr) [1] TRUE > is.data.frame(ex1) [1] TRUE > tr* ex1[1,] a...

...;- c(rep(1970,3),rep(1980,3)) > x <- c(20,30,40,25,35,45) > ex <- data.frame(id=id,year=year,x=x) > ex id year x 1 A 1970 20 2 B 1970 30 3 C 1970 40 4 A 1980 25 5 B 1980 35 6 C 1980 45 I want to draw a subset of "ex" by selecting only the A and B units: > ex1 <- subset(ex[which(ex$id=="A"|ex$id=="B"),]) Now I want to do some computations on x for each unit: > tapply(ex1$x, ex1$id, mean) A B C 22.5 32.5 NA But this gives me an NA value for the unit C, which I thought I had already left out. How do I ensure that t...

...ar <- c(rep(1970,3),rep(1980,3)) > x <- c(20,30,40,25,35,45) > ex <- data.frame(id=id,year=year,x=x) > ex id year x 1 A 1970 20 2 B 1970 30 3 C 1970 40 4 A 1980 25 5 B 1980 35 6 C 1980 45 I want to draw a subset of "ex" by selecting only the A and B units: > ex1 <- subset(ex[which(ex$id=="A"|ex$id=="B"),]) Now I want to do some computations on x for each selected unit only: > tapply(ex1$x, ex1$id, mean) A B C 22.5 32.5 NA But this gives me an NA value for the unit C, which I thought I had already left out. How do I en...

...r the fact that expression objects behave somewhat ``wierdly'', IMHO. We can have expressions of length(.) > 1, and access elements of these like list elements. However, other natural things don't work with expressions; some example code: ##-- `Simple' length--1 expressions: ex1 <- expression(hat(beta) == (X^t * X)^{-1} * X^t * y) ex2 <- expression(bar(x) == sum(over(x[i], n), i==1, n)) ex3 <- expression(paste(over(1, sigma*sqrt(2*pi)), " ", plain(e)^{over(-(x-mu)^2, 2*sigma^2)})) ##-- A length--2 expression: ex12 <- expression(bar(x) == sum(over...

...bit easier. I thought it might be that the shell doens't exist, but I tried making a user with a invalid shell and get back error code 1, so its not that. Ideas? Cheers Miles -----Original Message----- From: Miles Roper Sent: Monday, 3 February 2003 08:54 a.m. To: 'MCCALL,DON (HP-USA,ex1)' Cc: 'samba-technical@lists.samba.org'; 'samba@lists.samba.org'; Esh, Andrew; Ronan Waide; STEFFENS,MICHAEL (HP-Germany,ex1); 'Richard Sharpe'; 'John H Terpstra' Subject: RE: [Samba] RE: Winbind on HPUX11, Totally Stuck, Please Help Thanks for your help, still...

I've noticed with certain datasets that when I try to do an anova and test for main effects and interaction for two explanatory variables, sometimes the main effect results are given but not the interaction results. For example, ex1 = aov(Score ~ var1*var2, data=myData) summary(ex1) gives me only the main effects for var1 and var2, but not the interaction. I also tried ex1 = aov(Score ~ var1+var2+var1:var2, data=myData) summary(ex1) but it still only gives the main effects. The only way I can get the interaction results is...

...r seperator for winbind) B) Is it case sensitive? in other words do I specify JWAD+dantest rather then jwad+dantest? c) Can I specify an NT user group such as all accounts in the JWAD Domain Users account and how would I do this? Thanks, -Dan -----Original Message----- From: MCCALL,DON (HP-USA,ex1) [mailto:don_mccall@hp.com] Sent: Thursday, February 14, 2002 10:49 AM To: 'Thomas, Daniel J.'; MCCALL,DON (HP-USA,ex1) Subject: RE: [Samba] Winbind - Why won't you authenticate??? Hi Daniel, Does Solaris have a getpwnam or getpwent command, that you can see if your libnss_winbindd.s...

Warning: This is only for those interested in R language minutiae A recent post on this list asked if there was a simple way to change the R language object: ex1 <- expression(x < a) ## just the part to the right of the <- assignment to the object expression( x < a & y < b) ## or something like this Phil Spector showed how to do this by essentially deparsing and reparsing the string: parse(text =paste(ex1, "& y < b"...

Hi, I am trying to reconcile anova table in R (summary(lm)) with individual t.test. datafilename="http://personality-project.org/R/datasets/R.appendix1.data" data.ex1=read.table(datafilename,header=T) #read the data into a table summary(lm(Alertness~Dosage,data=data.ex1)) gives: Call: lm(formula = Alertness ~ Dosage, data = data.ex1) Residuals:    Min     1Q Median     3Q    Max -8.500 -2.437  0.250  2.687  8.500 Coefficients:             Estimate Std. E...

I've got two tables.... first one(table1): ID chrom start end Ex1 2 152 180 Ex2 10 2000 2220 Ex3 15 3000 4000 second one ( table2): chrom location name 2 160 Alv 2 190 GNN 2 100...

...and store it in a numeric vector (of length 100). The for loop should be placed in a function that allows the user to vary the sample size, the simulation size, the integration limits and the input function. In addition construct a histogram of the estimates and include a red vertical line at pi. ex1.fcn<-function(x){ h<-4/(1+x^2) return(h) } n=1000 a=0 b=1 my.rand.x=runif(n,min=a,max=b) pi.MC = ((b-a)/n)*sum(ex1.fcn(my.rand.x))

...0,0.5]. The function follows a parabola and has two local maxima located at the margins of the interval [0,0.5]. Please try the following code which performs this optimisation, produces a plot of the function within the given interval and adds the computed maximum returned by 'optimise()' ex1 <- function(x) log(4 * x^2 - 2 * x + 0.5) opt1 <- optimise(ex1, lower=0, upper=0.5, maximum=T) x <- seq(0,0.5,0.001) plot(x, ex1(x), type="l") abline(v=opt1$maximum, col="red") The returned value of the maximum is 0.309017 which is definitely wrong. Example 2: The ne...

Well, I managed to get Samba 2.2.3 up and running on our Solaris 8 machine. I installed with the winbind option and everything went though just find. I was able to join the NT domain and now I can do a wbinfo -u "and get a domain user list as well as a "wbinfo -g and get a group list. For some reason though, the authentication isn't working. I tried to "wbinfo -a" and

...rscores in the middle. We hope this works - there is a strong assumption here on just how complicated your indices are, for example. We are assuming they are either identifiers (not using two successive underscores in their names) or numbers. If they are not, you need to get even craftier. > Ex1 <- D(Ex, "A__2") > Ex1 A__1 * (exp(A__2 * X) * X) > Ex1s <- deparse(Ex1) > Ex1s [1] "A__1 * (exp(A__2 * X) * X)" > pat1 <- paste(Ident, "__", Subsc, sep = "") > rep1 <- "\\1\\[\\2\\]" > Ex1ss <- gsub(pat1, rep1, Ex1s...

Hello I'm currently using the dbFD function of the FD package and i'm having some things that I can't do. Is there any way to check the relations between dbFD indexes? Function cor for example? I can't manage to put the informations correctly dbFD function gives a lot of output (indexes - nbsp, sing.sp, FRic, FEve, FDiv, FDis and RaoQ). I want to see the relationships between

...h/Sweave/Sweave-manual.pdf and i have tried to replicate the example on pages 4-5, but i encounter the following problem: IT DOES NOT WORK. I can run Sweave, create the .tex file and then when building the .tex file in TeXnikCenter i get all sorts of error messages. The syntax: >Sweave("ex1.Rnw") Writing to file ex1.tex Processing code chunks ... 1 : echo term verbatim 2 : term verbatim eps pdf You can now run LaTeX on 'ex1.tex' I have attached the log file with the error messages. Initially i had a problem with \usepackage {sweave} but i think i have solved (as it...

...figured ical for zimbra, caldav for google mail and ews for exchange 2010 calendar. ical and caldav setup working fine and i am getting my calendar events perfectly. But for exchange 2010 calendar i am getting following error. "Unable to communicate with Exchange Web Service at ' https://ex1.domain.com/EWS/Exchange.asmx': Could not authenticate to server: ignored NTLM challenge, GSSAPI authentication error: Unspecified GSS failure. Minor code may provide more information: Credentials cache file '/tmp/krb5cc_0' not found" my calendar.conf is as follows [calendar3] ty...

I am a new Samba user so forgive my ignorance. I am using Samba version 2.0.6 on HP UX 11 and I need to use the "security mask" parameter to prevent even the owner of a file from changing the read only attribute (by right clicking on the file from a Windows NT 4.0 desktop) if the file is read only (Unix permissions: 0440). This is what I have in the "smb.conf" file: [global]

...th ms applications hanging on trying to open a file. Turns out it was an nfs locking issue. Hope this helps, Don -----Original Message----- From: Partenheimer Nate [mailto:npartenh@butler.edu] Sent: Monday, March 26, 2001 1:23 PM To: "MCCALLDON "@envelope.hp.com; HP-USA@envelope.hp.com; ex1 " Subject: RE: RE: Explorer & MS Office hang Sensitivity: Confidential Thanks for all your help - I finally solved it! We were, coincidentally, having some trouble with sas on our enterprise server too so after speaking to sas support we determined that our nfs lock data had become corr...

Hi, I apparently was asleep at the wheel using scp, and accidentally copied a file onto itself. Scp generated an Input/Output error and did not perform the copy. However, now the file is corrupt. Is this a but? Is there anyway to fix the file I messed up? Thanks, Don Gathman 208.396.6675