cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 
Try the Materials Informatics Toolkit, which is designed to easily handle SMILES data. This and other helpful add-ins are available in the JMP® Marketplace
Choose Language Hide Translation Bar
joaquimgui
Level II

Plot 3D Diagram

Hello,

I am trying to plot a 3D diagram but unsuccessful. I have the data in Excel like this:

 Y4,1581575710,5410784
X Z1Z2
0,00038225 0,047405620,09242527
0,00032384 0,13448339

0,34250103

But I have a lot of rows and columns. This is just an example. I would appreciate any help.

Thanks.

 

1 ACCEPTED SOLUTION

Accepted Solutions
txnelson
Super User

Re: Plot 3D Diagram

I created a sample Excel file, based upon your description and have attached it to this entry.  I then developed the following JSL to read it in, and to created the 3D plot of the data.  I hope this is what you want

3dplot.PNG

Names Default To Here( 1 );

// Read in Excel file
dt = Open(
	"<your path to the excel file>\example.xlsx",
	Worksheets( "example" ),
	Use for all sheets( 1 ),
	Concatenate Worksheets( 0 ),
	Create Concatenation Column( 0 ),
	Worksheet Settings(
		1,
		Has Column Headers( 0 ),
		Number of Rows in Headers( 1 ),
		Headers Start on Row( 1 ),
		Data Starts on Row( 1 ),
		Data Starts on Column( 1 ),
		Data Ends on Row( 0 ),
		Data Ends on Column( 0 ),
		Replicated Spanned Rows( 1 ),
		Replicated Spanned Headers( 0 ),
		Suppress Hidden Rows( 1 ),
		Suppress Hidden Columns( 1 ),
		Suppress Empty Columns( 1 ),
		Treat as Hierarchy( 0 ),
		Multiple Series Stack( 0 ),
		Import Cell Colors( 0 ),
		Limit Column Detect( 0 ),
		Column Separator String( "-" )
	)
);

// Get Y Matrix
Ymat = {};
For( i = 3, i <= N Col( dt ), i++,
	If( Column( i )[1] != "",
		Ymat = Ymat |/ Matrix( Num( Column( i )[1] ) )
	)
);

// Get rid of unrequired cells, and get the X and Z matrices
dt << select rows( [1, 2, 3] );
dt << delete rows;
dt << delete column( 2 );

// Convert all columns to numerics
For( i = 1, i <= N Col( dt ), i++,
	Column( dt, i ) << data type( numeric )
);

// Get X matrix
Xmat = Column( dt, 1 ) << get as matrix;
dt << delete column( 1 );

// Get Z matrix
Zmat = dt << get as matrix;

// Create new data table, and output the data

dtFinal = New Table( "Converted", New Column( "X" ), New Column( "Y" ), New Column( "Z" ) );

For( Xcnt = 1, Xcnt <= N Rows( Xmat ), Xcnt++,
	For( Ycnt = 1, Ycnt <= N Rows( Ymat ), Ycnt++,
		For( Zrowcnt = 1, Zrowcnt <= N Rows( Zmat ), Zrowcnt++,
			For( Zcolcnt = 1, Zcolcnt <= N Cols( Zmat ), Zcolcnt++,
				dtFinal << Add Rows( 1 );
				dtFinal:X[N Rows( dtFinal )] = Xmat[Xcnt];
				dtFinal:Y[N Rows( dtFinal )] = Ymat[Ycnt];
				dtFinal:Z[N Rows( dtFinal )] = Zmat[Zrowcnt, Zcolcnt];
			)
		)
	)
);
dtFinal << Scatterplot 3D(
	Y( :X, :Y, :Z ),
	Frame3D(
		Set Grab Handles( 0 ),
		Set Walls( 1 ),
		Set Rotation( -54, 0, 38 ),
		Set Text Scale( 2 ),
		Set Line Scale( 1 ),
		Set Marker Transparency( 1 ),
		Set Marker Quality( 0.25 ),
		Set Marker Scale( 1 ),
		Set View Perspective( 0.1 ),
		Set View Zoom( 1 )
	)
);

 Be sure to change the path in the JSL to point to where ever you save the attached Excel file

Jim

View solution in original post

11 REPLIES 11
ian_jmp
Level X

Re: Plot 3D Diagram

Generally, JMP Platforms expect the values for each variable (in this case 'x', 'y' and 'z') to be in their own column. So for your data I think this would be:

Screen Shot 2017-11-23 at 15.27.55.png

You could then try, for instance 'Graphs > Scatterplot 3D'.

You can also use JMP to rearrange the data into this shape, either by hand or using a script. But the details will depend on exactly how you read it in from Excel.

txnelson
Super User

Re: Plot 3D Diagram

More explanation of your data needs to be made.

Are you expecting all of the data to be displayed in one plot, or does each Z column, Z2,Z2, etc. denote a separate plot to be made?

Since each displayed point in a 3D plot is defined by the 3 dimensions, can you please show how your values in the displayed table would be displayed if the table was in the form of

X Y Z
     
     
     
     
Jim
joaquimgui
Level II

Re: Plot 3D Diagram

