similar to: Peculiar behavior of attached objects

I''m a physicist working on fusion energy and dabble in statistics only occasionally, so please excuse gaps in my statistical knowledge. I''d appreciate any help that a real statistics expert could provide. Most people in my field do only very simple statistics, and I am trying to extend some work on multivariate linear regression to account for significant between-group

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)

Any thoughts on the following? Am I missing something? [I am running R for Windows 0.63.1 on a Win NT 4.0 (workstation, sevice pack #4) box.] Thanks. FC. R : Copyright 1998, The R Development Core Team Version 0.63.1 Beta (Dec 5, 1998) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type "?license" or

Hello Group, I am having trouble with the sqldf package on unix. The same code works fine on windows. Silly Example script: # Load the package library(sqldf) # Use the titanic data set data(women) colnames(women) head(women) sqldf('select height, count(*) from women where height is not null group by weight') Unix Output and error: bash-3.00$ R --vanilla <testR.R

I am wanting to use the "[" operator in an S-chunk, e.g. <<>>= str(women) women$height women[,1] "["(women,1) @ to show the equivalence of three methods of extracting an element from a data.frame. However Sweave returns the last of these as women[1] in the S input chunk How can I force it not to do this and return "["(women,1)

Goodmorning, This is a documentation related question about the B-spline function in R. In the help file it is stated that: "df degrees of freedom; one can specify df rather than knots; bs() then chooses df-degree-1 knots at suitable quantiles of x (which will ignore missing values)." So if one were to specify a spline with 6 degrees of freedom (and no intercept) then a basis

I generally run 'make; make check' (with more settings) when building the Debian package. Running 3.2.4 rc from last night, I see a lot of package loading issues during 'make check'. Here is splines as one examples: checking package 'splines' * using log directory '/build/r-base-3.2.3.20160303/tests/splines.Rcheck' * using R version 3.2.4 RC (2016-03-02 r70270) *

Hi With ggplot2 I can very easily create beautiful histograms but I would like to put two histograms on the same plot. The histograms may be over-lapping. When they are overlapping the bars are shown on top of each other (so that the overall height is the sum of the two). Is there any way to get them to display overlapping (with smaller value in front, larger value behind) so that the overall

Hi, I would be grateful for help with a problem which is irritating me. I am quite sure that I am doing something stupid, but I can't see what it is. I am running R 1.7 on Windows 2000. The graphics device is the PC screen. The graphics from the nlme demonstration in Bates an Pinheiro's manual work just as advertised. The CO2 data and the Orthodont data dsiplay

----- Mike Hammett Intelligent Computing Solutions http://www.ics-il.com -------------- next part -------------- An HTML attachment was scrubbed... URL: http://lists.digium.com/pipermail/asterisk-users/attachments/20070905/c62f4465/attachment.htm

As part of carrying out a complicated maximum likelihood estimation, I am trying to learn to program likelihoods in R. I started with a simple probit model but am unable to get the code to work. Any help or suggestions are most welcome. I give my code below: ************************************ mlogl <- function(mu, y, X) { n <- nrow(X) zeta <- X%*%mu llik <- 0 for (i in 1:n) { if

My provider migrated from an old EOL softswitch to Trixbox. I have a number (8159093011) on a different server on a different network. It appears as though the incoming calls are trying to authenticate against that number, which isn't present on the box. Could someone help me decode this debugging output? I was calling 8159911010. My server is 208.100.1.33. Theirs is 208.1.87.235. I

I need to setup a small mail server on a local network. It only needs SMTP ability as it's just so Asterisk can send out emails. The machine has sendmail installed. My primary mail server seems to be rejecting the messages. Some research says something isn't configured properly. What do I have to do so the outside world accepts emails from my Asterisk box? It is behind a NAT.

One of my providers has a different SIP account for each number. I have all of my users in one outbound context (caller ID passes fine). How do I ensure that the callers get routed down their correct SIP account with my provider for e911 purposes without each having their own context? ----- Mike Hammett Intelligent Computing Solutions http://www.ics-il.com -------------- next part

My provider has a Coppercom switch. I have included the authentication information they gave me. How would I structure this in Asterisk to the registration and the entry in sip.conf? User Name - 8159093010 Password - XXXXX No Pin Proxy - sip.essex1.com (10.1.3.2) Outbound Proxy - proxy.essex1.com (63.164.210.14) Change setting to use "outbound Proxy" ---------- Mike Hammett

I'm writing a script to sync some mp3 files. Due to a limitation in the number of destination files that can be read from my thumb drive, I'm not looking to preserve the original file structure (actually, I'm looking to sync *only the files* to the new destination directory). The source files are all subfolders under /backup/Music: ./Adrian Legg/Mrs. Crowe's Blue Waltz/Paddy

I have two Asterisk servers. One is my primary server that I link to all of my providers and the other is at an office building with multiple tenants. If I tell Asterisk to dial an entry in the iax.conf that is for one customer off that second box, why does it use a different account for a different customer? It still ends up at the correct box, but it is hard to troubleshoot issues when

I thought I read somewhere that you no longer have to do a special make command for the 2.6 kernel. Is this true, or should I still make linux26? I'm having problems getting anything zaptel to load properly. ---- Mike Hammett Intelligent Computing Solutions http://www.ics-il.com -------------- next part -------------- An HTML attachment was scrubbed... URL:

I have multiple upstreams in my office. The primary upstream is having some issues with latency\jitter. I want to move the VoIP traffic to another interface. I have the router set to send all traffic destined for "local" networks out the respective interfaces. Traffic destined to the Internet goes out one of the upstreams. I can do this on a per-IP basis and have successfully done

Hello there, I have two groups (men and women) and I know per group how many of them smoke or don't smoke (women 40 of 200; men 100 of 300). I would like to know how I can compare in R if men and women differ significantly in their smoking. However, because there are more men in the sample than women I cannot just compare the number of smokers and non-smokers in both groups, right?! (I would