Questions on the use of Mathematica to construct models for approximating empirical data.

learn more… | top users | synonyms (3)

3
votes
1answer
124 views

Ordinary Least Squares (OLS) with RLink

I need to calculate the parameters b0, b1 and b2 in: ...
0
votes
1answer
102 views

Binomial Regression where n > 1

I am not well versed in Generalized Linear Models, but am trying to use this approach to model some discrete data. I have binomial data where the sample size for each observation is the same size, ...
0
votes
1answer
91 views

Regression with 3 variables

I'm trying to find a function starting from its values but I can't find anything related to functions with 2 variables and I don't know how to go on. So, given this sample data, I'd like to know how ...
0
votes
1answer
70 views

How to fit this data to be used in NIntegrate later?

I have a list of data points that resemble a chi-squared distribution (long tailed gaussian) that I am having trouble fitting or finding an equation with. The goal is to use the fitted equation in an ...
1
vote
1answer
227 views

FindFit::nrlnum: error with FindFit

I am trying to write a code to fit actual data of the orbit of a star with Mathematica 8. First I tried to write a code to fit a mock orbit, using the Application with differential equations (the ...
3
votes
1answer
497 views

Test Goodness of fit for FindFit

Is there any way to get an estimate of the goodness of fit from FindFit, for instance getting the final error? The documentation doesn't provide any hints, and a quick Google search didn't return any ...
2
votes
2answers
253 views

FindFit working with Integrate but not NIntegrate - I am using “?NumericQ”

I'm having some syntax problems with finding a fit to an equation that doesn't have an analytical expression for it's integral. Because I think the problem is one of syntax, I am substituting a ...
2
votes
0answers
83 views

Is there a way to produce a FittedModel object as the solution to a regularized regression?

I can get the Ridge or Tikhonov regression solution to Y = \[Beta].X easily enough as ...
3
votes
0answers
159 views

The CorrelationMatrix is equal to corrcoef of Matlab?

How can I generate a matrix which is equivalent to the results of the "correlation coefficients" command (corrcoef) in Matlab? When I make a LinearModelFit, I can get a CorrelationMatrix, but is this ...
11
votes
0answers
345 views

Implementation of Partial Least Squares (PLS)?

I'm looking for a PLS implementation for Wolfram Mathematica. I just can't seem to find any implementation for it - do you know one? In general, I'm a bit confused that there exist only so few ...
1
vote
0answers
135 views

Fitting curves with NonLinearModelFit not working with charge and discharge like data

I'm trying to get a reasonable fit out of mathematica, but i'm not very practical with it. I'm not that good with mathematica at all actually(like cleaniness of the code...). btw, here's my problem: ...
2
votes
1answer
124 views

FindDistributionParameters and Hypothesis Testing

Please excuse me if I am making any logical error in posting this question since I am not a math major. However, let me post my question. I am trying to fit multivariate parametric distribution to ...
2
votes
1answer
158 views

How do I find out what a Mathematica function is doing (which numerical method it uses and how)?

Long story short, I am working on a project which entails fixing a bug I have found in Mathematica (8 and 10). In general, I would like to know what a specific Mathematica function is actually doing ...
1
vote
1answer
100 views

“Power series” in Cos[n*x] for X coordinate of planetary ellipse

Consider the functions below: ...
1
vote
2answers
231 views

How to obtain the best possible NonlinearFit

Here is some sample data ...
4
votes
1answer
229 views

Find GeometricTransform from 3D object external surfaces to a positional point cloud data

I am looking for another way to approach this problem. Any suggestion on how to optimize the computing speed is welcome. I have a BSplineFunction which accurately ...
1
vote
1answer
184 views

Mixture coefficients and Parameter estimation for moving Normal distributions

I am searching for a way to do the following steps below. I am stumped and I have stumped a few profs too. Any help would be appreciated. I would like to import data from, say, excel or a .txt file ...
0
votes
2answers
69 views

Handling msg for fitting

I have a trivial question: I do Monte Carlo with some curve fitting. Now always when a random run creates a complex fitting solution it should write 1 in a Matrix. How can I handle such msgs in a ...
4
votes
1answer
120 views

Mystifying failure of NonlinearModelFit

Fitting the following two time-series (two repeats of the same exponential growth process) ...
3
votes
1answer
164 views

Troubles with FindFit for ODE

I need to use FindFit for ODE with arbitary number of parameters. This code (from documentation) works: ...
0
votes
1answer
123 views

Fitting data with inclusion of constraints

How do I find a fit for the dataS (see below) with the constraint that the parameter np is an integer, and that the parameters la+lc = 100. A weak attempt at this is provided below for a specific ...
16
votes
2answers
667 views

Curve fitting by running through hundreds of models and return the one with best fit

In the past, I used to solve a lot of my regression analysis (curve fitting) problems with a program called DataFit which runs on Windows. It has hundreds of regression models which it runs through in ...
3
votes
1answer
543 views

Fitting data to a circle with FindFit

I'm currently trying to use the Mathematica's FindFit command to fit a curve of some data to the bottom half of a circle. The data is of the form ...
5
votes
1answer
255 views

Forcing NonLinearModelFit through a point

Like the title says, is there a way to make an NLM fit pass through a specific point.
1
vote
1answer
103 views

Using List/Histogram Data As A Weighted Function?

I have a list of data points that I would like to use as a weighted function. The data represents the number of particles within a bin of a particle distribution. Ultimately I would like to multiply ...
2
votes
1answer
70 views

Is there a way to show some “weighted” residuals after LinearModelFit?

I am fitting some functions to some data and then I want to analyze the residuals of that fit. I'd like to see how those residuals compare to the values of the original data, hence my "weighted" in ...
0
votes
2answers
135 views
2
votes
3answers
173 views

BSpline or other curve through arbitrary ordered sequence of two-dimensional points

How does one create a BSplineCurve or other smooth curve passing through each point in an ordered list, where the curve may be non-simple, i.e., can double-back ...
4
votes
2answers
557 views

Best fit equation from a scanned image

I am trying to extract a curve from a scanned graph and find a best fit equation. For example, starting from: How could one find an equation starting from the image file ?
1
vote
0answers
161 views

Error associated with parameters in Minimize (or NMinimize) in discrete fitting

my problem goes as follows: I need to a discrete fit of a discrete function. Say I have the following made up data points, which are generated by the parabola ...
0
votes
0answers
323 views

Multi-peak fitting

Suppose we have a data set that looks like this: I have to fit this data to 9 peaks knowing μ and σ beforehand & have to ...
1
vote
0answers
86 views

Updating procedural fitting algorithm to more efficient style?

For some time now, I have used this procedural programming approach to fit data in my research that is too cumbersome for the built in Mathematica functions. (Although this simple example works fine ...
0
votes
1answer
253 views

Equation of line from loglog plot

Given some data that looks linear on a ListLogLogPlot, how can I find the correct equation (and particularly, the correct slope of the line)? I have a sample code here ...
0
votes
0answers
127 views

How to perform a distribution fit test to censored data

With following code I will get an error: ...
4
votes
2answers
164 views

Remove redundant parameters from equation

I have random expressions like (b1 x + b2 x)/(b3 + b4/b5 x) + Sin[b6 x] where $x$ is a variable and {b1, ..., b6} are parameters to be fitted. This expression is ...
0
votes
1answer
173 views

Non linear model fit

I need to fit this, but for some reason I cannot... data = {{0.617,0.8}, {0.605,0.6}, {0.5997,0.4}, {0.5972,0.2}, {0.5985,0.1}}; soln = NonlinearModelFit[data, {y - 2478.82*wlog[kc + 1]}, {y,w,k}, ...
10
votes
0answers
274 views

PenaltyFunction Option in NMinimize

I am curious if anybody has some experience in using the built-in Option "PenaltyFunction" in any of the numerical optimization functions like NMinimize. According to the documentation the user ...
7
votes
1answer
201 views

How to derive SinglePredictionBands from a NonlinearModelFit model

I have found NonlinearModelFit to be extremely useful in my work. Here I have fit 4 points to a function and am being asked to explain how the 95% SinglePredictionBands associated with the plot below ...
1
vote
0answers
90 views

Simultaneous fitting to multiple data sets with multiple variables and some common parameters [duplicate]

I am trying to fit two data sets to a function model described below which share a few common parameters. Say I have two functions with some common parameter(s) ...
2
votes
1answer
114 views

Struggling with FindFit

I'm trying to fit a sinusoidal wave to some numerical data. I can do it by hand using manipulate, so I know roughly the sort of parameter values that I should get back, but functions such as FindFit ...
6
votes
2answers
425 views

2D smoothing spline interpolation

Does Mathematica have 2D smoothing spline interpolation built in? I requires an interpolation method with smooth first derivatives and cubic bivariate splines fulfill this nicely. In python I would ...
6
votes
3answers
3k views

How can I obtain a smooth curve from my data but preserve two inversion points

I have a list of data points, {x, y}, my data. When I plot them, the curve is jagged. I want to smooth the curve and preserve the two sharp corners. This is a plot ...
13
votes
2answers
870 views

What is the difference between FindFit and NonlinearModelFit

I cannot get the difference between these two fitting functions in Mathematica. I even read in MMA that the result of these two are the same.
13
votes
1answer
890 views

Curve fit based on minimal data

Inspired by the fantastic answers here and here, I would like to ask a questions along similar lines. I have 6 curve types that follow these rules: Type A: If point 1 is positive and point 2 is ...
14
votes
3answers
2k views

Fitting ellipse to 5 given points on the plane

Five points are required to define a unique ellipse. An ellipse has five degrees of freedom: the $x$ and $y$ coordinates of each focus, and the sum of the distance from each focus to a point on the ...
3
votes
1answer
149 views

Saving data from NonlinearModelFit

If you use NonlinearModelFit the results are put into some sort of output package. How do you save these packages so that you can go back to them later? My particular problem is that I generating ...
1
vote
0answers
172 views

Comparing two ways of fitting systems of ODEs to data

I'd like to ask you if you see any difference between the two following ways of estimating parameters. Given my system of ordinary differential equations that I want to use to find the values of all ...
1
vote
1answer
72 views

Why does my Dynamic NonliearModelFit keep running?

I have complex data and which I am fitting to a complex formula. In order to use NonlinearModelFit I split the formula and the data into real and imaginary parts. ...
5
votes
3answers
906 views

Finding exponential model from data [duplicate]

I'm trying to find an exponential model from data for a homework problem. There is an accompanying video explaining how to do the bigger problem but it omits any instructions on the step where the ...
1
vote
1answer
330 views

General F statistic and p value in LinearModelFit

Many statistical packages provide a "F statistic" and the corresponding p value in their output of a multiple linear regression. I don't find that in the parameters of the Mathematica function ...