Thank you for replying. I will show you how my data looks like in Excel:

 Y4,1581575710,541078414,645384719,76265229,588463340,103585152,772737964,566401471,633986680,179196392,5639742108,325272114,711013121,815695142,186092160,139901175,938007189,208612200,409857209,675249221,03791241,182479256,870356
X Z1Z2Z3Z4Z5Z6Z7Z8Z9Z10Z11Z12Z13Z14Z15Z16Z17Z18Z19Z20Z21Z22Z23
0                        
-3,8225E-05 0,047405620,092425270,232599291,382111412,041801783,068162310,276730710,088082520,857564830,115119720,397465550,222907370,068690820,484059791,492759330,649649970,852558622,70090820,924209380,0712172120,39353912,75357372,59458619
0,00032384 0,134483390,342501030,379900330,384103631,706275090,893493470,421779430,206233381,519176790,262837270,680637520,752848430,100525930,30084610,981382041,023368370,318973931,473025680,785446060,16804080,33508951,471713417,46166212

 

Below you can see how my data looks like in form X,Y,Z but I have 23 points from Y and hundreds rows.

yy1y2y3y4
x1z11z21z31z41
x2z12z22z32z42

 

I would like to display the coordinates in one 3D plot.

That's the problem, because to show you in a displayed table X Y Z, I have to expand this matrix. For example:

x1y1z11
x1y1z12
x1y1z13
   
x1y2z21
x1y2z22
x1y2z23

And so on.. But I have no clue how to do this in JMP and to plot it. I know there is an automatic way to do this using Origin, but I need to do using JMP. Could you understand what I mean?

Thanks.

txnelson
Super User

Re: Plot 3D Diagram

I created a sample Excel file, based upon your description and have attached it to this entry.  I then developed the following JSL to read it in, and to created the 3D plot of the data.  I hope this is what you want

3dplot.PNG

Names Default To Here( 1 );

// Read in Excel file
dt = Open(
	"<your path to the excel file>\example.xlsx",
	Worksheets( "example" ),
	Use for all sheets( 1 ),
	Concatenate Worksheets( 0 ),
	Create Concatenation Column( 0 ),
	Worksheet Settings(
		1,
		Has Column Headers( 0 ),
		Number of Rows in Headers( 1 ),
		Headers Start on Row( 1 ),
		Data Starts on Row( 1 ),
		Data Starts on Column( 1 ),
		Data Ends on Row( 0 ),
		Data Ends on Column( 0 ),
		Replicated Spanned Rows( 1 ),
		Replicated Spanned Headers( 0 ),
		Suppress Hidden Rows( 1 ),
		Suppress Hidden Columns( 1 ),
		Suppress Empty Columns( 1 ),
		Treat as Hierarchy( 0 ),
		Multiple Series Stack( 0 ),
		Import Cell Colors( 0 ),
		Limit Column Detect( 0 ),
		Column Separator String( "-" )
	)
);

// Get Y Matrix
Ymat = {};
For( i = 3, i <= N Col( dt ), i++,
	If( Column( i )[1] != "",
		Ymat = Ymat |/ Matrix( Num( Column( i )[1] ) )
	)
);

// Get rid of unrequired cells, and get the X and Z matrices
dt << select rows( [1, 2, 3] );
dt << delete rows;
dt << delete column( 2 );

// Convert all columns to numerics
For( i = 1, i <= N Col( dt ), i++,
	Column( dt, i ) << data type( numeric )
);

// Get X matrix
Xmat = Column( dt, 1 ) << get as matrix;
dt << delete column( 1 );

// Get Z matrix
Zmat = dt << get as matrix;

// Create new data table, and output the data

dtFinal = New Table( "Converted", New Column( "X" ), New Column( "Y" ), New Column( "Z" ) );

For( Xcnt = 1, Xcnt <= N Rows( Xmat ), Xcnt++,
	For( Ycnt = 1, Ycnt <= N Rows( Ymat ), Ycnt++,
		For( Zrowcnt = 1, Zrowcnt <= N Rows( Zmat ), Zrowcnt++,
			For( Zcolcnt = 1, Zcolcnt <= N Cols( Zmat ), Zcolcnt++,
				dtFinal << Add Rows( 1 );
				dtFinal:X[N Rows( dtFinal )] = Xmat[Xcnt];
				dtFinal:Y[N Rows( dtFinal )] = Ymat[Ycnt];
				dtFinal:Z[N Rows( dtFinal )] = Zmat[Zrowcnt, Zcolcnt];
			)
		)
	)
);
dtFinal << Scatterplot 3D(
	Y( :X, :Y, :Z ),
	Frame3D(
		Set Grab Handles( 0 ),
		Set Walls( 1 ),
		Set Rotation( -54, 0, 38 ),
		Set Text Scale( 2 ),
		Set Line Scale( 1 ),
		Set Marker Transparency( 1 ),
		Set Marker Quality( 0.25 ),
		Set Marker Scale( 1 ),
		Set View Perspective( 0.1 ),
		Set View Zoom( 1 )
	)
);

 Be sure to change the path in the JSL to point to where ever you save the attached Excel file

Jim
joaquimgui
Level II

Re: Plot 3D Diagram

Thank you very much. That is exactly what I wanted to do and it works perfectly.

joaquimgui
Level II

Re: Plot 3D Diagram

I have just one more question. Is there a way to make this code run faster, because I'm dealing with a big data and it takes hours to run and sometimes it doens't. 

My notebook is a Asus, Windowns 7, 64 bit, 4GB (RAM) and 197GB memory. 

So I don't think that the problem is with my computer, right?

 

Thanks.

 

txnelson
Super User

Re: Plot 3D Diagram

My first suggestion is to move the Add Rows() up higher into the looping.  Add Rows() is a pretty slow function, so if you could add all of the rows you need for the data table, to right after the initial New Table(), that would make the program much faster.  That would also mean that you would have to add in an overall row counter, to keep track of what row you are working on, instead of using N Rows().

Names Default To Here( 1 );

// Read in Excel file
dt = Open(
	"c:\users\txjim\documents\read in issue\example.xlsx"
	/*"<your path to the excel file>\example.xlsx"*/,
	Worksheets( "example" ),
	Use for all sheets( 1 ),
	Concatenate Worksheets( 0 ),
	Create Concatenation Column( 0 ),
	Worksheet Settings(
		1,
		Has Column Headers( 0 ),
		,
		Number of Rows in Headers( 1 ),
		Headers Start on Row( 1 ),
		Data Starts on Row( 1 ),
		Data Starts on Column( 1 ),
		Data Ends on Row( 0 ),
		Data Ends on Column( 0 ),
		Replicated Spanned Rows( 1 ),
		Replicated Spanned Headers( 0 ),
		Suppress Hidden Rows( 1 ),
		Suppress Hidden Columns( 1 ),
		Suppress Empty Columns( 1 ),
		Treat as Hierarchy( 0 ),
		Multiple Series Stack( 0 ),
		Import Cell Colors( 0 ),
		Limit Column Detect( 0 ),
		Column Separator String( "-" )
	)
);

// Get Y Matrix
Ymat = {};
For( i = 3, i <= N Col( dt ), i++,
	If( Column( i )[1] != "",
		Ymat = Ymat |/ Matrix( Num( Column( i )[1] ) )
	)
);

// Get rid of unrequired cells, and get the X and Z matrices
dt << select rows( [1, 2, 3] );
dt << delete rows;
dt << delete column( 2 );

// Convert all columns to numerics
For( i = 1, i <= N Col( dt ), i++,
	Column( dt, i ) << data type( numeric )
);

// Get X matrix
Xmat = Column( dt, 1 ) << get as matrix;
dt << delete column( 1 );

// Get Z matrix
Zmat = dt << get as matrix;

// Create new data table, and output the data

dtFinal = New Table( "Converted", New Column( "X" ), New Column( "Y" ), New Column( "Z" ) );
dtFinal << add rows( N Rows( Xmat ) * N Rows( Ymat ) * N Rows( Zmat ) * N Cols( Zmat ) );

rowcnt = 0;

For( Xcnt = 1, Xcnt <= N Rows( Xmat ), Xcnt++,
	For( Ycnt = 1, Ycnt <= N Rows( Ymat ), Ycnt++,
		For( Zrowcnt = 1, Zrowcnt <= N Rows( Zmat ), Zrowcnt++,
			For( Zcolcnt = 1, Zcolcnt <= N Cols( Zmat ), Zcolcnt++,
				rowcnt++;
				dtFinal:X[rowcnt] = Xmat[Xcnt];
				dtFinal:Y[rowcnt] = Ymat[Ycnt];
				dtFinal:Z[rowcnt] = Zmat[Zrowcnt, Zcolcnt];
			)
		)
	)
);
dtFinal << Scatterplot 3D(
	Y( :X, :Y, :Z ),
	Frame3D(
		Set Grab Handles( 0 ),
		Set Walls( 1 ),
		Set Rotation( -54, 0, 38 ),
		Set Text Scale( 2 ),
		Set Line Scale( 1 ),
		Set Marker Transparency( 1 ),
		Set Marker Quality( 0.25 ),
		Set Marker Scale( 1 ),
		Set View Perspective( 0.1 ),
		Set View Zoom( 1 )
	)
);
Jim
joaquimgui
Level II

Re: Plot 3D Diagram

Thank you very much for your help, I really appreciated it.

However when I run the script, it shows an error in the line 70 saying that the value in the columm X can not be determined, because the number of rows (0) is invalid.

I attached the Excel data that I am working on. Would you mind to take a look? I am a new user in JSL, that's why I don't have much experience.

Thank you.

 

txnelson
Super User

Re: Plot 3D Diagram

I downloaded your excel file and saved it.  I then pointed the JSL script to the new file.  I also had to change the line in the code

Worksheets( "example" ),

to

Worksheets( "JMP" ),

because you had changed the name of the spreadsheet within the workbook.

I did not run the program to the full 104,284,944 rows that the JSL will produce if run to completion, but I did run it to a lesser number to prove to myself that the code would work.

If the above change is not the issue, the next question I would have, is what version of JMP are you using?

 

Jim